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This  report  has  two  principal  goals:  First,  to  present  a  general  review  of  magnetic  field 
gradiometers  and,  second,  to  provide  new  data  concerning  these  gradiometers,  including  new 
information  about  their  response  to  magnetic  dipole  fields.  A  system  of  nomenclature  is  in¬ 
troduced  that  is  consistent  with  the  mathematical  concept  of  gradient  and  which  provides  a 
basis  for  discussions  of  the  different  functions  of  magnetic  field  gradiometers  and  differential 
magnetometers.  The  distinction  between  component  gradiometers  and  total  field  gradiometers 
is  also  stressed.  An  historical  review  provides  an  opportunity  to  describe  the  different  charac¬ 
teristics  of  the  many  kinds  of  magnetic  field  gradiometers  that  have  been  developed  since 
the  first  report  of  such  a  gradiometer  in  1925:  rotating  induction  loop,  fixed  induction  loop, 
fluxgate,  proton  precession,  optically  pumped,  and  superconducting  gradiometers  are  discussed. 
It  is  pointed  out  how  the  great  sensitivity  of  superconducting  gradiometers,  and  possibly  other 
varieties  of  modern  magnetic  field  gradiometers,  may  invalidate  the  popular  ‘source-free’  a.s- 
sumption  undec  particular  circumstances.  Further,  because  these  high  sensitivities  will  make 
the  gradiometers  more  susceptible  to  the  geomagnetic  field  gradient,  expressions  are  derived 
for  the  components  of  this  gradient  and  some  representative  numerical  values  are  calculated. 
The  response  of  both  component  and  total  field  gradiometers  to  dipole  sources  is  considered 
for  a  number  of  different  source-gradiometer  configurations.  On  a  more  speculative  note,  two 
varieties  of  rotating  component  gradiometers  are  discussed,  with  particular  attention  being 
given  to  their  possibly  unique  characteristics.  The  report  ends  by  recapitulating  the  many 
applications  of  magnetic  field  gradiometers,  particularly  in  such  important  areas  as  medicine, 
energy  production,  and  defense,  and  by  stressing  the  need  for  gradiometer-related  basic  re- 
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Introduction 


You  boil  it  in  sawdust:  you  salt  it  in  glue: 
You  condense  it  with  locusts  and  tape: 
Still  keeping  one  principal  object  in  view 
To  preserve  its  symmetrical  shape. 

The  Hunting  of  the  Snark 
LEWIS  CARROLL 


§1.1  Preamble 

For  many  years  magnetic  held  gradiometers  have  had  potential  application  in  a  surprisingly 
wide  variety  of  activities:  medical  research  and  di^nosis,  prospecting  Tor  oil  and  ot  her  minerals, 
ship  and  submarine  detection,  the  detection  of  small  arms  (hijack  prevention),  arch(H)log,v 
(<letection  of  artifacts),  studies  of  electric  current  flow  in  the  ionosphere,  and  studies  of  ocean 
waves  and  currents.  These  applications  follow  from  the  magnetic  field  gradiometer’s  primary 
function,  which  is  to  measure  the  space  rate  of  change,  or  gradient,  of  a  magnetic  field.  The 
measurements  can  be  used  in  two  rather  different  ways.  First,  given  a  magnetic  field  from  a 
single  time-varying  source,  a  magnetic  held  gradiometer  can  give  information  about  the  source 
that  would  be  unobtainable  with  a  single  magnetometer.  Thus  the  gradiometer  applicat  ion  to 
studies  of  ionospheric  electric  currents  and  to  studies  of  ocean  waves.  Second,  when  two  or 
more  magnetic  fields  are  present,  any  or  all  of  which  may  be  time-varying  (and  time-varying 
magnetic  fields  are  ubiquitous,  even  though  they  are  usually  of  small  amplitude),  a  magnet  ic 
field  gradiometer  can  sometimes  be  used  to  measure  the  characteristics  of  one  of  t  he  fields, 
even  a  very  weak  field,  while  ignoring  the  ot.>ers.  Thus  the  application  to  medicim*  and  to 
source  detection.  Unfortunately,  .ntil  rece-  magnetic  field  gradiometers  have  lieen  either 
too  bulky  or  too  insensitive  (or,  in  le  Crt—ss,  lioth  too  bulky  and  too  insensit  ive)  to  receive 


2 


I.  lilt rcKliKM  ion 


wide  use  in  any  of  their  areas  of  application.  However,  remarkable  advance.s  in  the  technology  of 
magnetometers  over  the  last  two  decades  have  now  made  possible  the  construct  ion  of  nmgiK't  ic 
field  gradiometers  of  truly  unprecedented  sensitivity  that  are  also  reasonably  compact  even 
portable.  These  new  gradiometers  are  not  yet  widely  available,  but  it  is  likely  that  they  will  find 
increasing  use  over  the  next  few  years,  particularly  in  the  areas  I  have  listed  above.  It  is  even 
possible,  in  my  view,  that  the  gradiometers  will  finally  begin  to  achieve  some  of  the  promist' 
they  have  offered  for  so  long  and  that  their  use  will  lead  to  substantial  new  capabilities  and 
increa.ses  of  knowledge  not  only  in  their  obvious  areas  of  application  but  in  other  unexpected 
areas  as  well.  The  purpose  of  this  report  is  to  provide  a  convenient  reference  for  their  use, 
and  to  provide  some  extra  insight  into  their  possible  applications  and  the  problems  that  might 
ari.se  in  these  applications. 

Magnetic  field  gradiometers  are  not  necessarily  complicated  devices,  as  I  will  show  in  the 
next  section.  Complications  can  arise  in  their  use,  however,  due  to  several  different  factx)rs. 
In  my  experience  the  most  important  of  these  factors  are  intrinsic  practical  limitations  of  the 
equipment,  ambiguous  gradiometer  nomenclature,  and  the  wide  variety  of  diflerent  gradiomet<'r 
nnsponses  to  the  fields  produced  by  typical  sources.  1  will  discuss  each  of  these  factors  in  the 
following  chapters.  I  will  also  discuss  a  particular  kind  of  magnetic  field  gradiometer  t  hat  is 
not  well  known  but  which  has  a  most  desirable  characteristic:  its  sensitivity  is  continuously 
variable  over  a  very  wide  range,  at  least  in  principle.  The  references  I  give  do  not  coniprise 
a  complete  list,  but  they  should  provide  a  good  start  to  any  thorough  search  of  the  magnetic 
field  gradiometer  literature,  which  is  not  extensive. 

§1.2  A  Simple  Magnetic  Field  Gradiometer 

My  acquaintance  with  magnetic  field  gradiometers  began  in  1967,  when  I  was  working  with 
Lee  R.  Tepley  at  the  Lockheed  Palo  Alto  Research  Laboratory.  We  were  studying  the  properties 
of  geomagnetic  pulsations  small  fluctuations  of  the  earth’s  magnetic  field  with  frequencies 
less  than  about  5  Hz  (these  are  called  ultra-low  frequencies,  or  ULF)  and  for  this  purpose  Lee 
had  designed  and  built  a  number  of  steel-cored  multi-turn  solenoids  to  sense  the  geomagnetic 
field  changes.  Since  it  was  these  solenoids  that  provided  the  basis  for  the  simple  gradiometer 
that  I  am  about  to  describe,  the  following  practical  details  of  the  solenoids’  construct  ion  and 
electromagnetic  properties  may  be  of  interest. 

For  portability,  the  solenoids  were  each  made  in  three  separable  parts:  the  central  steel 
core  (a  special  high-permeability  steel),  and  two  identical  coil  .sections.  The  core  had  a  lengt  h 
of  1.84  m,  a  diameter  of  3.7  cm,  and  a  mass  of  about  16.6  kgm;  the  coil  sections  each 
contained  a  nominal  9200  turns  of  aluminum  wire  and  they  were  0.80  m  long,  with  internal  and 
external  diameters  of  4.5  cm  and  7.1  cm,  respectively,  and  their  masses  were  about  10.9  kgm. 
Remembering  that  a  kilogram  is  2.2  pounds,  it  can  be  .seen  that  the  complete  solenoid  units 
were  quite  heavy  (I  learned  about  gradiometer  portability  t  ^  hard  way).  The  coil  sections  wc'n* 
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slipped  over  the  steel  core,  without  being  fastened  to  it,  and  they  were  connected  together 
in  series  to  form  the  assembled  solenoid  unit.  As  might  be  expt'cted,  the  inductance  of  the 
solenoids  was  high:  on  the  order  of  10  H  at  1  Hz.  Their  resistance,  on  the  other  hand,  was 
comparatively  low  (about  100  Q)  due  to  the  u.se  of  heavy  gauge  wire. 

The  electrical  signals  from  the  solenoids  were  first  atnplified  using  low-noise  galvanometc'r 
phototube  amplifiers  (voltage  gain  approximately  500,000)  after  which  they  were  filtc'red  and 
further  amplified  before  being  recorded  on  paper  charts  and  magnetic  tape.  In  general,  two 
solenoids  and  their  associated  equipment  were  installed  at  each  HLF  geomagnetic  pulsation 
observatory.  Both  solenoids  were  horizontal,  with  one  oriented  in  a  geomagnetic  Kast-VVesi 
direction  and  the  other  in  a  geomagnetic  North-.South  direction.  When  fully  adjusted  and 
calibrated,  each  solenoid  system  had  an  approximately  flat  frequency  response  over  the  range 
0.3  7.2  Hz,  with  a  rapid  decline  of  response  outside  t  his  range;  their  internal  noise  was  such 
that  a  2  pT  (peak-peak)  sinusoidal  magnetic  field  variation  wit  h  a  frequency  in  t  he  indicated 
pass  band  could  Just  be  distinguished  on  the  chart  records  (in  more  convent  ional  noise  lenns, 
the  measured  system  noise  was  about  0.4  pT/Hz5  throughout  the  frequency  range  0.3  7.2  Hz.) 

Because  of  their  high  sensitivity  to  magnetic  field  changes,  the  solenoid  syst('nis  were 
always  installed  at  locations  that  were  a  large  distance  away  from  roads  railroa<ls.  and  any 
other  sources  of  man-made  magnetic  fields  that  could  be  identified.  However,  even  with  the 
greatc-st  care,  spurious  magnetic  field  fluctuations  would  be  detected  on  occasion.  For  example, 
during  heavy  winds  the  wire  fences  near  one  pulsation  observatory  would  vibrate  and  produce 
weak  high-frequency  (5  7  Hz)  magnetic  field  changes  at  the  .solenoids.  Although  we  ourselves 
were  not  affected,  we  also  learned  that  cattle  can  create  surpri.singly  strong  magnetic  field 
disturbances  in  their  vicinity  on  occasion:  some  ranchers  deliberately  feed  small  magnets  to 
their  cattle,  which  lodge  harmlessly  in  their  stomachs  and  prevent  any  small  pieces  of  wire 
they  might  ingest  from  progressing  further  through  their  digestive  systems. 

During  the  major  outbreak  of  aircraft  hijacking  that  took  place  in  the  latter  half  of  t  he 
1960’s,  we  began  a  .series  of  experiments  to  see  if  our  sensitive  magnetic  field  measurement 
systems  could  be  used  to  help  prevent  hijacking  by  detecting  the  magnetic  fields  from  concealed 
weapons.  The.se  experiments  showed  that  weapons  with  ferromagnetic  components  coiild  be 
detected  at  useful  ranges  on  many  occa.sions,  but  the  sensitivity  of  the  .solenoid  .sy.stems  was 
limited  by  their  simultaneous  response  to  whatever  natural  magnetic  field  fluctuations  wen' 
in  progress.  Worse,  during  magnetic  storms  or  other  large  naturally-occurring  dist  urbances  of 
the  earth’s  magnetic  field,  it  was  difficult,  and  sometimes  impossible,  to  detect  ferromagnet  ic 
objects  at  any  range.  Faced  with  this  situation,  and  being  aware  of  exploratory  work  at  Varian 
A.s.sociates  on  magnetic  field  gradiometers  constructed  from  their  recently-developed  cesium 
and  rubidium  vapor  magnetometers  [e.g.,  Slack  et  ai,  1967:  Staff  Rept,  1967;  Grice,  1968],  we 
decided  to  combine  two  of  our  solenoids  into  a  single  gradiometer  system  and  to  Ic'st  its  ability 
to  detect  localized  sources  of  magnetic  fields. 


Figure  1.1.  Top  and  side  views  of  the  two-solenoid  gradient  .sensor  used  in  t.h('  simple  magiu'iic 
field  gradiometer  described  in  the  text. 

The  gradiome(.er  was  constructed  as  follows.  The  two  solenoids  were  placed  4  m  apart 
and  orient.ed  in  as  close  to  the  same  horizontal  direction  as  po.ssible  (Figure  l.l).  Next,  their 
wire  leads  were  connected  together  in  such  a  way  that  the  electrical  output  of  the  combination 
was  zero.  or.  significantly,  v(>ry  nearly  zero,  whenever  t  he  two  solenoid  sen.sors  were'  exposed  to 
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the  same  magnetic  field  fiuctuations.  A  small  switching  unit  was  constructed  for  t  his  purpose'; 
it  had  a  two-position  switch  that  connee'ted  together  the  two  pairs  of  input  leads  from  tlu' 
solenoids  in  their  two  possible  combinations.  Then,  with  the  same  magnetic  fit'ld  fluctuations 
applied  to  each  solenoid,  one  switch  setting  would  give  an  output  signal  that  was  twice  the 
signal  produced  by  one  solenoid  alone  (switch  in  the  ‘aiding’  position)  and  the  other  switch 
setting  would  give  the  minimum  output  (switch  in  the  ‘bucking’  position).  The  signals  produced 
by  this  solenoid  combination  were  amplified  and  displayed  using  the  same  equipment  as  was 
used  with  the  individual  solenoids. 

Figure  1.2  (bottom  panel)  shows  three  examples  of  the  respon.se  of  the  solenoid  magnet¬ 
ometer-grad  iometer  system.  The  first  of  these  examples,  1,  is  part  of  a  typical  chart  recorder 
trace  taken  during  an  interval  of  weak  geomagnetic  pulsation  activity.  Both  solenoids  were' 
connecUxl  in  the  aiding,  or  magnetometer,  mode  for  this  record.  The  next  example,  2,  shows 
the  .system  response  when  the  solenoids  were  connected  in  the  bucking,  or  gradiorneter,  mode. 
The  pulsations  shown  in  trace  1  were  still  occurring,  but  the  electrical  output  of  the  two- 
solenoid  combination  was  reduced  almost  to  zero.  The  final  example,  3,  shows  the  response 
of  the  gradiorneter  when  a  small  stool  constructed  largely  of  steel  was  carried  past  the  two 
solenoids  at  a  distance  of  6  m  from  the  nearest  solenoid.  The  experimental  configuration  for 
this  last  experiment  is  detailed  in  the  top  panel  of  the  figure.  Two  large-amplitude  impulsive 
signals  were  recorded,  one  during  the  outward  passage  and  the  other  during  the  return,  and 
if  the  recorder  trace  is  examined  carefully  it  can  be  seen  that  there  is  an  interesting  diflen'nce 
in  the  form  of  the  signals:  a  positive  excursion  occurs  first  on  the  outward  pa.s.sage,  and  a 
negative  excursion  occurs  first  on  the  return.  (The  amplitude  difference  is  not  significant  ;  it 
was  caused  most  probably  by  a  small  decrease  in  the  pa(;e  at  which  the  stool  was  carried 
past  the  solenoids  on  the  return  path.)  The  impulsive  signals  occurred  during  comparatively 
short  time  intervals  centered  on  the  times  of  closest  approach  to  the  near  solenoid,  which 
suggests  that  the  response  of  the  two-.solenoid  gradiorneter  system  is  highly  directional.  This 
high  directionality  is  confirmed  by  theory,  as  will  be  shown  in  Chapter  5.  In  fact,  the  solenoid 
gradiorneter  has  very  little  response  to  the  changing  magnetic  field  of  a  dipole  source  whenever 
the  angle  <!>  shown  in  the  figure  (top  panel)  has  a  magnitude  greater  than  2GP.  When  account 
is  taken  of  the  fact  that  the  response  of  the  gradiorneter  also  varies  roughly  with  distance  r  as 
r“'*,  where  r  is  measured  from  the  center  of  the  gradiorneter  to  the  source  of  magnetic  field, 
the  directionality  becomes  obvious. 

None  of  the  data  I  have  shown  give  any  indication  of  the  maximum  range  for  detection  of 
the  steel  object.  However,  it  can  be  seen  that  this  range  will  depend  primarily  on  two  factors: 
(1)  the  amplitude  of  the  impulsive  signals  produced  in  the  output  of  the  gradiorneter  {e..g., 
trace  3),  and  (2)  the  amplitude  of  the  noise  background  (e.g.,  trace  2).  In  other  words,  it  is  a 
signal-to-noise  ratio  that  determines  the  maximum  range  for  detection,  bet  us  now  consider 
how  this  signal-to-noise  ratio  can  be  increased. 

The  amplitude  of  the  signal  will  depend  in  general  on  the  sensitivity  of  the  solenoid  sensors 
and  on  the  distance  between  the  .source  and  the  solenoids,  i.e.,  the  gradiorneter  sen.sor.  The 
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distance  between  the  source  and  this  gradiometer  sensor  can  be  easily  controlled  under  the 
conditions  I  have  described,  but  in  general,  under  realistic  conditions,  it  is  not  a  variable  I  hat 
the  operator  of  the  gradiometer  can  regulate.  An  effort  can  be  made  to  locate  the  gradiometc-r 
sensor  as  close  as  possible  to  the  location  of  potential  sources,  but  there  is  little  else  t  hat  can 
be  done  to  control  the  distance,  particularly  if  the  sources  are  being  moved  to  avoid  detection. 
This  leaves  the  sensitivity  of  the  solenoid  sensors  under  the  control  of  the  gradiometer  designer, 
and  it  is  plain  that  the  sensitivity  should  be  as  high  as  possible  if  the  gradiometer  is  to  b('  used 
to  detect  sources  at  large  distances.  The  noise  in  the  gradiometer  consists  of  a  combination 
of  the  internal  noise  in  the  solenoids  (basically  thermal  noise)  and  the  residual  geomagnc'lie 
noise.  There  is  not  a  lot  that  can  be  done  about  the  int.ernal  noise.  It  can  be  minimized  by 
keeping  the  resistance  of  the  solenoids  low,  but  this  means  that  the  length  of  the  wire  in  the 
solenoids  must  be  kept  as  small  as  possible,  and  this  can  only  be  done  by  reducing  the  number 
of  turns  and  the  sensitivity  of  the  solenoids.  Thick  wire  can  be  used,  but  it  increases  th(' 
size,  weight,  and  cost  of  the  solenoids.  In  practice,  the  resistance  of  the  solenoids  involves  many 
compromises  and  it  can  only  be  varied  within  strict  limits.  This  leaves  reduction  of  t  he  residtial 
geomagnetic  noise  as  a  means  of  controlling  the  gradiometer  noise  background,  and  it  is  Iktc 
that  considerable  progress  is  possible. 

It  was  never  anticipated  that  the  solenoids  I  have  been  describing  would  be  paired  up 
to  form  a  gradiometer  sensor,  and  no  attempt  was  made  during  their  construction  to  make 
them  absolutely  identical.  As  a  result,  it  is  unlikely  that  the  number  of  turns  on  each  solenoid 
was  the  same.  Further,  the  construction  of  the  solenoids  was  not  designed  with  accurate 
alignment  in  mind.  One  particular  weakness  in  this  regard  was  the  lack  of  a  rigid  connect  ion 
between  the  individual  coils  and  the  steel  cores,  which  made  exact  alignment  of  the  coils 
relative  to  the  cores  and  to  each  other  almost  impossible  to  obtain.  These  two  weaknesses  are 
representative  of  two  of  the  three  major  difficulties  encountered  in  the  design  and  constr>iction 
of  any  magnetic  field  gradiometer,  which,  following  Morris  and  Pedersen  [1961],  I  will  refer  to 
as  (1)  electrical  or  magnetic  mismatch  ,  and  (2)  physical  misalignment.  It  is  essential,  if  the 
greatest  possible  gradiometer  sensitivity  is  to  be  achieved,  that  there  be  complete  symmet  ry 
between  the  constituent  sensors:  they  should  be  identical  electrically  and  magnetically,  and.  if 
component  sensors,  they  should  be  capable  of  being  aligned  in  the  exactly  the  same  direct  ion 
with  the  greatest  possible  accuracy.  The  third  major  difficulty  that  is  encountered  hardly 
needs  separate  listing,  but  it  is  the  most  basic  one  in  gradiometer  design;  it  is  (3)  obtaining 
magnetometers  with  sufficient  sensitivity  to  give  a  gradiometer  with  adequate  range.  The 
history  of  magnetic  field  gradiometer  development  is  essentially  a  history  of  how  increasingly 
more  advanced  technology  has  been  brought  to  bear  on  the  problems  caused  by  the  abovt* 
difficulties.  I  will  review  this  history  in  the  following  chapter,  but  first,  in  the  following  section. 
I  will  provide  some  background  information  on  the  magnetometers  that  are  presently  being 
used,  or  which  might  be  used,  in  modern  gradiometer  construction. 


1.2.  A  Simple  Magnetic  Field  Gradiometer 


(a)  GRADIOMETER  CONFIGURATION 
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Figure  1 .2.  Examples  of  the  magnetic  field  data  recorded  with  the  two-solenoid  gradiometer 
sensor  shown  in  Figure  1.1.  The  first  two  traces  in  the  bottom  panel  show  the  magnetometer 
and  gradiometer  responses  of  the  sensor,  and  the  bottom  trace  shows  what  happens  to  the 
gradiometer  trace  when  a  small  steel  object  is  carried  by  the  solenoids,  as  depicted  in  t  he 
upper  panel. 
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§1.3  Technology  and  Modern  Gradiometers 

If  the  gradiometer  technology  described  in  the  previous  section  were  representative  of 
its  ultimate  possible  development,  it  would  be  safe  to  say  that  magnetic  field  gradiometers 
would  be  limited  to  a  few  specialized  applications.  Howev<‘r.  as  I  have  indicated  in  Section 
l.l,  recent  technological  developments  have  created  new  opportunities  for  these  gradioniet<Ts, 
since  it  is  now  possible  to  construct  magnetic  field  gradiometers  that  combine  remarkabl(' 
compactness  with  extraordinary  sensitivity.  The  technological  developments  under  discussion 
involve  magnetometers  the  basis  for  any  magnetic  field  gradiometer. 

This  century  has  established  an  unprecedented  record  for  technological  development ,  atui 
the  changes  in  technology  have  been  so  many,  and  in  some  cases  so  remarkable  for  their 
impact  on  the  public,  that  a  few  of  the  changes  have  gone  largely  unrecognized  except  among 
comparatively  small  groups  of  specialists.  The  changes  in  (he  technology  of  magnet ometc'rs 
provide  one  good  example.  These  instruments  have  been  the  subject  of  scientific  re.search  for 
some  hundreds  of  years,  and  by  the  end  of  the  last  century  they  had  reached  a  high  level  of 
refinement,  particularly  those  that  were  used  in  magnetic  observatories  (sw  C’hapter  XXVI 
of  Chapman  and  Bartels  [1940],  for  an  authoritative  discus.sion  of  the  early  history  of  these 
observatories).  Considering  this  long  history  of  development,  it  may  be  difficult  t«  believe 
that  there  are  now  two  to  three  more  different  kinds  of  magnetometers  in  u.se  than  there 
were  at  the  end  of  last  century.  Nevertheless,  it  is  a  fact  that  magnetometers  ba.sed  on  (he 
Hall  effect  and  Sample,  1980],  fiber  optics  jyartt;  and  Windsor,  1980;  Lewis,  1980], 

proton  precession  [Packard  and  Varian,  1954;  Waters  and  Francis,  1958;  Hartmann,  1972], 
thin  films  [West  et.ai,  1963;  Irons  and  Schwee,  1972],  optical  pumping  [Grivet  and  Malnar, 
1967;  Hartmann,  1972,  Aleksandrov,  1978],  fluxgates  [Fromm,  1952;  Serson,  1957;  Gordon  and 
Brown,  1972],  and  the  the  Josephson  effect  in  superconductivity  [Josephson,  1962;  Webb,  1972; 
Clarke,  1974]  are  all  new,  most  of  them  having  been  developed  only  over  the  last  three  decades. 
Any  pair  of  these  magnetometers  can  be  combined  in  a  spaced  configuration  to  give  a  part  icular 
kind  of  magnetic  field  gradiometer,  with  its  own  individual  characteristics.  This  is  one  rea-son 
for  the  new  opportunities  that  have  arisen  in  gradiometer  research. 

The  new  magnetometers  differ  in  many  respects,  but  in  general  they  provide  new  accuracy, 
sensitivity,  reliability,  and  convenience  of  use  for  magnetic  field  measurements.  However,  it  is 
the  increase  of  sensitivity  that  is  their  most  remarkable  feature.  To  illustrate  the  magnitude 
of  this  increase,  consider  the  change  that  has  taken  place  in  studies  of  the  earth’s  magnetic 
field.  As  pointed  out  by  Grivet  avd  Malnar  [1967],  the  traditional  magnetic  field  unit  in  these 
studies  for  many  years  was  the  gamma  (7),  which  was  es.sentially  the  smallest  measurable 
variation  of  the  field.  In  present  MKS/SI  units,  the  gamma  is  equivalent  to  10“”  Tesla,  i.c.. 
1  nanoTesla  (nT).  Practically  all  of  the  newly-developed  magnetometers  used  for  geomagnet  ic 
measurements  have  sensitivies  exceeding  1  nT,  and  the  sensitivities  of  some  of  them  exceed  t  his 
limit  by  many  orders  of  magnitude.  The  magnetometers  based  on  the  .Josephson  effect ,  which 
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are  commonly  referred  to  as  superconducting  magnetometers  (even  though  only  a  very  small 
part  of  them  is  superconducting),  are  the  most  sensitive  of  all  at  the  present  time:  the  first 
commercial  models  could  easily  measure  magnetic  field  variations  with  amplitudes  in  t  he  range 
0.1  1  picoTesla  (pT;  1  pT  =  10"^  nT)  [Fraser-Smith  and  Buxton,  1975]  and  recent  models 
can  measure  variations  as  small  as  10“^  pT  if  required.  In  other  words,  there  has  lanm  an 
improvement  by  a  factor  of  approximately  a  million  in  the  sensitivity  that  can  be  achieved 
with  magnetometers  over  the  last  few  decades. 

Magnetometers  based  on  proton  precession  and  optical  pumping  measure  the  total  field, 
whereas  the  others  I  have  listed  measure  components  of  the  field.  It  is  of  course  always  po.ssible 
in  principle  to  assemble  an  array  of  the  component  magnetometers,  measure  t  hree  orthogonal 
components  of  the  magnetic  field,  and  compute  the  total  field.  However,  difficulties  aris<' 
with  this  procedure  in  practice  and  it  appears  likely  as  a  result  that  gradiometers  constructed 
from  total  field  magnetometers  will  form  a  special  class  of  gradiometer  for  a  long  time  to 
come,  with  their  own  special  applications.  Of  the  component  magnetometers,  superconduci  ing 
magnetometers  are  the  most  sensitive,  as  I  have  already  pointed  out.  They  have  an  additional 
advantage  in  that  their  sensing  elements  are  .small  plane  superconducting  loops  \e.g.,  Goodman 
et  al.,  1973;  Clarke,  1974],  which  can  be  constructed  with  considerable  precision;  they  can 
also  be  assembled  in  pairs  on  a  single  compact  supporting  structure,  balanced  magnetically 
(with  great  accuracy),  and  connected  in  opposition  to  produce  a  gradiometer  sensing  element  . 
These  elements  are  so  exceptionally  compact  and  integrated  that  .superconducting  gradiometers 
could  well  appear  to  provide  an  exception  to  the  rule  that  magnetometers  are  the  basis  for  any 
magnetic  field  gradiometer.  However,  as  noted,  the  individual  small  plane  superconducting 
loops  can  also  be  used  as  magnetometer  sensing  elements  once  this  fact  is  recognized,  there  can 
be  no  doubt  that  superconducting  gradiometers  conform  to  the  rule.  Because  .sensitivity  and 
symmetry  are  the  two  most  desirable  characteristics  for  the  pair  (or  pairs)  of  sensors  comprising 
a  gradiometer,  and  superconducting  gradiometers  are  distinguished  by  the  sensitivity  and 
symmetry  of  their  sensing  elements,  it  seems  likely  that  superconducting  gradiometers  will 
provide  the  most  innovative  gradient  mea.surements  for  some  years  to  come. 


(liupIcT  2 

The  Gradient  Concept 


§2.1  Introduction 

Suppose  we  wish  to  measure  a  simple  .scalar  quantity,  temperature  T  for  example,  which 
varies  over  a  known  region  of  space  but  which  is  not  varying  with  time.  We  can  start  with 
a  single  measuring  instrument,  in  this  case  a  thermometer,  and  measure  the  temperature  at 
representative  points  throughout  the  region,  thus  obtaining  a  plot  of  temperature  with  po.sition. 
Figure  2.1  shows  the  kind  of  result  that  might  be  obtained  for  a  two-dimensional  distribution 
of  temperature.  If  the  representative  points  are  closely  spaced,  the  picture  obtained  of  (he 
temperature  variation  throughout  the  region  will  be  comparatively  complete  and  the  contours 
that  are  drawn,  as  in  Figure  2.1,  or  other  representations  of  the  temperature  distribution,  will 
be  comparatively  accurate.  On  the  other  hand,  if  the  points  where  measurements  are  made 
are  widely  spaced,  the  picture  obtained  of  the  temperature  distribution  will  be  comparatively 
incomplete  and  the  contours  or  other  representations  of  the  field  will  be  inaccurate. 

A  difficulty  arises  if  the  temperature  distribution  is  varying  with  time  and  only  a  single 
thermometer  is  available  for  measurement.  If  the  time  variation  is  harmonic,  a  plot  similar  to 
that  shown  in  Figure  2.1  can  be  prepared  by  measuring  the  maximum  value  of  the  temperat  ure 
at  each  point  and  then  using  these  maximum  values  as  original  data  for  the  plot.  However, 
if  the  time  variation  is  not  harmonic,  or  contains  a  non-harmonic  component  (or,  in  practice, 
if  it  is  harmonic  over  too  long  a  time  scale),  arrays  of  thermometers  must  be  Uvsed  to  make 
simultaneous  measurements  of  the  temperature.  Plots  similar  to  those  in  Figure  2. 1  can  then 
be  prepared  showing  the  temperature  distribution  at  particular  in.stants  of  tim--*. 

The  temperatures  shown  in  the  figure  comprise  a  scalar  field:  there  is  a  value  of  temperat  ure 
at  each  point  (either  indicated  specifically  or  implied),  but  no  direction  is  as.sociated  with  the 
values.  We  will  mostly  be  concerned  with  vector  fields  in  this  work,  where  each  jwint  in  the 
region  of  interest  has  associated  with  it  both  the  magnitude  and  direction  of  some  quant  ity,  but 
scalar  fields  can  always  be  constructed  from  these  vector  fields  either  by  neglect  ing  t  he  direct  ion 
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Figure  2.1.  Two-dimensional  temperature  distribution  surrounding  a  source  of  heat,  S.  Four 
temperature  gradient  vectors  Gi,  ..,04  are  shown  in  the  figure  and  discussed  in  the  text. 


and  considering  only  the  magnitude  of  the  quantity  or  by  taking  and  using  these  components 
either  singly  or  in  various  combinations.  Similarly,  vector  fields  can  be  constructed  from  .scalar 
fields  by  a  variety  of  techniques.  One  of  these  techniques  is  to  compute  the  gradient  of  the 
scalar  quantity  at  each  point. 

Without  going  into  great  detail,  an  initial  non-mathematical  picture  of  the  meaning  of  the 
gradient  operation  and  of  the  significance  of  the  vector  quantities  that  are  obtained  by  this 
operation  can  be  obtained  by  considering  the  temperature  contours  in  Figure  2.1.  Suppose 
arrows  are  drawn  in  the  direction  of  increasing  temperature  and  perpendicular  to  the  contours 
at  points  A,  B,  C,  and  D.  It  is  immediately  apparent  that  the  arrows  tend  to  point  toward 
the  source  of  heat  S  that  is  producing  the  temperature  distribution.  We  obvio'  .ly  should  not 
expect  this  directivity  at  all  times,  particularly  when  many  different  sources  and  sinks  of  lu'at 
are  present,  but  it  is  an  interesting  feature  of  the  arrows  that  suggests  they  may  have  some 
physical  significance. 


2.2.  Mathematical  Definition  of  Gradient 


1.1 


To  see  this  physical  significance,  consider  the  flow  of  heat  in  our  hypothetit^al  example. 
We  know  from  the  laws  of  thermodynamics  that  heat  will  not  flow  spontaneously  bet  ween  t  wo 
points  at  the  same  temperature.  It  follows  that  the  direction  of  heat  flow  in  the  figure  must 
everywhere  be  perpendicular  to  the  constant  temperature  contours,  just  as  water  on  the  ground 
tends  to  flow  (when  it  is  moving  slowly)  in  directions  perpendicular  to  the  altitude  contours. 
Thus  the  arrows  in  Figure  2.1  and  the  direction  of  heat  flow  share  the  property  being 
perpendicular  to  the  constant  temperature  contours;  they  differ  in  that  their  actual  din'd  ions 
are  antiparallel,  since  heat  flows  toward  lower  temperatures  and  the  arrows  are  directed  loward 
higher  temperatures. 

Suppose  we  now  consider  other  constant  temperature  contours  close  to  those  drawn  in 
the  figure  and  necessarily  differing  from  the  original  contours  only  slightly  in  (emf)erat  ure.  It 
is  readily  seen  that  the  maximum  rate  of  change  of  temperature  with  distance  occurs  along 
the  perpendiculars  to  the  contours.  This  circumstance  suggests  that  we  should  quantify  our 
arrows  by  making  their  magnitudes  (i.e.,  their  lengths)  equal  to  the  maximum  rate  of  change 
of  temperature  with  distance  at  the  points  where  they  are  located.  The  final  result  linking 
heat  flow  with  the  arrows  in  the  figure  is  then  provided  by  experiment:  it  is  found  that  the  rate 
of  heat  flow  between  two  points  at  different  temperatures  in  a  conducting  medium  is  directly 
proportional  to  the  difference  in  temperature  and  inversely  proportional  to  the  distance  between 
the  two  points.  This  means  that  the  heat  flow  at  points  A,  B,  C,  and  D  in  Figure  2.1  takes 
place  in  directions  antiparallel  to  the  arrows  and  that  the  rate  of  heat  flow  is  proportional  to 
the  magnitude  of  the  arrows.  The  physical  significance  of  the  arrows  is  now  clear;  they  relate 
heat  flow  in  the  figure  to  the  temperature  distribution. 

The  vector  quantities  represented  by  the  arrows  are  referred  to  as  the  gradient  of  the 
temperature  at  the  applicable  points.  More  generally,  given  a  scalar  field  such  as  T,  a  vector 
field  denoted  by  grad  T  (or,  equivalently,  by  VT)  can  be  constructed  from  the  scalar  field  by 
determining  the  direction  and  magnitude  of  the  maximum  spatial  rate  of  increase  of  the  scalar 
quantity  at  each  point  and  using  these  directions  and  magnitudes  to  define  the  vector  field. 


§2.2  Mathematical  Definition  of  Gradient 

Given  a  three  dimensional  scalar  field  f{x,y,z),  the  gradient  of  /  is  a  vector  field  defined 

where  i,  j,  and  k  are  unit  vectors  directed  along  the  positive  x,  y,  and  z  axes  of  the  Cartesian 
coordinate  system.  It  is  not  immediately  apparent  that  this  gradient  quantity  is  identical  to 
the  gradient  discussed  in  the  preceding  paragraphs.  However,  the  following  argument  u.sed  by 
Skilling  [1974],  shows  that  the  two  gradients  are  identical. 


Suppose  we  have  an  infinitesimal  length  dl  given  by 


dl  =  dx  i  +  dys  +  dz  k, 


(2.2) 


which  lies  in  the  scalar  field  /.  Taking  the  scalar  product  of  grad  /  and  dl,  we  obtain 


grad  /  •  dl  =  ^dx  +  ^dy  +  ^dz  =  df, 
dx  ay  dz 


where  d/  is  the  total  differential  of  /,  i.e.,  the  change  in  /  over  the  distance  dl.  Now  consider 
the  scalar  product  in  greater  detail.  First,  grad  /  •  dl  =  0  generally  only  when  dl  and  grad  / 
are  at  right  angles.  This  condition  also  implies  that  d/  =  0,  i.e.,  there  is  no  change  in  /  over 
the  distance  dl,  which  will  only  occur  in  general  when  dl  lies  in  a  surface  of  constant  /.  It 
follows  that  grad  /  must  always  be  normal  to  surfaces  of  constant  /.  This  result  also  implies 
that  grad  /'dl  =  d/  is  a  maximum  when  dl  is  similarly  oriented  in  a  direction  perpendicular 
to  the  surfaces  of  constant  /.  Clearly,  the  direction  of  the  gradient  is  also  the  direction  of 
the  greatest  spatial  rate  of  change  of  the  scalar  function  /  and  its  magnitude  is  actually  equal 
to  this  greatest  rate  of  change.  This  magnitude  will  be  written  |grad  f  \  and  it  follows  from 
Equation  2.1  that 


Igrad 


It  is  interesting  to  note  the  growth  in  complexity  that  occurs  when  the  gradient  operation 
is  used  to  convert  a  scalar  field  to  the  corresponding  gradient  vector  field.  The  initial  .scalar 
field  f{x,y,z)  is  completely  defined  by  a  single  value  of  the  scalar  quantity  at  each  fwint, 
whereas  the  corresponding  gradient  field  requires  knowledge  of  three  quantities  at  each  point 
if  it  is  to  be  completely  defined.  If  measurements  are  made  at  a  single  point,  the  scalar  field 
requires  just  one  measurement,  while  at  least  three  measurements  and  perhaps  as  many  as 
six  measurements  must  be  made  to  obtain  the  gradient.  The  larger  number  arises  because  t  he 
variations  of  /  with  distance  in  each  of  the  i,  y,  and  z  directions  can  only  be  determined,  at  least, 
in  principle,  by  measurements  made  at  two  or  more  locations  spaced  a  small  distance  apart  . 
It  is  possible,  of  course,  to  construct  a  single  sensor  containing  two  or  more  measuring  units 
that  will  automatically  provide  measurements  of  the  spatial  rate  of  change  of  /  in  a  part  icular 
direction.  A  combination  of  three  such  sensors  oriented  at  right  angles  to  each  other,  and  three 
measurements  (one  by  each  sensor),  would  then  suffice  to  determine  the  gradient.  Such  an 
instrument  is  very  naturally  called  a  gradiometer.  The  somewhat  confusing  and  inconsistent 
nomenclature  used  for  gradiometers  at  the  present  time  will  be  described  in  the  following 
chapter  and  a  consistent  nomenclature  introduced,  but  an  inconsistency  immediately  becomes 
apparent  in  the  title  to  this  report.  If  a  gradiometer  measures  the  gradient  of  a  scalar  quantity, 
how  can  there  exist,  .such  an  instrument  as  a  magnetic  field  gradiometer? 
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§2.3  Gradient  of  Vector  Quantities 

The  gradient  of  a  vector  quantity,  such  as  the  magnetic  field,  can  be  defined  by  a  simple 
extension  of  the  basic  definition  given  above  for  a  scalar  field.  Suppose  the  vector  quantity  is 
A,  with  components  {Ax,Ay,A2).  Each  of  the  three  components  designates  a  scalar  field,  and 
the  following  three  gradients  can  be  defined: 

j  ^  dAx.  dAx.  ,  dAx^ 
dAy.  dAy.  dAy 
^  dAx,  .  dAx,  ,  dAx^ 

We  can  then  speak  of  the  gradient  of  the  vector  A  on  the  understanding  that  this  gradient 
consists  of  the  three  independent  parts  grad  Az,  grad  Ay,  and  grad  Az.  In  matrix  notation, 
grad  A  can  be  written 

grad  A  =  [grad  Ai  grad  Ay  gradAz], 


where 


‘dAx  dAz' 

dx  dx  dx 

[grad  Ai  grad  Ay  grad  Az]  =  [i  j  k]  ^  ^  ^ 

dAx  dAy  dAz 
■  dz  dz  dz  ■ 


The  3x3  square  matrix  will  be  referred  to  as  the  gradient  matrix  [G].  In  general,  this  matrix 
contains  9  independently  varying  quantities,  each  of  which  must  be  measured  if  a  complete 
measure  of  the  gradient  of  the  vector  quantity  is  to  be  obtained. 

In  the  special  case  of  the  magnetic  field  B,  not  all  of  the  9  terms  in  the  gradient  matrix 
are  independent.  The  restrictions  on  the  terms  are  of  a  fundamental  nature  and  they  may  be 
derived  directly  from  two  of  Maxwell’s  equations.  The  first  restriction  arises  because  of  the 
requirement  that  div  B  =  0,  which  may  be  written 

S-g-S-o 

Thus  the  sum  of  the  three  diagonal  terms  in  the  gradient  matrix  for  the  magnetic  field  must 
be  zero,  or,  given  any  two  of  the  terms,  the  third  is  restricted  to  the  negative  sum  of  the  ot  her 


T 


2.  The  Gradioiil  ('oncept 


IS 

The  se(‘ond  restriction  of  the  terms  in  the  gradient  matrix  only  applies  in  a  ‘source- free’ 
region,  where  another  of  Maxwell’s  equations  can  be  written  curl  B  =  0.  (A  source-free  region 
here  implies  that  both  the  impressed  current  density  J  and  the  displacement  current  density 
ffDIdt  are  zero).  Equating  the  individual  components  of  curl  B  to  zero,  we  obtain 

dy  dz  '  dx  dz  ’  dx  dy  ' 

These  three  relations  imply  that  the  gradient  matrix  in  a  source-free  region  is  symmetric. 
Combining  Equations  2.7  and  2.8,  it  can  be  seen  that  only  5  of  the  9  terms  in  the  gradient  matrix 
for  the  magnetic  field  can  be  truly  independent  in  such  a  region.  Although  there  is  no  difficulty 
visualizing  .source-free  regions  in  theory,  there  may  be  considerable  practical  difficulty  obtaining 
such  regions  for  gradient  measurements,  particularly  if  measurements  of  high  sensitivity  are 
involved.  The  validity  of  the  source-free  assumption  is  examined  in  Appendix  A. 

§2.4  Gradient  of  the  Total  Field 

As  was  noted  in  Section  2.1,  it  is  possible  to  construct  a  scalar  field  from  a  vector  field 
by  neglecting  the  direction  and  considering  only  the  magnitude  of  the  quantity.  For  example, 
given  the  vector  field  A  =  {Ax,  Ay,  Ax)  considered  in  the  previous  section,  we  have  at  each 
point  a  magnitude  |A|  given  by 

\A\  =  [Ax^  +  Ay^  +  Ax^]^.  (2.9) 

In  many  contexts  this  scalar  quantity  |A|  is  referred  to  as  the  total  field  in  distinction  to  t  he 
components  Ax,  Ay,  and  A*  of  the  field. 

Given  the  scalar  field  represented  by  (A|,  we  can  obtain  the  gradient  by  using 


dx 


dy 


dz 


(2.10) 


Then,  since 
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we  have 


or 


I  A  I  I  I  A  I  /A  A  A  \  • 

|A|gr«l|A|  ). 

AfA  ?^+A  ^)i 

+  Qy  By  +  dy 

/  A  A  A  ^^2  \  1 


grad  |A|  =  |A|  ‘  {Ax  grad  Az  +  Ay  grad  Ay  +  A^  grad  .4j). 
In  matrix  form,  this  expression  can  be  written 


(2.1 i: 


grad  I  A|  =  |A|  '  [  *  j 


kl[G] 


(2.12) 


where  [G]  indicates  the  3X3  gradient  matrix  defined  in  the  previous  section. 

It  is  easily  seen  that  the  gradient  of  the  magnitude  of  a  vector  (i.e.,  the  gradient  of  the 
total  field)  is  not  a  simple  quantity.  In  practice,  measurement  of  the  gradient  of  the  total  field 
can  be  quite  involved  if  only  component  sensors  are  available,  or  relatively  simple  if  sensors 
are  available  that  provide  a  measure  of  the  total  field.  Let  us  suppose  that  only  com|)onent 
sensors  are  available.  In  this  case  a  single  measurement  of  the  gradient  of  the  total  field 
involves,  in  general,  nine  separate,  simultaneous  measurements  of  the  nine  spatial  rates  of 
change  of  the  three  field  components  as  well  as  three  separate,  simultaneous  measurements  of 
the  actual  components  themselves,  i.e.,  12  simultaneous  measurements  in  all.  On  the  other 
hand,  if  total  field  sensors  are  available,  only  three  measurements  of  the  spatial  rates  of  change 
d\A.\/dx,  d\A.\/dy,  d\A.\ldz  need  be  made  and  the  gradient  can  then  by  calculated  directly 
using  Equation  2.10.  For  measurements  of  the  magnetic  field,  it  turns  out  that  both  component 
and  total  field  magnetometers  are  available  and  it  would  appear,  at  least  on  the  surface,  as  if  (  he 
gradient  of  the  total  field  could  be  measured  most  easily  with  total  field  sensors,  even  though 
there  would  be  some  simplification  in  the  component  measurements  due  to  the  relationships 
between  some  of  the  terms  involved  (as  discussed  in  the  previous  section).  Unfortunately,  total 
field  magnetometers  are  relatively  insensitive.  There  is  no  ambiguity  for  electric  fields;  all 
electric  field  sensors  in  common  use  measure  components  of  the  electric  field. 
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The  Magnetic  Field  Gradiometer 

Nomenclature 


The  prrsent  nomenclature  for  magnetic  field  gradient  measurements  has  developed  largely 
on  an  ad  hoc  basis  and  it  contains  terms  that  are  not  standard  and  which  are  sometimes 
used  differently  by  different  investigators.  In  addition,  there  is  at  least  one  term  (torning  into 
informal  use,  namely  ‘long-baseline  gradiometer,’  which  does  not  in  fact  usually  describe  a 
gradiometer;  the  related  term,  ‘short-bastdine  gradiometer’  is  not  quite  as  undesirable,  but  it 
contains  a  redundancy.  In  this  section  a  .system  of  standardized  nomenclature  is  introduced, 
which  is  then  used  throughout  the  remainder  of  the  report.  This  might  all  sound  very  dry 
and  pedantic,  but  it  is  a  nec<'s.sary  step  for  us  to  take  if  we  are  to  avoid  a  lot  of  unnecessary 
confusion  wliile  discu.ssing  gradiometers,  and  it  gives  us  an  opportunity  to  consider  t  he  es.sent  iai 
characteristics  of  a  magnetic  field  gradiometer. 


§3.1  Gradient  and  Gradiometers 

One  of  the  first  difficulties  that  arises  in  gradiometer  nomenclature  is  the  use  of  the  tc'rm 
‘gradient  ’  for  different  quantities.  On  the  one  hand  we  have  the  well-established  gradient  of 
vector  analysis,  which  is  intrinsically  a  vector  quantity  with  both  magnitude  and  direction 
(Sect  ion  2.2).  On  the  ot  her  hand,  we  have  the  everyday  concept  of  gradient  as  a  slope,  or  spat  iaI 
rate  of  change,  with  direction  undefined  or  at  brat  specified  only  in  vague,  general  terms.  In 
between  these  two  concepts  of  gradient  we  have  the  gradient  defined  in  the  dictionary  as  t  he 
change  in  value  of  a  quantity  per  unit  distance  in  a  specified  direction.  This  definition  is  close 
to  the  one  used  in  vector  analysis,  but.  it  differs  in  one  important  feature:  the  vector  atialysis 
gradient  has  a  direction  that  is  completely  determined,  whereas  the  dictionary  definition 
appears  to  leave  open  the  po.ssibility  of  choice  of  direction.  The  vector  analysis  and  dictionary 
definitions  can  be  made  to  correspond  by  simple  additions  to  the  latter  to  give  the  following: 

The  gradient  is  the  change  in  value  of  a  scalar  (giantity  per  unit  di.stancc  in  a  s|)ocili('d 
<lirection.  i.c..  the  direction  in  which  the  increase  in  value  of  the  (juantity  per  unit  disiancc' 
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III  this  work  (hi'  Icrm  ‘gradient’  will  he  used  in  its  striet  veetor  analysis  sense,  as  described 
in  the  previous  chapter.  The  gradient  is  therefore  a  vector  quantity  derived  iiasically  from  ji 
scalar  function  of  position  according  to  Equation  2.1.  The  gradient,  of  a  vector  (|uant  ity  is  .dso 
meaningful,  provided  it  is  derived  according  to  the  procedure  described  in  Section  ’I'otal 
field  gradients  are  potentially  troublesome  to  measure  unless  total  field  sen.sors  are  available, 
but  there  is  no  inconsistency  in  the  concept  of  a  total  field  gradient  in  the  sen.se  discus.sed  in 
Section  2.4.  We  can  therefore  use  the  term  ‘gradiometer’  without  mi.sgiving  provided  it  is  used 
solely  to  describe  an  instrument  for  the  measurement  of  the  gradient  of  .scalar  or  vi'ctor  fields, 
with  gradient  being  used  in  its  strict  vector  .sim.si*.  Simple  descriptive  terms  can  then  be  added 
to  indicate  t  he  type  of  gradiometer.  For  example,  we  can  speak  of  temperature  gradiornetcTs, 
magnetic  field  gradiometers,  and  so  on. 

§3.2  Number  of  Axes 

A  second  difficulty  that  arises  with  nomenclature  is  a  moderately  simple'  one  involving  the' 
number  eif  measurements  of  the  components  of  the  graelient  that  a  give'ti  graeliomete'r  can  make 
simultaneously.  The  simplest  gradiometers,  such  as  the  one  de'.scribed  in  (’hapte'r  1,  e-onsist 
basically  of  two  sen.sors  space'el  a  small  distaiu'c  apart  anel  the'y  only  have  the  capability  e>f 
me'asuring  eine'  spatial  rate  of  change  at  any  erne  time.  If  it  is  the  gradient  of  a  se  alar  fie'lel 
that  is  being  nu'asure'd,  the  gradiometer  must  in  general  be  reorienti'd  two  more  linu's  along 
orthogonal  directions  in  order  for  all  three  components  of  the  gradient  to  be  obtained.  This 
pro<*edure  may  not  be  a  problem  if  the  scalar  field  is  stationary,  but  if  it  is  varying  with  tinu' 
it  may  be  inqKxssible  to  obtain  a  meaningful  gradient  measurement.  The  difficulty  with  vector 
fields  is  even  more  acute,  since  at  least  five  and  passibly  as  many  as  eight  measurements  have 
to  la'  mad<'  to  obtain  the  gradient  of  a  magnetic  field,  and  nine  measurements  are  re(|uir('d 
in  general.  It  is  possible  of  course  to  argue  that  the  gradient  measurements  only  havi'  to  la' 
made  along  tlu'  diri'ction  of  maximum  spatial  rate  of  change  of  the  .scalar  quant  ity  or  vector 
compoiK'nt  that  is  of  interest,  and  that  this  will  cut  down  on  the  measurements  reiiuiri'd. 
However,  this  is  not  a  valid  argument,  because  the  direction  of  maximum  spatial  rate  of  change 
can  only  la'  dc'tected  with  a  simple  two-sensor  gradiometer  by  varying  its  orientation  until  a 
maximum  measurement  is  obtained.  If  the  fields  are  time  varying  the  gradiinit  measurement 
obtained  by  t  his  means  will  once  again  probably  not  be  meaningful. 

It  is  imfairtant  t  herefore  to  be  able  to  sja'cify  t  he  capability  of  a  gradiometer  to  nu'asun' 
simultaneou.sly  a  given  number  of  orthogonal  spatial  variations.  The  neces.sary  nonu'iiclature 
is  already  part  ly  in  use:  a  gradiometer  that  can  measure  just  one  spatial  variation  at  a  time 
is  sometime's  rc'ferred  to  as  a  single-axis  gradiometer  and  gradiometers  intended  for  magiu'tic 
field  measurements  are  occasionally  descriln'd  as  five-axis  gradiometers.  I  suggest  that  the 
descriptive  term  ‘ra-axis,’  where  n  is  the  numbe'r  of  spatial  variations  that,  can  be  measured 
simultaneously,  be  appended  to  any  gradiometer  reference,  at  least  initially,  so  that  the 
capability  of  the  gradiometer  is  immediately  apparent.  In  summary,  the  general  gradiometi'r 
nomenclature  I  would  recommend  for  use  at  the  pn'sent  tinu'  has  t  he  form 

n-axi.s  (optional:  nuKlc  of  ofieration)  (kind  of  field)  gradiometer 

with  typical  examples  Iw'ing  single-axis  temjM'rature  gradiometer,  fivi'-axis  magiu't  ic  lii'ld  gradio- 
rneter.  and  five-axis  siqjercoiKlucting  magnetic  field  gradiorni'ti'r. 
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§3.3  G radiometers  or  Differential  Magnetometers? 

Lot  us  now  consider  a  third  difficulty  with  terminology,  namely  the  appropriateness  of  lh(' 
term  ‘long-baseline’  in  the  description  of  a  gradiometer.  (Notice  that  I  have  left  a  plac(>  for 
this  term  in  my  nomenclature  designation  above!)  As  should  be  obvious  by  now,  I  hav<'  Immmi 
recommending  terminology  that  is  faithful  to  the  mathematical  basis  of  vector  analysis,  :ind 
this  is  the  first  instance  of  a  description  that  is  frequently  used  in  a  manner  that  is  inconsistent 
with  this  mathematical  basis.  Before  passing  any  |)o.ssibly  arbitrary  judgement  on  the  term,  I 
will  briefly  outline  how  it  arises. 

Sup|)ose  the  solenoid  sensors  in  my  simple  single-axis  magnetic  field  gradiometer  (Section 
1.2)  were  to  be  moved  further  apart  than  the  4  m  that  was  specified  for  the  gradiometer,  but 
with  no  change  made  in  their  orientation.  It  would  .soon  be  found  t  hat  the  spacing,  or  baseline. 
could  be  made  very  large  indeed  before  the  natural  magnetic  field  fluctuations  measured  In- 
each  solenoid  became  so  different  that  the  signals  from  each  solenoid  would  no  longer  cancel. 
The  measurements  that  have  been  made  so  far  of  the  spatial  rates  of  change  of  fluctuations 
with  frequencies  in  the  range  0.3  7.2  Hz  are  not  yet  adequate  enough  for  me  to  (luote  distances 
for  cancellation  with  confidence,  but  I  am  reasonably  sure  that  substantial  cancellat  ion  would 
be  observed  in  the  solenoid  system  for  spacings  of  up  to  1  10  km.  Now,  suppose  we  wished 
to  count  the  number  of  cars  driven  along  a  particular  stretch  of  road  by  placing  one  of  tlu' 
solenoid  sensors  beside  the  road  and  using  the  large  output  signals  produced  every  t  irnt'  a  car 
went  by  [Johnston  and  Stacey,  1968]  to  operate  an  automatic  counter.  If  we  were  to  .set  up 
stich  a  system  we  would  soon  find  that  it  was  quite  unreliable,  because  every  tinu'  there  was 
a  magnetic  storm,  or  other  large  natural  magnetic  disturbance,  the  signal  produced  by  a  car 
would  either  (1)  be  lost  in  the  noise  produced  by  the  natural  disturbance  or  (2)  it  could  not  be 
distinguished  from  other  comparable  signals  produced  by  the  disturbance.  The  solution  to  t  his 
noise  problem  is  to  place  another  .similarly-oriented  .solenoid  sensor  some  distance  away  from 
t  he  road  sensor,  at  a  lo(!ation  where  the  signals  produced  by  the  cars  have  beconu*  too  small  tf) 
be  measured,  and  to  use  its  output  to  cancel  the  natural  noise  component  in  the  out  put  of  t  he 
monitoring  solenoid.  Then,  ideally,  the  only  output  of  the  two-solenoid  system  would  be  the 
signals  produced  by  the  cars  on  the  road,  and  they  could  be  proces.sed  automatically  without 
concern  for  the  effects  of  natural  magnetic  field  fluctuations. 

Monitoring  systems  analogous  to  the,  one  I  have  Just  described  have  bwm  studied  for  a 
number  of  years  [e.g.,  Boger  and  Bostick,  1968;  Smith  et  ai,  1969;  Fowler  et  ai,  1973;  Ko 
et  ai,  1976;  Gouhau  et  ai,  1978;  Cornwell  and  Hart,  1979;  M^ore,  1979]  and  they  provide 
valuable  information.  It  is  systems  such  as  these  that  are  now  being  informally  referred  to  ii.s 
long-baseline  gradiometers.  However,  as  I  will  now  show,  they  are  not  really  gradiometers  at 
all. 

As  pointed  out  '  i  .  previous  chapter,  an  es.sential  feature  of  gradient  mea.surements  is 
the  determination  ot  .-^patial  rates  of  variation  of  the  form  dffdx,  where  /  is  the  scalar  field 
quant  ity  whose  gradient  is  required  arul  x  is  a  measure  of  distanee.  In  practice,  t  he  gradient 
measurement  involves  a  determination  of  ratios  of  the  form  Af  j Ax.  where  AJ  is  the  change 
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in  /  over  the  small  distance  Ax.  Strictly,  the  gradient  term  under  consideration  is  oblained  in 
the  limit,  since 


Of 

—  =  lull  -X-. 
ax  Ar-^O  Ax 


where  I  assume  the  other  independent  variables  are  kept  constant  as  Ax  is  varied.  An  ex¬ 
traordinarily  c  onscientious  experimenter  would  therefore  choose  a  small  distance  over  which  to 
measure  the  change  in  /;  he  (or  she)  would  then  halve  the  distance  and  measure  the  change  in 
/  again.  If  t  he  limit  is  being  approached,  the  ratio  Af  j Ax  will  lie  very  closely  t  he  same  in  the 
two  cases.  Such  a  procedure  would  be  im|Kxs.sibly  tim<*-consuming  if  many  gradient  measuri'- 
ments  were  required,  and  it  is  not  a  regular  proisHlure  in  any  event.  However,  it  does  provide' 
us  with  a  useful  practical  criterion  for  determining  when  a  pair  of  magnetometers  is  sjiaeed 
correctly  for  gradiometer  operation.  The  basic  criterion  of  course  is  that  the  spacing  bet  wee'ti 
the  magnetometers  mast  be  ‘small,’  but  the  question  arisi's  small  relative  to  what?  We  can 
now  answer  this  question  by  stating  the  following  (wo  practical  criteria  for  gradiometer 
operation: 


The  spacing  between  two  magnetometers  is  suitable  for  the  jiair  (o  be  used  as  a  gradiometer 
if  (1)  the  spacing  is  small  relative  to  the  applicable  gradiometer-lo-source  distance.  If  this 
criterion  leaves  any  doubt,  we  can  further  specify  that  (2)  the  spacing  is  small  when  the  response 
of  the  magnetometer  pair  (when  connected  in  opposition)  to  the  magnetic  field  of  the  source 
varies  linearly  with  the  spacing.  The  latter  criterion  follows  from  Equation  3.1. 

Having  developed  criteria  for  gradiometer  spacing,  we  can  now  establish  whethc'r  a  par¬ 
ticular  magnetometer  pair  is  capable  of  functioning  as  a  gradiometer  under  a  specified  s('(  of 
experimental  conditions.  It  follows  from  the  criteria  that  the  use  of  magnetometer  pairs  for 
monitoring  under  the  specific  conditions  I  have  described  above  for  the  two-.solenoid  .system 
is  not  a  legitimate  gradiometer  application:  the  .spacing  between  the  magnetometers  is  com¬ 
parable  to  the  distance  from  the  .source  to  the  center  of  the  line  joining  the  magnetometers. 
Note,  however,  that  the  magnetometer  pair  could  function  as  a  gradiometer  if  the  source  that 
was  being  monitored  was  moved  far  enough  away  from  the  magnetometers  for  the  gradiomc'ti'r 
criteria  to  apply. 

The  nomenclature  I  would  suggest  for  monitoring  systems  based  on  pairs  of  widely-spaced 
magnetometers  is  once  again  derived  from  descriptions  us«!d  by  some  authors  {e.g.,  Chapman 
and  Nelson  [1957],  and  Breiner  [1973]).  Inst«id  of  Iong-ba.selinc  gradiometer,  1  recommend  u.se 
of  the  term  differential  magnetometer.  Thus  a  magnetic  field  gradiometer  can  be  considered  to 
be  a  particular  kind  of  differential  magnetometer,  i.e.,  one  distinguished  by  the  fact  t.hat  the 
spacing  between  its  sensors  is  small  with  respect  to  the  distance  to  the  sources  whose  gra<iient,s 
are  to  be  measured  [Breiner,  1973).  Short-baseline  gradiometers,  of  course,  are  best,  refi'rred  to 
simply  as  gradiometers  under  moat  circumstances  because  of  the  redundancy  in  terminology, 
but,  there  may  be  .special  .situations,  such  as  the  one  mentioned  at  the  end  of  this  sect  ion,  wlu're 
‘short-baseline’  is  a  meaningful  gradiometer  description. 
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The  reader  is  now  asked  to  refer  iiaek  to  tlie  description  of  tiu'  simple  gradiomete'r  given 
in  Section  1.2,  and  in  particular  to  the  top  panel  of  Figure  1.2,  and  to  question  whether  the 
two-solenoid  system  is  being  used  strictly  as  a  gradiotneter  when  it  detects  the  magnetic  field 
of  the  steel  object.  Correct!  The  system  is  being  used  as  a  differential  magnetometer,  and  not 
as  a  gradiometer,  because  the  distance  from  the  source  of  magnetic  field  to  the  center  of  t  he 
line  Joining  the  solenoids  (8  m)  is  comparable  to  the  length  of  the  ba.seline  (4  m). 

Before  concluding  this  discussion  of  pairs  of  widely-spaced  magnetometers  and  their  triin- 
sition  into  gradiometers  as  their  spacing  is  reduced,  or  as  the  distance  to  the  source  of  magnetic 
fields  is  increased,  I  would  like  to  draw  attention  to  an  important  feature  of  a  magnet  ic  field 
gradiometer ’s  response.  When  the  spacing  between  the  two  magnetometers  is  small  according 
to  the  criteria  I  have  listed,  it  is  still  possible  to  vary  the  spacing  between  zero  and  some'  max¬ 
imum  value  that  will  depend  on  the  particular  source-gradiometer  configuration,  and  there  is 
an  advantage  to  having  the  spacing  close  to  the  maximum  value.  Since  we  are  still  in  the  linear 
regime,  by  definition,  the  difference  signal  produced  by  the  pair  of  magnetometers  is  largest 
when  their  spacing  is  at  its  maximum  value.  A.s.suming  that  the  noi.se  in  the  magn('tometers 
(or  in  the  system  as  a  whole)  is  independent  of  the  spacing,  which  is  likely  to  be  a  gornl  as¬ 
sumption  for  most  systems,  it  can  be  seen  that  the  largest  signal-to-noise  ratio  is  obtained  for 
the  largest  spacing.  There  is  therefore  considerable  incentive  to  have  the  spacing  belwevn  lh<' 
magnetometers  as  large  as  possible  while  still  operating  in  the  gradiomet.er  mode.  For  this  hist 
situation,  and  for  others  where  the  spacing  is  smaller  but  still  comparatively  large.  us('  of  the 
term  ‘long-baseline  gradiometer’  may  well  be  appropriate  to  di.stinguish  the  gradiomeU'r  from 
one  with  a  comparatively  small  baseline  (‘short-baseline  gradiomeU^r’). 


§3.4  The  Order  of  a  Gradiometer 

The  fourth  and  final  difficulty  with  nomenclature  that  1  wish  to  consider  here  arises  in 
the  design  and  application  of  superconducting  gradiometers.  These  gradiometers  are  being 
used  increasingly  for  biomedical  research  and  as  clinical  instruments.  For  these  applications, 
the  gradiometers  are  operated  in  what  Goodman  et  al.  [1973]  call  the  ‘near-to-source  field 
measurement’  mode:  because  of  the  extreme  weakness  of  the  fields,  the  gradiometer  must 
be  located  very  close  to  their  sources  {e.g.,  the  heart,  lungs,  or  brain  of  a  living  person). 
Unfortunately,  the  gradiometers  must  also  be  operated  in  laboratories  or  clinics,  where  sources 
of  electromagnetic  noise  are  abundant  and  located  clo.se  to  the  gradiometers.  Il  is  usually 
impractical  to  construct  a  completely  .shielded  enclosure  for  the  gradiometer  measurements, 
although  at  least  one  such  enclosure  has  been  built  for  research  purposes  [Cohen  et  al.,  1970; 
Cohen,  1975).  A  particularly  effective  solution  to  this  problem  is  to  construct  more  complicated 
coil  configurations  for  the  gradiometer  sensor,  which  make  the  gradiometer  less  sensitive  to  t  he 
magnetic  fields  from  nearby  noise  sources  while  retaining  high  sensitivity  to  the  fields  from 
sources  close  to  the  gradiometer.  I  have  not  yet  described  the  sensors  used  in  superconducting 
magnetometers  and  gradiometers;  they  consist  of  small  superconducting  coils,  usually  of  a 
single  turn,  that  are  connected  to  the  input  coil  of  a  Superconducting  QUant  um  Interfenmc*' 
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Figure  3.1.  Sensing  coils  used  in  superconducting  magnetometers  and  gradiometers.  Tht'  fol¬ 
lowing  configurstions  are  shown:  a)  magnetometer,  b)  simple,  or  first-order,  gradiotiieter,  c) 
ofT-diagonal  gradiometer,  d)  second-order  gradiometer,  e)  asymmetrical  first-order  gradio- 
meter,  f)  asymmetrical  second-order  gradiometer  [Williamson  and  Kaufman,  1981].  All 
of  the  gradiometer  configurations  are  single-axis,  i.e.,  they  only  make  measurements  on  a 
single  magnetic  field  component. 


Device,  or  SQUID  [Clarke  1974]  a  Josephson  junction  device  that  is  the  electronic  heart  of 
any  superconducting  magnetometer  and  gradiometer.  Although  the  sensing  coils  are  small  and 
superconducting,  their  configurations  are  conventional:  a  .single  coil  is  used  as  the  basic  semsing 
element  in  superconducting  magnetometers,  two  coils  connected  in  series  opposition  arc*  used 
as  the  sensing  element  in  superconducting  gradiometers.  Figure  3. 1 ,  taken  from  a  review  by 
Williamson  and  Kaufman  [1981],  shows  these  two  coil  configurations  labelled  a)  and  b).  The 
other  coil  configurations  in  the  figure  represent  most  of  the  varieties  of  the  more  complicated 
configurations  that  are  used  in  practice. 

The  configuration  of  immediate  interest  to  us  is  the  one  labelled  d)  and  referred  to  as 
a  second-order  gradiometer.  It  will  also  be  noted  that  the  basic  (single-axis)  gradiometer 
is  referred  to  by  Williamson  and  Kaufman  as  a  first-order  gradiometer.  We  have  not  so 
far  encountered  these  descriptive  terms  for  gradiometers  in  this  report,  nor  has  t  he  term 
asymmetric  been  used.  In  fact,  this  new  terminology  is  presently  confined  largely  to  thos«> 


3.1.  The  Order  of  h  (Jradiometer 


25 


X  x  +  Ax  x-<-2Ax 

Figure  .3.2.  The  seeond-order  gradiomet<‘r  shown  in  I'igure  3.1  d),  with  an  x-axis  added  and 
labi'lling  to  denote  the  positions  of  the  coils  on  the  axis.  The  quantities  /  shown  abov(' t  h<' 
coils  indicate  what  the  coils  are  measuring. 


superconducting  gradiometers  used  for  biomedical  re.s('arch  and  as  clinical  instruments.  I 
dhscuss  it  here  not  merely  for  completeness  but  for  the  more  positive  reason  that  it  is  u.s<'ful 
terminology  and  the  instrumentation  that  is  b<*ing  develop<*d  could  have  application  in  other 
areas. 

The  two  examples  of  asymmetrical  gradiometers  shown  in  Figure  3.1  ad(>qua(ely  illustrate 
the  meaning  of  the  term  and  we  need  not  Im?  concerned  further  with  it  here.  The  origin  of 
the  descriptions  first-  and  second-order  may  not  be  so  obvious,  unless  it  is  recognized  that  the 
second-order  gradiometer  sensor  shown  in  the  figure  is  simply  two  of  the  first-order,  or  standard, 
gradiometer  sensors  moved  a  small  distance  apart  along  their  common  axis  (in  th<'  examph* 
shown  (his  distance  is  equal  to  the  spacing  of  the  coils)  and  connected  in  series  opposition. 
The  procedure  used  to  construct  a  second-order  gradiometer  is  therefore  wholly  analogous  in 
principle  to  the  procedure  used  to  construct  a  (first-order)  gradiometer  from  a  magnetometer. 
It  should  be  no  surprise  to  learn  that  the  second-order  gradiometer  measures  the  second  spatial 
derivative  of  the  magnetic  field  component  directed  along  the  gradiometer  axis.  The  following 
derivation  shows  how  the  second-derivative  response  can  be  obtained  from  first  princi|>les  using 
only  our  basic  knowledge  of  gradiometers. 

Let  us  assume  that  we  have  two  single-axis  double-loop  gradiometers  of  the  kind  illustrated 
in  Figure  3.1  a),  both  sharing  the  same  axis  and  identical  in  every  way.  Take  this  common  axis 
to  be  the  x-axis,  and  suppose  one  of  the  gradiometers  is  moved  along  the  axis  relative  to  tlu' 
other  gradiometer  until  two  of  the  coils  coincide,  thus  giving  the  configuration  shown  in  Figure 
3.1  d)  once  the  gradiometers  are  connected  in  series  opposition.  The  spacing  between  the  coils 
in  the  gradiometers  is  Ax,  which  is  as.sumed  to  be  small  relative  to  the  distance  to  t  he  source 
of  the  magnetic  field  component  that  they  are  measuring.  For  generality,  and  to  maintain  a 
notation  consistent  with  the  one  used  in  Section  3.3,  I  will  use  /  to  denote  the  quantity  being 
measured.  Looking  more  closely  at  these  measurements,  the  coils  in  the  l-coil  system  (which  I 
am  not  referring  to  as  a  gradiometer  at  this  stage)  make  three  simultaneous  measurements  of 
/.  as  indicated  in  Figure  3.2.  With  the  notation  shown  in  the  figure,  the  gradiorneliT  on  tlu' 
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left  measures  the  quantity 


fix)  -  f{x  +  Ax) 
Ax 

and  the  gradiometer  on  the  right  measures  the  quantity 


fix  +  Ax)  —  /(x  +  2Ax) 

Ax 

The  4-coil  system  measures  the  difference  between  the  measurements  made  by  the  individual 
gradiometers.  Assuming  that  the  coil  spacing  Ax  is  included  in  this  measurement  in  the  same 
way  as  it  is  in  the  measurements  made  by  the  gradiometers,  the  output  of  the  4-coil  sysU'rri 
can  be  written 

1  /  fjx)  -  fjx  Ax)  fjx  +  Ax)  —  fjx  +  2Ax) 

Ax  V  Ax  Ax 

We  now  use  Taylor  series  expansions  to  convert  the  functions  /(x  +  Ax)  and  /(x  -I-  2Ax)  into 
the  following  approximate  forms 


fix  +  Ax)  =  fix)  +  1“  ^  i^x)^^ 

fix  +  2Ax)  =  fix)  (2^^)'- 


The  approximation  involved  in  these  expansions  is  the  neglect  of  higher-order  terms,  atid  as 
Ax  tends  to  zero  the  expansions  become  increasingly  accurate.  In.serting  the  expanded  forms 
of  fix  -t-  ^x)  and  fix  +  2 Ax)  into  the  expression  for  the  output  of  the  4-coil  system  and  taking 
the  limit,  we  obtain 


/fix)- fix  +  Ax)  _  fix  +  Ax)-fix  +  2Ax)\  ^  ^ 

4i— oV  Ax  Ax  )  dx^ 

This  shows  that  the  4-coil  system  does  indeed  measure  the  second  spatial  derivative  of  the 
variable  of  interest.  It  is  therefore  appropriate  to  refer  to  it  as  a  second-order  (single  axis) 
magnetic  field  gradiometer  an  instrument  that  measures  the  second  spatial  derivative  of  the 
magnetic  field  component  directed  along  its  axis. 

The  advantage  of  the  second-order  single-axis  gradiometer  is  its  increased  discrimination 
against  noise  from  nearby  sources.  If  even  greater  discrimination  is  needed,  it  is  po.ssible  to 
repeat  the  gradiometer  generation  process  with  two  second-order  single-axis  gradiometers  and 
thus  produce  a  third-order  single-axis  gradiometer  that  measures  d^f  fdx^  [Vr6a  et  al.  1982). 
The  procedure  could  be  repeated  ad  infinitum,  but  probably  with  a  rapidly  diminishing  return. 
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The  Magnetic  Field  Gradiometer 
Historical  Development 


§4.1  Early  G radiometers 

The  first  published  article  describing  a  magnetic  field  gradiometer  appears  to  be  one 
authored  by  Haalck  [1925],  who  based  his  gradiometer  on  what  was  then  a  well-known  measur¬ 
ing  instrument  in  Physics  laboratories-  the  earth  inductor.  Haalck’s  article  is  written  in 
German,  but  an  English  summary  is  given  by  C.  A.  Heiland  in  a  later  article  by  Roman  and 
Sermon  [1934]  describing  a  similar  form  of  gradiometer.  Both  of  the  gradiometers  described  by 
these  authors  consist  of  two  identical  coils  spaced  a  small  distance  (1-2  m)  apart  and  rotated 
together  at  the  same  speed,  either  on  parallel  shafts  [Roman  and  Sermon,  1934]  or  on  a  single 
common  shaft  [Haalck,  1925].  Haalck’s  gradiometer  was  more  typical  in  that  his  coils  were 
connected  in  series  and  in  opposition,  with  the  pair  producing  no  output  signal  when  the  two 
coils  were  rotated  in  the  same  magnetic  field.  When  the  two  coils  were  not  located  in  the 
same  magnetic  field,  the  output  signal  consisted  of  the  difference  of  the  emf’s  induced  in  the 
coils.  It  is  worth  noting  that  this  output  signal  was  proportional  to  the  speed  of  rotation  of 
the  coils.  Roman  and  Sermon  measured  the  ratio  of  the  emf’s  produced  in  their  coils,  which  is 
independent  of  the  speed  of  rotation  of  the  coils  and  which,  in  addition,  is  not  a  direct  measure 
of  the  gradient  of  the  geomagnetic  field.  (However,  Roman  and  Sermon  point  out  that  the 
ratio  measures  a  relative  gradient  that  is  nearly  proportional  to  the  actual  gradient.) 

The  work  of  Haalck  [1925]  and  Roman  and  Sermon  [1934]  broke  new  ground,  but  it  was 
never  followed  up,  as  one  would  expect,  by  further  publications  describing  improvements  to 
their  rotating  coil  gradiometers.  The  reason  for  this  apparent  dead-end  to  their  line  of  research 
is  not  completely  clear.  However,  the  development  of  fluxgate  magnetometers  and  re.search 
during  World  War  II  on  magnetic  anomaly  detection  (MAD;  the  description  is  commonly  used 
synonymously  with  submarine  detection)  that  led  to  more  sensitive  gradiometers  based  on 
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fluxgate  sensors  appear  to  have  played  a  major  role,  as  indicated  by  the  following  comments 
taken  from  an  article  by  Fromm[  1952]: 

British  MAD  developments  were  such  that  by  early  1941  a  two-coil  gradiometer  was  being 
tested  operationally.  This  gradiometer  consisted  of  coils  about  1  ft  in  diameter,  mounted 
coaxially  in  a  supporting  framework  and  separated  by  about  8  ft.  Because  this  gradiometer 
system  measured  the  space  rate  of  change  of  the  magnetic  gradient,  which  varies  inversely 
as  the  fifth  power  of  the  distance  from  the  submarine,  its  performance  was  poor  and  its 
possibilities  very  limited.  However,  its  development  was  continued  in  the  United  States 
until  late  in  1941,  at  which  time  it  became  apparent  that  the  saturable-core  magnetometer 
was  more  promising.  Development  of  other  magnetometers,  including  the  gradiometer,  was 
therefore  terminated. 

According  to  Fromm,  the  coils  in  these  gradiometers  were  not  rotated.  It  was  int^'nded 
that  they  would  be  flown  in  aircraR  and  that  signals  would  be  induced  in  them  because  of 
their  translational  motion  through  the  geomagnetic  and  any  other  magnetic  fields  that  might 
be  present.  The  lack  of  rotation  of  the  coils  may  have  been  a  convenience,  but  it  meant  that 
the  amplitudes  of  the  signals  induced  in  the  coils  were  proportional  to  the  speed  of  the  aircraft, 
which  was  not  very  great  in  the  early  1940’s  (and  which  probably  can  never  be  very  great  in 
any  case  during  a  search  for  submarines).  The  combination  of  slow  aircraR  speed,  giving  small 
induced  signals  in  the  coils,  and  the  inverse  fiRh  power  decline  with  distance  of  the  response 
of  the  gradiometer  to  submarine  fields,  was  a  severe  handicap  to  overcome,  and  in  retrospect 
it  is  no  surprise  that  the  performance  of  the  gradiometer  was  poor. 

• 

§4.2  The  Fluxgate  Gradiometer 

Fromm  [1952]  indicates  that  gradiometer  development  for  MAD,  and  presumably  for  other 
uses  as  well,  was  abandoned  during  World  War  II  in  favor  of  work  on  the  fluxgate,  or  saturable- 
core,  magnetometer.  Although  research  on  the  fluxgate  principle  had  been  undertaken  during 
the  1930’s,  it  appears  that  V.V.  Vacquier  of  the  Gulf  Research  and  Development  Company  was 
the  first  to  suggest  its  application  in  sensitive  magnetic  field  detectors  and  the  development 
of  fluxgate  magnetometers  for  prospecting  purposes  was  begun  by  Gulf  in  1940.  Hood  [1965] 
reports  that  Vacquier  applied  for  a  (fluxgate)  gradiometer  patent  in  1941,  at  the  same  time  as 
he  applied  for  a  patent  on  the  fluxgate  sensor:  both  patents  were  awarded  in  1946.  As  might 
be  expected,  application  of  fluxgate  magnetometers  to  submarine  detection  followed  rapidly 
on  their  initial  development  for  geophysical  prospecting,  and  Fromm  describes  the  resulting 
instrumentation  in  considerable  detail.  It  appears  that  Vacquier’s  fluxgate  gradiometer  concept 
was  not  pursued  further  during  the  war.  Following  the  war,  an  article  by  Wurm  [1950]  appeared 
describing  wartime  research  in  Germany  on  both  a  fluxgate  instrument  (Feldstarkemesser)  and 
a  gradiometer  based  on  it  (Feldstarkedifferenzmesser).  However,  as  pointed  out  by  Fromm,  the 
instrumentation  that  was  developed  does  not  seem  to  have  been  put  into  operational  use,  and 
the  lack  of  later  publications  concerning  it  suggests  that  the  idea  of  a  fluxgate  gradiometer 
was  not  pursued  further  in  Germany.  In  fact,  by  the  early  1950’s,  when  Wurm  and  Fromm 
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published  their  articles,  the  brief  life  of  ftuxgate  gradiometers  was  rapidly  drawing  to  a  close. 
An  attempt  by  Gulf  to  use  one  for  airborne  prospt'cting  was  a  failure  [Wickerham,  19r)4j, 
although  the  fluxgate  magnetometer  continued  to  be  actively  developed  for  ‘aeromagnetics’ 
[Reford  and  Sumner,  1964].  A  final  paper  describing  a  fluxgate  gradiometer  constructed  for 
laboratory  measurements  was  published  by  Morris  and  Pedersen  [1961].  This  latter  instrument 
had  a  base  length  of  4.6  m  and  a  noise  level  of  about  5  nT/m  peak-to-peak  for  a  bandwidth  that 
was  unspecified  but  which  was  probably  of  the  order  of  0  1  Hz  (i.e.,  the  bandwidth  mentioned 
by  Serson  [1957]  for  his  fluxgates,  on  which  the  gradiometer  was  based). 

One  of  the  driving  forces  behind  the  development  of  magnetic  field  sensors  aher  the  war 
was  the  interest  in  using  the  sensors  for  geophysical  prospecting,  and  fluxgate  gradiometers 
did  not  prove  to  be  useful  h  ^ely  because  of  the  difficulty  ensuring  that  the  sensor  elements 
were  accurately  aligned  at  all  times.  As  pointed  out  at  the  end  of  Section  1 .2,  the  alignment 
problem  is  one  of  the  crucial  problems  that  must  be  faced  when  using  gradiometers  based  on 
component  magnetometers  (in  distinction  to  total  field  magnetometers),  and  it  is  particularly 
difficult  to  ensure  alignment  when  the  gradiometer  is  subject  to  considerable  handling  and 
substantial  vibration,  which  can  hardly  be  avoided  on  an  aircraft.  Wickerham ’s  comments  on 
this  problem  are  still  pertinent  today: 

The  requirements  for  alignment  of  one  element  with  respect  to  the  other  are  impossibly 
high.  If  operated  in  the  most  favorable  position,  namely,  parallel  to  the  total  vector,  the 
total  ambient  field  observed  by  each  element  may  be  of  the  order  of  50,000  gamma.  Then, 
from  an  orientation  standpoint,  a  misalignment  of  1  gamma  represents  cos''  0  =  0.99998, 
or  0  =  1/3°.  Hence  a  transient  misorientation  of  one  element  axis  with  respect  to  the  other 
of  l/^  would  be  recorded  as  a  very  significant  apparent  anomaly. 

Partly  because  of  this  alignment  problem,  there  was  interest  among  the  geophysical  prospecting 
community  in  what  was  called  ‘time-derivative  gradiometers’  [Wickerham,  1954;  Glicken,  1955]. 


§4.3  The  Time-Derivative  Gradiometer  Technique 

The  time-derivative  gradiometer  technique  is  best  illustrated  by  reference  to  Figure  4.1, 
which  is  taken  from  the  article  by  Wickerham  [1954].  Shown  in  the  bottom  portion  of  the 
figure,  beneath  the  line  designating  the  earth’s  surface,  is  an  idealized  intrusion  of  igneous 
rock,  which  is  comparatively  highly  magnetic.  In  this  idealized  example  the  igneous  rock  has 
a  typical  narrow  zone  of  mineralization  at  its  interface  with  the  surrounding  weakly  magnetic 
or  non-magnetic  rock.  Locating  this  mineralization  is  the  primary  goal  of  the  aeromagnetic 
search.  The  upper  part  of  the  figure  shows  an  example  of  the  kind  of  record  of  the  total  field 
that  might  be  obtained  by  an  aircraft  equipped  with  a  total  field  magnetometer  as  it  flies  over 
the  intrusion  at  constant  speed  in  a  horizontal  straight  line.  The  variation  of  the  total  field 
that  is  shown  is  a  variation  with  distance,  but  it  could  as  well  be  considered  a  variation  with 
time  since  the  distance  travelled  by  the  aircraft  is  proportional  to  the  time  measured  from  the 
start  of  the  pass  over  the  intrusion.  The  dashed  line  in  the  figure  shows  the  spatial  rate  of 
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change  of  the  total  field  along  the  line  travelled  by  the  aircraft;  it  is  obtained  by  differentiat  ing 
the  upper  curve  with  respect  to  distance.  With  the  exception  of  one  feature  that  I  will  shortly 
discuss,  the  dashed  curve  is  essentially  the  response  that  would  be  obtained  from  a  single- axis 
total  field  gradiometer  if  it  was  flown  along  the  path  of  the  aircraft  with  its  axis  aligned  along 
the  path.  It  is  for  this  reason,  and  because  the  differentiation  with  respect  to  distance  is 
equivalent  to  differentiation  with  respect  to  time  insofar  as  the  shape  of  the  gradient  curve 
is  concerned,  that  the  system  is  referred  to  as  a  time-derivative  gradiometer.  I  prefer  not  to 
call  it  a  gradiometer,  since  its  sensor  is  a  single  total  field  magnetometer,  but  ‘time-derivative 
gradiometer  technique’  is  descriptive  and  not  inaccurate.  The  defect  of  this  technique  is  not, 
discussed  either  by  Wickerham  or  Glicken,  but  the  reader  will  surely  be  aware  of  it  by  now.  ih(' 
total  field  variation  recorded  by  the  magnetometer  includes  time  variations  of  the  geomagnetic 
field  that  occur  during  the  flight.  As  a  result  of  these  time  variations,  the  small  oscillations  of 
the  gradient  curve  that  are  supposed  to  identify  the  minerali'^ation  could  as  well  be  caused  by 
a  magnetic  storm  or  other  strong  disturbance  of  the  geomagnetic  field. 

The  time  derivative  gradiometer  technique  was  effectively  rendered  obsolete  by  one  of  (he 
two  major  advances  that  have  taken  place  in  gradiometer  technology  since  World  War  II.  In 
this  case  the  major  advance  was  the  invention  and  development  of  compact  and  reasonably 
rugged  total  field  magnetometers — not  fluxgate  sensors  oriented  and  stabilized  by  mechanical 
means  along  the  direction  of  the  magnetic  field  (as  were  used  by  Wickerham  and  Glicken), 
but  magnetometers  with  an  intrinsic  response  to  the  total  field.  These  magnetometers  could 
easily  be  flown  in  pairs  and  their  output  signals  combined  in  opposition  to  produce  a  true 
single-axis  total  field  gradiometer,  thus  avoiding  the  ambiguity  caused  by  natural  geomagnetic 
field  fluctuations  in  the  time-derivative  technique. 


§4.4  Total  Field  Magnetometers  and  Gradiometers 

The  first  of  the  new  total  field  magnetometers  appeared  in  1953  [Packard  and  Varian, 
1954]  and  they  are  now  generally  referred  to  either  as  proton  precession  magnetometers  or, 
more  generally,  as  one  variety  of  a  class  of  resonance  magnetometers  [Grivet  and  Malnar,  1967; 
Hartmann,  1972]  Their  construction  is  simple  and  rugged  and  they  are  therefore  much  used  for 
geophysical  measurements  [Breiner,  1973],  where  these  characteristics  and  their  easy  portability 
are  particular  advantages.  They  usually  do  not  have  a  continuous  output  and  readings  are 
provided  every  few  seconds;  the  actual  rate  depends  to  some  extent  on  the  sensitivity  and 
higher  reading  rates  generally  imply  lower  sensitivity.  Gradiometers  have  been  constructed 
using  these  magnetometers  [Waters  and  Francis,  1958;  Aitken  and  Tite,  1962]  and  commercial 
units  are  available.  Aitken  and  Tite  [1962]  used  a  3.05  m  (10  ft)  baseline  and  achieved  a 
sensitivity  of  about  600  pT/m,  i.e.,  the  smallest  gradient  that  could  be  distinguished  above 
the  system’s  internal  noise  (and  possibly  also  external  noise)  was  about  600  pT/m.  Assuming 
linearity  of  respon.se  (as  did  Aitken  and  Tite),  this  sensitivity  transforms  into  a  value  of  2000 
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Figure  4.1.  Illustrative  profiles  showing  the  variation  with  horizontal  distance  of  the  total  mag¬ 
netic  field  (solid  line)  and  of  the  horizontal  component  of  the  gradient  of  the  total  magnetic 
field  (dashed  line)  associated  with  an  idealized  intrusion  of  igneous  rock  [Wickerham,  1954]. 


pT/m  for  a  1  m  baseline.  The  sensitivity  of  one  recent  commercial  unit  is  advertised  as  being 
0.82  pT/m,  but  the  magnetometer  spacing  required  to  achieve  this  sensitivity  is  a  high  152  m. 
Once  again  assuming  linearity  of  response,  the  sensitivity  transforms  into  125  pT/m  for  a  1  rn 
baseline— a  greatly  improved  sensitivity  compared  with  that  of  Aitken  and  Tite’s  pioneering 
instrument. 

The  second  variety  of  the  new  total  field  magnetometers  made  its  first  appearance  in 
1957.  Operation  of  these  magnetometers  is  based  on  the  principle  of  optical  pumping  of 
electrons  in  a  gas  or  vapor;  the  principle  was  well  known  in  the  1950’s,  but  the  feasibility 
of  the  magnetometers  only  became  apparent  following  the  work  of  Dehmelt  [1957].  Alkali 
metal  vapors  (usually  rubidium  or  cesium  vapors)  or  helium  gas  are  now  most  commonly 
used  as  the  active  medium  in  a  small  absorption  cell,  and  the  magnetometers  are  referred 
to  variously  as  rubidium- vapor,  cesium-vapor,  and  helium  magnetometers  depenain^  an  the 
active  medium.  Because  the  same  general  resonance  theory  applies  to  the  processes  involved  in 
these  magnetometers,  the  optically-pumped  magnetometers  are  sometimes  included  with  the 
proton  precession  magnetometers  under  the  general  title  of  resonance  magnetometers  [Grivet 
and  Malnar,  1967;  Hartmann,  1972].  However,  considering  the  very  substantial  differences 
in  their  construction  and  operation,  it  seems  more  appropriate  to  separate  them,  and  ot  her 
systems  based  on  them,  into  two  classes  as  I  have  done  here. 
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Optically-pumped  magnetometers  have  a  continuous  output,  a  high  signal-to-noise  level, 
and  high  sensitivity,  but  they  are  complicated  instruments  particularly  when  compared  with 
proton  precession  magnetometers.  The  complication  depends  very  much  on  the  intended  use  of 
the  instrument.  If  it  is  to  be  used  in  a  geomagnetic  observatory,  for  example,  where  no  motion  is 
involved,  a  single  sensing  cell  may  be  used  and  the  magnetometer  will  be  comparatively  simple 
[Usher  and  Stuart,  1970].  For  most  field  operations,  however,  and  for  aeromagnetic  surveys 
and  airborne  ASW  in  particular,  additional  sensing  cells  are  required  and  the  complexity  of  t  he 
instrument  is  increased.  There  are  two  reasons  for  having  the  additional  cells,  both  relating 
to  the  dependence  of  the  output  signal  of  an  optically-pumped  cell  on  its  orientation  relative 
to  the  earth’s  magnetic  field  [Reford  and  Sumner,  1964;  Slack  et  al.,  1967).  The  first  reason 
for  adding  cells  is  to  avoid  a  comparatively  minor  error  in  the  measured  value  of  the  field  t  hat 
occurs  when  an  individual  cell  is  rotated  through  I8flP:  a  dual  cell,  single  axis  sensor  .solves  t  his 
problem.  The  other  reason  for  adding  cells  is  to  avoid  dead  zones  in  the  response  of  the  cells. 
These  zones  result  from  the  changes  that  occur  in  the  amplitude  of  the  output  .signal  of  an 
individual  cell  as  its  orientation  changes  with  respect  to  the  earth’s  field.  No  error  is  involved, 
but  the  signal  amplitude  can  drop  below  a  level  required  for  operation  of  the  signal  processing 
electronics.  The  solution  to  this  problem  is  to  triple  the  number  of  cells  and  arrange  for  their 
axes  to  be  mutually  perpendicular.  The  net  result  of  these  two  .solutions  is  a  sensor  assembly 
containing  six  cells.  Despite  the  need  for  a  number  of  cells,  the  resulting  optically  pumped 
magnetometers  can  be  made  rugged  enough  for  airborne  operation,  and  the  large  number 
of  articles  in  the  scientific  literature  describing  measurements  made  by  airborne  instruments 
attests  to  their  reliability  and  sensitivity.  It  might  also  be  added  that  the  standard  airborne 
magnetometer  for  ASW  in  the  U.S.  is  the  Model  ASQ-81  helium  magnetometer,  containing  six 
helium  cells. 

Considerable  use  has  been  made  of  optically  pumped  magnetometers  in  gradiometers, 
usually  for  airborne  use.  In  most  cases  the  two  mj^netometers  comprising  the  gradiometer  are 
flown  one  above  the  other  in  a  predominantly  vertical  configuration,  thus  measuring  the  vertical 
component  of  the  gradient  of  the  total  field  (more  precisely,  because  the  gradiometer  is  being 
flown,  it  will  measure  the  variation  along  some  horizontal  line  of  the  vertical  comfjonent  of  the 
gradient)  [Slocfc  et  al.,  1965).  To  assist  in  the  interpretation  of  the  gradient  data,  theoretical 
studies  have  been  made  of  the  variation  of  the  vertical  component  over  various  representative 
magnetic  anomalies  [e.g.,  Hood  and  McClure,  1965;  Hood,  1965).  Marine  application  of  the 
magnetometers  appears  to  have  been  minimal.  However,  a  cesium-vapor  gradiometer  for  marine 
use  was  constructed  and  tested  in  1966.  Incorporated  into  a  catamaran,  and  intended  to  be 
towed  behind  another  vessel,  the  gradiometer  consisted  of  two  magnetometers  about  4  m  (12 
ft)  apart  and  its  noise  level  in  the  absence  of  swell  was  quoted  as  being  about  4  pT/m  [Staff 
Rept,  1967]. 

The  sensitivities  of  gradiometers  constructed  from  optically  pumped  magnetometers  are 
reasonably  high.  The  sensitivity  of  the  gradiometer  described  by  Slack  et  al.  [1965],  for 
example,  was  about  0.7  pT/m  for  a  30.5  m  (100  ft)  baseline;  the  corresponding  value  for  our 
reference  1  m  baseline  is  20  pT/m.  The  catamaran  gradiometer  had  a  similar  sensitivity  in 
the  absence  of  wave  noise.  It  is  very  likely  that  the  sensitivities  have  been  improved  in  recent 
years. 
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§4.5  Superconducting  Magnetometers  and  Gradiorneters 

We  come  now  to  the  most  recent  major  magnetometer  development:  superconducting  mag¬ 
netometers.  The  development  of  these  magnetometers  followed  the  discovery  of  the  Jo.sephson 
effect  [Josephson,  1962;  Anderson^  1970]  a  discovery  in  basic  science  that  has  since  had  a  major 
impact  on  applied  science  and  engineering.  The  remarkably  sensitive  new  instruments,  includ¬ 
ing  magnetometers,  that  utilize  the  Josephson  effect  are  well  documented  [e.^.,  Webb,  1972; 
Goodman  et  ai,  1973;  Clarke,  1974],  so  I  will  not  review  the  development  of  the  magnetometers 
here.  Instead,  I  will  list  the  following  properties  of  the  magnetometers  that  are  relevant  to 
their  gradiometer  applications;  (1)  their  sensitivities  are  high  compared  with  all  other  mag¬ 
netometers,  and  may  be  several  orders  of  magnitude  better  than  the  sensitivities  of  fluxgat(>  or 
optically-pumped  magnetometers,  (2)  they  measure  components  of  the  magnetic  field,  not  t  he 
total  field  (they  are  therefore  more  clo.sely  comparable  with  fluxgate  magnetometers  than  with 
proton  precession  or  optically  pumped  magnetometers),  (3)  they  have  a  flat  frequency  response 
from  dc  to  frequencies  in  the  kilohertz  range,  (4)  their  sensing  elements  are  small  supercon¬ 
ducting  loops,  as  described  in  the  previous  section,  (5)  their  sensing  loops  can  be  completelv 
isolated  from  the  natural  magnetic  noise  background  by  means  of  superconducting  shields, 
thus  making  possible  direct  measurement  of  their  internal  noise,  and  (6),  a  disadvantage,  their 
sensing  assembly  must  be  maintained  at  liquid  helium  temperatures  while  they  are  in  operation, 
which  means  that  a  liquid  helium  dewar  and  a  regular  supply  of  liquid  helium  are  required. 

The  first  superconducting  gradiometer  of  which  I  have  any  record  is  one  discussed  by  Ottala 
[1969],  but  it  is  unclear  whether  it  was  ever  put  into  operation.  It  appears  therefore  that  the 
first  superconducting  gradiometer  to  be  built  and  tested  is  the  one  described  by  Zimmerman 
and  Frederick  [1971].  This  gradiometer  was  extraordinarily  small  by  conventional  gradiometer 
standards:  the  spacing  between  the  two  sensing  elements  (single-turn  superconducting  loops 
with  a  radius  of  1.25  cm)  was  only  10  cm.  Despite  its  small  size,  Zimmerman  and  Frederick 
inferred  that  its  sensitivity  was  of  the  order  of  1  pT/m  (0.1  pT/m  for  a  standard  1  m 
baseline).  There  has  since  been  much  improvement  both  in  the  internal  noise  characteristics 
of  SQUID  systems  and  in  the  alignment  and  balancing  procedures  used  in  the  construction  of 
multi-loop  gradiometer  elements  [e.g.,  Overweg  and  Walter  Peters,  1978;  Jaworski  and  Crum, 
1980]  and  it  is  probable  that  sensitivities  much  greater  than  0.01  pT/m  (for  a  1  m  baseline) 
will  be  achieved  in  the  near  future.  Indeed,  recent  instruments  have  already  demon.strated 
substantially  improved  sensitivities.  For  example,  Goodman  et  ai  [1973]  measured  a  peak-to- 
peak  magnetic  field  gradient  noise  of  less  than  10 pT/m  (bandwidth  0  40  Hz)  with  a  20  cm 
baseline  single-axis  gradiometer  at  a  remote  test  site.  Wynn  et  al.  [1975]  reported  a  sensitivity 
of  3  X  10“ pT/m  for  a  single-axis  gradiometer  having  sensing  loops  of  2.5  cm  radius  and  a 
30  cm  baseline  (equivalent  to  a  sensitivity  of  9  X  10“'^  pT/m  for  a  1  m  baseline)  and,  more 
recently,  Ketchen  et  al.  [1977]  obtained  a  .sensitivity  of  about  2  pT/m  (at  frequencies  near 
1  Hz,  in  a  1  Hz  bandwidth)  with  an  innovative  thin-film  superconducting  gradiometer  having 
a  ba.selength  of  the  order  of  2.4  cm  (equivalent  to  a  sensitivity  of  5  X  10  pT/m  for  a  I  m 
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baseline).  These  sensitivities  appear  to  be  close  to  the  value  of  3  X  10“^  pT/m  described  as 
‘today’s  target’  in  1971  by  Nicol  [1971].  A  caution  is  needed,  however,  in  interpreting  the 
sensitivities  too  literally.  As  pointed  out  by  Gillespie  et  al.  [1977],  the  sensitivites  reported 
by  Zimmerman  and  Frederick  [1971]  and  by  Wynn  et  al.  [1975]  do  not  take  adequate  account 
of  the  variation  of  the  gradiometer  noise  with  frequency.  Gillespie  et  al.  [1977]  carried  out 
a  careful  investigation  of  the  noise  characteristics  of  their  gradiometer,  which  was  ultimately 
used  to  measure  the  magnetic  field  gradients  associated  with  internal  waves  in  the  sea  [Podney 
and  Sager,  1979],  and  they  found  that  intrinsic  noise  of  its  SQUID  dominated  the  measured 
noise  spectrum  on  quiet  days  for  frequencies  less  than  0.1  Hz.  At  frequencies  above  0.1  11/. 
the  predominant  source  of  noise  appeared  to  be  thermal  noise  from  non-superconducting  metal 
components  located  near  the  SQUID.  Gillespie  et  al.’s  gradiometer  had  sensor  coils  2.225  cm  in 
radius  and  a  baseline  of  25  cm  and  its  measured  sensitivity  near  1  Hz  was  about  0.2  pT/m  in 
a  iHz  bandwidth,  or  about  0.05  pT/m  for  a  standard  1  m  baseline.  It  appears  therefore  that 
further  development  work  is  required  before  Nicol ’s  1971  target  can  be  met.  Nevertheless,  the 
sensitivities  currently  available  with  superconducting  gradiometers  are  remarkable  compared 
with  those  that  were  typical  just  a  decade  previously. 

Superconducting  gradiometers  differ  from  all  other  gradiometers  through  their  combina¬ 
tion  of  high  sensitivity  and  short  baseline.  This  means  that  they  can  not  only  perform  the 
gradiometer’s  traditional  task  of  detecting  and  measuring  the  magnetic  fields  from  distant 
sources  in  the  presence  of  natural  or  man-made  magnetic  noise  but  in  addition  they  can  detect 
and  measure  the  fields  from  sources  close  at  hand  while  rejecting  the  fields  from  distant  sources. 
To  the  best  of  my  knowledge,  this  dual  capability  was  first  explicitly  pointed  out  by  Goodman 
et  al.  [1973],  who  divided  gradient  measurements  into  two  categories:  (1)  away-from-source 
gradient  measurements,  which  are  directed  at  characterizing  the  source,  and  (2)  near-to-source 
measurements,  which  require  the  high  common  mode  noise  immunity  of  the  gradiometer.  In 
each  of  these  categories,  superconducting  gradiometers  are  now  being  developed  that  will 
provide  their  users  with  new  capabilities. 

Because  they  are  so  compact,  a  number  n  of  single-axis  superconducting  gradiometers  can 
be  combined  in  a  single  instrument  to  produce  an  n-axis  gradiometer  that  is  still  moderately 
portable.  If  the  number  n  is  five  or  more,  indicating  that  at  least  five  independent  simultaneous 
measurements  can  be  made  on  orthogonal  spatial  variations  of  a  magnetic  field,  and  if  provision 
is  made  to  measure  the  actual  components  of  the  magnetic  field  as  well  (which  is  easy  to  do  in 
a  superconducting  instrument),  a  new  kind  of  magnetic  field  instrument  is  obtained.  As  was 
pointed  out  in  increasing  detail  by  Nicol  [1971],  Goodman  et  al.  [1973],  and  Wynn  et  al.  [1975], 
such  an  instrument 

can  provide  from  a  single  location  and  a  single  set  of  measurements:  1 )  direction  to  the  source; 

2)  distance  to  the  source;  3)  magnitude  of  the  source  equivalent  dipole;  and  4)  orientation  of 

the  source  equivalent  dipole  [Goodman  et  al.,  1973). 
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The  combination  of  single-axis  gradiometers  and  magnetometers  can  therefore  track  a  moving 
magnetic  dipole  source  while  obtaining  information  about  the  size  and  orientation  of  its  mo¬ 
ment.  This  new  kind  of  instrument  is  sometimes  referred  to  as  a  tracker  and  classifier  (R.  .1. 
Dinger,  personal  communication,  1981).  I  discuss  some  possible  problems  with  the  use  of  thest' 
instruments  in  Appendix  A,  but  there  is  no  doubt  that  their  development  forms  a  unique  new 
chapter  in  the  history  of  magnetic  source  detection. 

The  close-to-source  capability  of  superconducting  gradiometers  has  introduced  an  entirely 
new  application  for  gradiometers:  biomedical  and  clinical  monitoring  of  the  various  magnetic 
fields  associated  with  some  of  the  body’s  important  organs  (primarily  the  heart,  lungs,  and 
brain).  These  biomedical  and  clinical  applications  are  discussed  in  an  excellent  review  by 
Williamson  and  Kaufman  [1981],  to  which  the  reader  is  referred  for  further  details. 


§4.6  Total  Field  versus  Component  Gradiometers 

A  number  of  different  magnetometers  and  gradiometers  are  described  in  the  foregoing 
sections  and  it  might  well  be  asked  why  so  much  effort  has  been  given  to  developing  new 
instrumentation  when  the  same  effort  devoted  to  the  improvement  of  some  of  the  existing 
magnetometers  and  gradiometers  could  have  led  to  a  substantially  enhanced  capability.  It  is 
important  to  recognize,  however,  that  the  apparently  large  number  of  magnetometers,  and  the 
gradiometers  based  on  them,  can  be  divided  into  two  groups  depending  on  whether  it  is  the 
total  field  that  is  being  measured  (proton  precession  and  optically  pumped  magnetometers  and 
gradiometers)  or  a  component  of  the  field  (stationary  and  rotating  induction  loop  gradiometers, 
fluxgate  and  superconducting  magnetometers  and  gradiometers).  Once  this  division  has  been 
made,  it  is  apparent  that  the  number  of  instruments  in  each  group  is  not  large,  but  the  question 
still  remains  of  why  there  should  be  two  groups.  Why,  in  other  words,  has  so  much  effort  gone 
into  developing  both  total  field  and  component  magnetometers  and  gradiometers? 

The  reason  for  the  division  is  an  important  practical  one.  Total  field  magnetometers  are 
comparatively  immune  to  the  effects  of  motion  in  the  earth’s  magnetic  field,  whereas  component 
magnetometers  are  comparatively  strongly  affected.  No  sophisticated  effects  are  involved.  A 
total  field  magnetometer  measures  a  quantity  that  does  not  vary  with  the  orientation  of  the  I 

instrument  (at  least  to  the  first  order;  as  I  have  described  above,  the  measurements  obtained 
with  a  single-cell  optically  pumped  magnetometer  may  vary  to  some  extent  with  orientation. 

However,  this  effect  is  eliminated  in  the  multi-cell  models  intended  for  use  while  in  motion); 

a  component  magnetometer,  on  the  other  hand,  measures  a  quantity  that  varies  substantially  | 

with  orientation.  For  example,  suppose  a  single-axis  component  magnetometer  is  used  to 

measure  the  horizontal  component  of  the  earth’s  field  at  a  location  where  the  component  is 

0.2  X  10““*  T  (0.2  Gauss).  If  for  some  reason  the  magnetometer  is  displaced  through  a  small 

angle  of  0.1°,  the  measurement  will  change  by  0.2  X  10~'’  cos(0.1°)  T,  or  by  30.5  pT.  This  change 

is  easily  measured  by  a  sensitive  component  magnetometer.  Obviously,  it  would  be  irnpo-ssible  | 
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to  interpret  the  measurements  of  a  component  m^netometer  in  motion,  since  the  data  would 
consist  almost  entirely  of  motion-induced  noise.  This  difference  in  motional  sensitivity  is  the 
basic  reason  for  the  effort  to  develop  both  total  field  and  component  magnetometers. 

The  same  motional  problems  associated  with  component  magnetometers  carry  over  to  the 
gradiometers  based  on  them.  If  the  two  magnetometers  comprising  single-axis  gradiometer 
could  be  completely  fixed  in  position  relative  to  each  other  there  would  be  no  problem,  but 
there  is  inevitably  some  elastic  yielding  in  the  supports  and  motional  noise  results.  This 
is  the  reason  for  the  failure  of  the  airborne  fluxgate  gradiometer  described  by  Wickerham 
[1954].  Because  of  these  problems  with  motionally-induced  noise,  component  gradiometers 
have  been  largely  confined  to  uses  where  they  can  be  kept  at  rest,  and  total  field  gradiometers 
used  for  measurements  where  motion  is  involved.  This  situation  may  change  in  the  future 
as  superconducting  gradiometers  come  into  more  widespread  use.  Because  their  coils  can  be 
made  of  thin  films  deposited  on  a  compact  and  rigid  substrate  \e.g.,  Goodman  et  al.,  1973; 
Ketchen  et  ai,  1977],  the  superconducting  magnetometer  sensor  elements  may  prove  resistant 
to  motion-induced  noise.  Unfortunately,  there  are  likely  to  be  other  motion-induced  noises, 
which  will  have  to  be  strictly  controlled  if  the  gradiometer  is  to  produce  useful  data.  Motion  of 
the  liquid  helium  in  the  dewar,  or  of  any  of  the  other  components  of  the  gradiometer,  including 
in  particular  the  superconducting  plates  or  other  devices  used  to  balance  the  gradiometer, 
could  all  lead  to  an  unacceptable  level  of  motion-induced  noise. 

§4.7  The  Future 

To  end  this  record  of  the  historical  development  of  magnetic  field  gradiometers,  the 
following  news  item  points  to  a  possible  new  application  of  gradiometers: 

Use  of  magnetic  gradiometry  technology  for  cruise  missile  terminal  guidance  is  being  evalu¬ 
ated  by  Sperry  Corp.  for  Defense  Advanced  Research  Projects  Agency  (DARPA).  The  study 
is  focused  on  the  types  of  fixed  and  mobile  military  targets  that  could  be  detected  by  this 
passive  sensing  technique  which  is  embodied  in  the  magnetic  anomaly  detectors  used  in 
antisubmarine  warfare  [Staff  Rept.,  1982]. 


There  seems  little  doubt  that  further  new  applications  will  be  forthcoming  in  the  near  future. 
However,  even  if  these  new  applications  do  not  eventuate,  the  present  applications  of  magnetic 
field  gradiometers  would  be  sufficient  to  keep  the  field  active  for  many  years  to  come. 


The  Magnetic  Field  Gradiometer 
Response  to  Magnetic  Dipole  Fields 


§5.1  Relevance  of  Dipole  Fields 

Most  inagiiotic  field  gradiometers  are  ii.>^ed  to  detect  and  then  to  incaaiire  ( he  clinracteriKt  ics 
of  magnetic  di()ole  fields.  This  concentration  of  effort  is  tlu'  result  of  two  factors.  I'irsl. 
gradionK'ti'rs  are  us<'d  most  frequently  to  detect  and  nu'asure  t  he  fields  from  distant  sourc<‘s: 
ships,  submarines,  ore  iKxlies.  buried  archeological  art  ifacts,  and  so  on.  Second,  th('  magm't  ic 
fields  outside  an  arbitrary,  bounded  source  approximate  more  and  more  closely  to  magnetic 
dipole  fi('lds  as  the  distance  from  the  source  is  increased.  It  is  the  combination  of  thes(>  two 
unrelatc'd  factors  that  makes  dipole  fields  of  particular  imj)ortance  in  magnetic  field  gradiometer 
st  udies. 

I  will  not  include  a  mathematical  proof  of  the  Unidency  for  magnetic-  fields  to  rc'sembli' 
dipole  fields  with  increasing  distance,  since  a  complete'  proof  is  given  by  Stratton  [1910]. 
However,  I  will  point  out  that  the  key  to  the  proof  is  the  demonstration  that  the  field  from 
the  arbitrary  .source  can  be  broken  up  into  a  series  of  terms  that  vary  with  distance'  r  as  r  ", 

where  n  takes  values  n  —  3,4,5 .  The  first  of  these  te'rms  (n  =  3)  is  the  dipoh'  tc'rm. 

the  next  term  (n  =  4)  is  a  quadrupole  tx'rm,  the  following  tertn  (n  =  5)  is  an  oe  tupolc'  tc'rrn, 
and  .so  on  to  higher  multipole  terms.  Theis,  at  distances  that  are  moderately  large  compared 
with  the  dimensions  of  the  source,  but  not.  very  large,  one  might  c'xpc'ct  to  sc'c'  a  prc'dominant 
dipole  fic'lel  together  with  some  admixture  of  a  quadrupole  field.  Closer  to  the  source,  furthc'r 
higher-order  terms  will  appear,  but  it  can  be  anticipated  that  thc'sc'  highe'r-orch'r  terms  will 
generally  be  weak  in  comparison  with  the  dipole  and,  to  a  l('s.ser  extent,  the  c|uadruiK>le'  terms. 
The  apfH'aranc'e  of  the  higher-order  non-dipole  terms  as  the  sourcc'-ree-eiver  distance  is  rc'cluced 
may  be  important  under  some  circumstances,  particularly  if  the  field  measurements  arc'  being 
prtx'c'sseci  to  derive  information  about  the  source  and  it  is  assumed  that  the'  lic'icl  is  purc'ly 
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di|>olar.  However,  it  appears  that  the  dipole  approximation  holds  up  very  well  in  practice,  as 
illustrated  by  the  Tollowing  additional  quote  from  Fromm  [1952]; 

TIk'  modern  s<>a-going  submarine  is  large  <>nough  so  (hat  its  various  sec  tions  may  have 
individual  magnetic  characteristics.  However,  at  •listanci's  greater  than  approximately  the 
length  of  the  submarine,  the  effect  of  each  individual  magnetized  (lortion  is  small  relative  to 
that  of  the  overall  magnetic  moment  so  that  for  all  normal  considc'rations  of  its  magnetic 
field  the  submarine  may  be  considered  as  a  single  magnet  ic  dipole. 

The  submarines  referred  to  by  Fromm  are  no  longer  modern,  of  course,  and  the  distance 
criterion  given  in  the  quote  may  no  longer  be  wholly  appropriate  for  the  larger  and  more  com¬ 
plex  submarines  now  entering  service.  Further,  the  increase  in  the  sensitivity  of  magnetometers 
over  the  last  two  decades  is  likely  to  have  made  the  distance  approximation  less  appropriate 
now  than  it  was  in  1952.  However,  once  the  distance  criterion  has  been  changed  suitably 
(distances  greater  than  approximately  twice  the  length  of  the  submarine?),  Fromm’s  (juote  is 
neces.sarily  as  applicable  today  as  it  was  in  1952,  since  it  is  firmly  based  on  physical  principle. 
Having  ascertained  the  relevance  of  magnetic  dipole  fields,  let  us  now  consider  how  magnet  ic 
field  gradiometers  respond  to  these  fields. 


§5.2  Response  of  Single-axis  Component  Gradiometers 

In  this  section  we  will  consider  the  various  responses  of  the  two  major  types  of  singh'-axis 
component  gradiometers  to  the  field  of  a  magnetic  dipole.  There'  are  a  wide  variety  of  po.s.sible 
dipole/gradiometer  configurations  and  thus  a  large  number  of  different  responses  so  many 
different  responses  that  I  will  not  be  able  to  consider  them  all  in  detail.  Fortunately,  there  is 
a  basic  method  for  deriving  the  response  for  a  given  configuration  that  I  will  cover  in  several 
worked  examples.  Thereafter,  I  will  merely  quote  results  and  pre^sent  curvra  showing  particular 
responses.  This  is  the  first  time  this  material  has  appeared  in  a  report  and  there  are  few  other 
published  results  with  which  it  may  be  compared.  The  only  other  similar  derivation  appears  in 
an  article  by  Karp  and  Duret  [1980],  but  the  approach  is  different  from  mine  and  the  resfxjnse 
data  that  are  presented  are  much  more  limited.  However,  Karp  and  Duret  [1980]  provide  some 
data  on  the  response  of  second-  and  higher-order  gradiometers,  which  I  will  not  consider  here, 
and  the  reader  is  referred  to  their  article  for  a  useful  general  account  of  gradiometer  response. 

The  two  major  types  of  single-axis  gradiomeU^rs  are  distinguished  by  the  configuration 
of  the  two  component  magnetometers  comprising  their  gradient  sensor  as.sembly.  These  mag¬ 
netometers  may  be  arranged  in  either  a  parallel  or  a  series  configuration:  in  the  parallel,  or 
paraxial,  configuration  the  baseline  of  the  gradiometer  is  perpendicular  to  the  axes  of  the  mag¬ 
netometers,  t.e.,  perpendicular  to  the  common  direction  along  which  the  magnetometers  are 
measuring  (he  magnetic  field,  whereas  in  the  series,  or  coaxial,  configuration  the  baseline  is 
aligned  along  the  shared  axis  of  the  magnetometers.  The  simple  gradiometer  descrilx'd  in 
the  first  chapter  is  an  example  of  paraxial  gradiometer;  it  could  be  converted  to  a  coaxial 


Figiin'  5.1.  Diagram.s  showing  the  coordinate'  systems  used  to  derive  the  dipole  magnet  ic  field 
responses  of  (a),  left  panel,  the  simple  magnetic  field  gradiometer  described  in  .'section  1.2 
(an  example  of  a  paraxial  gradiometer),  and  (b),  right  panel,  a  coaxial  gradiometer.  The 
arrow  at  the  origin  of  each  .system  represents  t  he  dipole  source. 


gradiometer  by  locating  the  two  solenoids  a  short  distance  apart  on  the  same  axis  (or.  putting 
it  another  way,  in  the  parallel  configuration  the  .solenoids  are  located  side-by-side,  wherea.s  in 
the  seric's  configuration  they  are  located  one  behind  the  other). 

I  will  use  a  short  solid  line  to  represent  the  component  magnetometers,  a.s  in  Fignn'  1 .2  or 
thus:  mmmm  .  The  axis  defined  by  the  line  is  the  axis  of  the  magnetometer.  When  sc'en  ‘end-on.' 
as  will  occur  in  .several  of  the  configurations,  the  magnetometer  will  be  repre.sent('d  by  a  small 
.solid  circle,  thus:  •  .  Every  gradiometer  configuration  will  necessarily  involve  two  component 
magnetometers.  However,  a  few  of  the  configurations  will  be  shown  with  what  ajipears  to 
be  only  one  magnetometer.  In  these  cases  the  magnetometers  comprising  the  gradioiiu'ter 
are  displaced  in  the  2-direction,  i.e.,  in  and  out  of  the  page,  and  the  magnetometer  in  front 
obscures  the  one  behind. 


5.2.1  Simple  .Vlagnetic  Field  Gradiometer 

F.el  us  now  begin  our  analysis  of  the  respon.se  of  magnetic  field  gradiorneters  by  taking 
the  simple  gradiometer  discussed  in  Chapter  1.2  as  a  first  ('xample.  Figure  5.1  (a)  shows  the 
dipole/gradiometer  geometry  that  will  be  used;  the  dipole  is  represented  by  an  arrow  at  the 
origin,  fiointing  in  the  j/-direction,  and  the  gradiometer  in  this  case  is  also  shown  oriented 
in  the  w-direction.  A  2  axis  perpendicular  to  the  x,  y  plane  in  the  figure  is  not  di.splayed. 
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but  its  presence  is  assumed.  Because  the  magnetometers  are  spaced  apart  in  the  i-direction, 
while  measuring  the  magnetic  field  components  in  the  y-direction,  the  gradiometer  is  paraxial 
and  it  mea.sures  ABy/Ax,  where  ABy  is  the  difference  in  the  By  measurements  made  by  t  he 
two  magnetometers  and  t  =  Ax  is  the  gradiometer  spacing.  Assuming,  as  I  will  throughout 
this  chapter,  that  the  spacing  is  small  compared  with  the  distance  r  from  the  dipole  to  t  he 
gradiometer.  the  gradiometer  also  measures  dByfdx. 

The  di|)ole  magnetic  field  components  Br  and  B0  can  be  written 

2noM  <’os0 


Be  = 


HqM  sin  0 
Atrr^ 


where  M  is  the  dipole  moment.  Starting  with  these  basic  equations,  we  now  need  to  derive  an 
expression  for  the  y  component  of  the  magnetic  field,  since  it  is  By  that  is  measuretl  by  the 
gradiometer ’s  .sensors.  Following  some  algebra,  we  obtain 

By  =  ^(2co8=*fl-sin2  0),  (5  2) 

where  K  —  hqM lilt  is  a  constant.  In  cartesian  coordinates,  Eejuation  5.2  becomes 


B^  =  K 


-  x^ 
(x2  +  y2)i 


Now.  to  obtain  the  response  of  the  gradiometer  to  the  dipole  field,  we  need  to  derive  dBy/dx 
from  Equation  5.3.  Carrying  out  the  differentiation,  we  obtain 

^By  Sxjx^  -  4y^) 

dx  (x*  +  y2)5  ■ 

Since  y/x  =  tan  4>,  this  equation  can  also  be  written 

(5.5) 

OX  I** 

Finally,  writing  dx  for  x  and  assuming  that  dx  t,  we  have 


where 


f{<t>)  =  cos^  ^(1—4  tan^  0). 
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Equations  5.6  and  5.7  are  the  ones  I  will  u.se  to  obtain  a  graphical  display  of  the  gradiometer 
response,  which  1  can  now  define  to  be  the  difference  in  the  magnetic  field  measurements  mad(' 
by  the  two  magnetometers  comprising  the  gradiometer.  In  the  present  case,  the  ((uantity  zi/fy 
is  the  response.  Other  quantities  could  just  as  well  be  defined  as  the  responst>,  but  in  practic(' 
it  is  the  difference  in  the  magnetic  field  measurements  that  is  most  likely  to  be  displayed  and 
understood  to  be  the  instrument’s  response  by  the  gradiometer  operator. 

To  investigate  the  variation  of  gradiometer  response  with  distance  to  the  source,  let  us 
consider  E(|uation  5.4  in  its  polar  coordinate  form: 


dx 


=  K 


(I  -  Icot"*  0). 


(5.8) 


This  equation  shows  that  the  gradiometer  response  varies  with  distance  r  to  the  .source  as 
r”'*.  This  should  be  compared  with  the  equivalent  distance  factor  in  the  response  of  a 
magnetometer.  Thus,  in  the  absence  of  a  background  of  natural  magnetic  noise  (and  t  herefore 
limited  only  by  magnetometer  internal  noise),  a  gradiometer  is  intrinsically  much  more  limited 
in  range  thari  a  magnetometer.  When  natural  noise  is  present  the  situation  changes,  as 
di.scussed  in  the  previous  chapters,  and  depending  on  the  level  of  the  natural  background  t  he 
range  of  the  gradiometer  can  exceed  that  of  a  magnetometer  despite  its  r~‘*  decline  of  response' 
with  distance. 


Returning  to  Equations  5.6  and  5.7,  let  tis  suppose  that  the  dipole  shown  in  F'igure  5.1  (a) 
begins  to  move  in  the  positive  y  direction.  The  only  (juantity  in  Equation  5.6  that  varies  during 
this  motion  is  the  angular  function  /(d>),  so  we  need  only  conisider  the  variation  of  /(<^)  wit  h 
position  of  the  dipole  along  the  y-axis  to  determine  the  corresponding  variation  in  the  response 
of  the  gradiometer.  Before  looking  at  this  variation,  however,  we  should  consider  the  act  ual 
physical  situation  that  is  being  analyzed.  Comparing  the  dipole/gradiometer  configuration  in 
Figure  5.1  (a)  with  the  experimental  arrangement  shown  in  the  top  panel  of  Figure  1.2,  an 
exact  equivalence  can  be  seen.  Note  the  direction  of  the  magnetic  field  in  the  experimental 
configuration:  it  will  magnetize  steel  bodies  to  give  a  dipole  moment  in  the  correct  direction 
(the  positive  y  direction)  for  the  theory  derived  above  to  apply. 

The  angular  function  /(0),  or  response  function,  is  plotted  in  Figure  5.2.  There  is  a 
moderately  sharp  peak  of  response  centered  on  (^  =  OP:  it  has  a  half-amplitude  width  of  '.i‘2A° 
and  zeros  at  ^  =  ±26.6°.  Outside  of  the  angular  range  — 20P<<^<20P  the  response  is  always 
less  than  a  third  of  the  maximum  value  for  <^  =  OP,  which,  as  .shown  by  Figure  5.1  (a),  defines 
the  point  of  closest  approach  of  the  dipole  to  the  gradiometer. 

The  simple  gradiometer  does  not  have  conventional  component  magnetometers  a.s  s('nsors. 
and  its  response  is  a  modified  version  of  what  I  have  just  de.scribed.  Since  induct  ion  loops  are 
u.sed  as  sensors,  the  rate  of  change  of  the  dipole  field  at  the  sensors  must  also  be  takc'n  into 
account.  The  change  from  an  increasing  to  a  decreasing  field  at  the  time  of  closest  approach 
int  roduces  a  sign  change  into  the  output  of  the  gradiometer,  which  then  consists  of  a  jwsitive 
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Figure  5.2.  Plots  of  the  response  function  /{<(>)  for  the  simple  magnetic  field  gradiometer 

(‘component  gradiometer’)  and  for  an  equivalent  total  magnetic  field  gradiometer. 

or  negative  pulse  while  the  dipole  moves  moves  from  <t>  =  20P  to  <i>  =  (F,  zero  output  at  the 
moment  of  closest  approach,  and  then  a  pulse  of  opposite  sign  to  the  initial  one  as  the  dipole 
moves  from  <f>  =  (f  to  ^  =  — 20P,  as  was  seen  in  Figure  1.2.  Motion  of  the  dipole  i.i  the  positive 
y-direction  is  equivalent  to  motion  of  the  gradiometer  in  the  negative  y-direction.  Once  the 
gradiometer  moves  below  the  z-axis,  which  is  the  same  as  the  dipole  moving  above  the  d,  line 
in  Figure  5.1  (a),  <t>  becomes  negative). 

Suppose  the  dipole  were  to  be  moved  in  the  pc»itive  i-direction  instead  of  in  the  positive  y- 
direction  while  maintaining  its  moment  in  the  y-direction,  i.e.,  it  is  moved  .sideways  in  Figure  5. 1 
(a)  and  not  upwards.  There  is  a  simple  technique  for  deriving  the  response  of  the  gradiometer 
for  this  new  situation.  All  that  is  necessary  is  to  substitute  dy  cot  for  dx  in  Equation  5.0 
for  the  response  and  then  to  incorporate  the  resulting  tan'*  ^  factor  into  the  angular  function. 
Since  dy  remains  constant  during  the  motion  of  the  dipole  along  the  i-axis,  the  substitution 
once  again  places  all  the  variables  into  the  angular  function.  When  the  sub.stitution  is  made, 
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WO  obtain 


where 


ABy 


3/ioMf 

4irdy^ 


=  cos^  <f)  tan'*  <^  ( 1  —  4  tan^  <!>). 


(5.9) 


(5.10) 


This  last  expression  does  not  look  very  different  from  the  one  in  Equation  5.7,  but  the  extra 
tan'*  <j>  factor  and  the  different  range  of  <{>  (now  OP  18QP)  ensure  that  the  variation  of  the  response' 
function  is  quite  different  from  the  previous  one.  The  new  variation  is  shown  in  Figure  5..‘i(c), 
(top  scale). 

Figure  5.2  also  shows  the  variation  of  the  angular  function  for  a  total  magnet  ic  field  grad- 
iometer  operating  under  exactly  the  same  conditions  as  the  simple  induction  loop  gradiometer. 
The  angular  function  is  derived  in  a  later  subsect  ion,  but  its  variation  is  shown  here  to  |)rovide  a 
contrast,  with  the  component  gradiometer  curve.  It  can  be  s('en  that  the  total  field  gradiometer 
is  much  less  directive  than  the  component  gradiometer. 


5.2.2  Coaxial  Ciradiometer 


The  analy.sis  in  the  previous  subsection  dealt  solely  with  a  paraxial  gradiometer.  We  will 
now  derive  the  response  of  the  coaxial  gradiometer  shown  in  Figure  5.1  (b).  because  of  t  heir 
orientation,  the  con\ponent  magnetometers  in  this  new  gradiometer  measure  the  .same  magnet  ic 
field  component  By  as  did  the  gradiometer  we  have  just  analyzed.  However,  Ix'cause  iIk' 
magnetometers  are  now  spaced  apart  in  the  y-direction,  the  gradiometer  mea.sures  ABylAy. 
where  ABy  is  the  difference  in  the  By  measurements  mad  '  by  the  two  magnetometers  and 
i  =  Ay  is  the  gradiometer  spacing.  .Since  Ax  is  assumed  to  be  small,  the  gradiometer  also 
measures  dByjdy, 

Using  the  same  cartesian  coordinate  expres,sion  for  By  that  was  derived  in  the  pn'vious 
subsection  (Equation  5.3),  it  is  easy  to  show  that 


dBy  _  ^  3y{3x^  -  2y'^) 

(x^  +  j/2)5 


(.5.11) 


where  K  =  fioMlin  as  before.  Introducing  <t>,  we  obtain 


dy  I'* 


(•5.12) 


Once  again  writing  d*  for  x  and  assuming  that  dj.  »  i,  we  have 


=  -4^ 


u 
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Figure  5.3.  Further  plots  of  the  response  function  /(0)  for  single-axis  component  gradiometers. 
Note  the  additional  0  scales  at  the  tops  of  the  upper  two  plots. 


which  is  formally  the  same  as  Equation  5.6,  but  where  we  now  have 

/(0)  =  cos^  <(>  tan  ^  (3  —  2  tan*  <l>).  (5.13) 

If  the  dipole  is  now  moved  in  the  positive  y-direction,  the  only  quantity  in  the  res|K)nse 
equation  that  varies  is  the  angular  function  given  by  Equation  5.13.  This  angular  funct  ion  is 
plotted  in  Figure  5.3  (c). 


dy 

• 

•  • 

• 

Figurf  5.4.  Single-axis  gradiometer  configurations  for  t  he  measurement  of  each  of  t  he  indicat  eii 
nine  spatial  variations  of  the  magnetic  field  components  Bx,  By,  and  B^. 


Let  the  dipole  now  be  moved  in  the  positive  x-direction.  Using  the  same  technique  as  was 
used  for  the  simple  gradiometer,  the  following  equations  for  the  response  are  obtained: 


ABy  = 


where 


f{<t>)  =  cos^ <i>  tan'^<^(3  —  2tan^<».  (•'*.14) 

This  last  resjxmse  function  is  plotted  in  Figure  5.3  (d);  it.  has  a  comparatively  large  mitiimiim 
when  (j>  =  9flP,  i.e.,  when  the  dipole  source  is  closest  to  the  gradiometer.  A  minimum  may 
be  iinexpectetl  in  this  case,  but  it  is  consistent  with  the  fact  that  By  is  fxxsitively  directed  at 
the  gradiometer  sensors  and  that  it  is  decreasing  as  y  increases:  dByfdy  is  negative  at  the 
gradiometer. 


5.2.3  Generalized  Two-Dimensional  Approach 

We  have  now  considered  two  different  dipole/component  gradiometer  configurations  atid 
obtained  four  different  response  curves.  All  of  these  curves  are  shown  in  Figure  5..'h  wIu'H' 
the  curve  in  panel  (a)  will  now  be  recognized  as  being  a  redrawn  version  of  the  compotu'iit 
gradiometer  curve  in  Figure  5.2.  I  will  now  generalize  the  approach  I  have  been  using  through 
one  step  to  include  the  nine  gradiometer  configurations  shown  in  Figure  5.4.  These'  configu¬ 
rations  include  the  two  previous  ones  and  like  them  they  are  still  basically  two-dimensional, 
since  the  centers  of  the  gradiometer  sensor  assemblies  all  lie  in  the  x,  y  plane.  Notice,  however, 
t  hat  the  magnetometer  axes  are  perpendicular  to  the  x,  y  plane  for  three  of  the  configurat  ions 
and  in  another  group  of  three  the  magnetometers  are  located  a  very  small  distance  Az/2  above 
and  below  the  x,  y  plane.  (One  of  the  configurations  belongs  to  both  of  these  groups). 

The  reason  for  the  choice  of  these  nine  configurations  is  indicated  by  the  part  ial  derivatives 
adjacent  to  each  of  the  gradiometers  shown  in  the  figure.  Each  of  the  single-axis  gradiomelers 
measures  one  of  the  spatial  variations  in  the  gradient  matrix.  If  all  nin('  of  the  gradiometers 
were  present  together  and  centered  on  the  one  point,  they  would  provide  a  complete  meastire 
of  the  gradient  of  the  magnetic  field  at  the  central  point.  We  know,  of  course,  that  one  of  the 
three  spatial  variations  dBifdx,  dBy/dy,  and  dBzjdz  will  be  dependent  on  the  other  two,  or 
zero,  because  of  the  requirement  that  div  B  =  0,  and  and  in  a  source-free  region  three  of  the 
other  variations  will  also  cease  to  be  independent. 

To  see  how  the  various  single-axis  gradiometers  respond  to  a  magnetic  dipole  field,  I  will 
once  again  assume  that  there  is  a  dipole  located  at  the  origin  of  their  coordinate  s>  stems  and 
directed  in  the  positive  y-direction.  The  responses  and  angular  functions  can  then  be  derived 
as  in  the  two  preceding  .subsections.  I  will  not  consider  dipole  motion  out  of  the  x,  y  plane  at 
this  stage,  since  it  requires  a  three-dimensional  approach. 

For  a  dipole  directed  along  the  y-axis  there  is  no  0-component  of  magnetic  field;  the  By 
component  was  derived  in  Section  5.2.1  (Equation  5.3)  and  the  component  can  be  deriv('d 
similarly.  When  collected  together,  the  components  have  the  following  cartesian  coordinate 
forms 


B.  =  K 


By  =  K 


(x2  -t-  y2)§ 
.  2y*  —  x^ 

(x2  -t-  y2)^ 


(5.15) 


B,=0, 


•■'.v  .  -j- 
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from  which  the  following  six  partial  derivatives  can  be  derived 


dBx 

3/ioM 

dx 

‘iTrx"* 

dBx 

3/ioM 

dy 

47rx'* 

dBy 

dBx 

dx 

dy  ' 

dBy 

^HqM 

dy 

~  inx^  ' 

dBx 

dBx  _ 

dx 

dy 

47 


(5.160) 


The  remaining  three  partial  derivatives  in  the  gradient  matrix,  dBx/dz,  dByfdz,  and  dB^jeJz, 
cannot  be  derived  from  Equations  5.15,  since  the  equations  provide  no  information  about  t  he 
field  values  out  of  the  x,  y  plane.  It  would  be  expected  that  dBxIdz  =  dBy/dz  =  0,  because 
of  the  rotational  symmetry  of  the  dipole  field  about  its  axis,  and  for  the  same  reaaon  it  would 
be  expected  that  dBg/dz  would  have  some  non-zero  value,  but  this  is  as  far  as  we  can  go  wit  h 
the  present  two-dimensional  approach.  I  derive  the  actual  partial  derivatives  in  the  following 
.section,  using  a  three  dimensional  approach,  and  they  have  the  form 


dz  dz 


dBx 

dz 


3/xoM 


e 

CO.S  <l>  tan  0. 


(5.166) 


Equations  5.16  show  that  there  are  only  four  independent  spatial  variations  under  the  rest  rict  ive 
two-dimensional  conditions  that  have  been  assumed.  Four  of  the  spatial  variations  an'  zero, 
and  when  they  are  matched  up  with  their  corresponding  configurations  in  Figure  5.'1  it  is  found 
that  the  only  configuration  among  the  five  in  the  bottom  row  and  right  column  that  has  any 
gradiometer  response  to  the  dipole  magnetic  field  is  the  one  measuring  dByfdz.  The  remaining 
configurations  include  the  simple  (paraxial)  gradiometer  and  the  coaxial  gradiometer  (hat  we 
have  already  studied.  To  help  in  the  description  of  the  responses,  I  will  begin  s|)ecifying  the 
configurations  in  Figure  bA  by  means  of  the  coordinate  axes  involved  in  the  partial  derivatives. 
Thus  the  configuration  giving  dE,  dy,  for  example,  will  be  referred  to  as  the  x,y  configuration, 
and  so  on.  Collecting  the  results  from  the  previous  .subsections,  and  adding  new  results  where 
needed,  the  responses  of  the  five  configurations  with  a  gradiometer  response  to  motions  of  (he 
dipole  source  along  the  x-  and  y-axes  are  as  follows; 


(1)  X,X  Configuration.  For  dipole  motion  along  the  x-axis  the  re.siH)n8<'  is  given  by 


ABx  =  Kyf{<t>)  =  Ky  co.s^0  tan''’<^(tan'^«/»  —  1),  (5.l7n) 


and  lor  dipolo  inotion  along  the  j/-axis  tho  responst?  is  given  by 

ABx  =  Kj-  =  Kx  cos^<^  tan<^(tan'^<^  —  4). 


(5.176) 


For  the  x-axis  niolion  /(</>)  is  plotted  in  b^giire  5.3  (a)  (top  scale),  and  for  the  y-axis  rnol  ion 
f{(^)  is  plotted  in  Figure  5.3  (b)  (bottom  .scale). 

(2)  X,  Y  Configuration.  For  dipole  motion  along  the  x-axis  the  response  is  given  by 


A  fix  —  Ky  /(</))  —  Ky  cos^  <i>  tan^  </>(!  —  4  tan'^  0), 
and  for  dipole  motion  along  the  y-axis  the  response  is  given  by 

ABx  =  Kx  /(</>)  =  Kx  eos^  ^(1—4  tan*^  <))). 


(5.17r) 


(5.17d) 


For  the  x-axis  motion  /(<^)  is  plotted  in  F’igure  5.3  (b)  (top  scale),  and  for  the  y-axis  mot  ion 
/(<A)  is  plotted  in  Figure  5.3  (a)  (bottom  scale). 

(3)  Y,X  Configuration.  Same  respon.s(>  as  for  the  i,y  configuration.  F'or  dipole  mot  ion 
along  the  x-axis  the  response  ABy  is  the  same  as  the  response  ABx  lor  the  x,y  configuration 
(E(|uation  5.17c),  and  for  dipole  motion  along  the  y-axis  the  respon.se  ABy  is  the  same  as  tlu' 
responsT'  ABy  for  the  x,y  configuration  (F^quation  5.l7d). 

(4)  Y,Y  Configuration.  For  dipole  motion  along  the  x-axis  the  response  is  given  l)y 


ABy  =  Ky  f{4>)  =  Ky  cos^  <j>  tan'’  <j>  (3  —  2  tan^  <i>), 
and  for  dipole  motion  along  the  y-axis  the  response  is  given  by 

ABy  =  Kx  f(<t>)  =  Kx  eos^ ^  tan  <>(3  -  2  tan'^  <^). 


(5.17c) 


(•'•)•  1 7/) 


For  the  x-axis  motion  f{<f>)  is  plotted  in  Figure  5.3  (d),  and  for  the  y-axis  motion  f(<i>)  is  plot  t  ed 
in  Figure  5.3  (c). 

(5)  Z,Z  Configuration.  For  dipole  motion  along  the  i-axis  the  resjionse  is  given  by 


ABx  —  Kyf(<t>)  =  Ky  sin’’0, 

and  for  dipole  motion  along  the  y-axis  the  response  is  given  by 

ABz  =  Kx  f{<t>)  =  Kx  <t>  tan  <^. 


(5.l7y) 


(5.176) 


The  angular  funct  ion  /(<^)  for  the  i-axis  motion  is  plotted  in  Figure  5.3  ((>),  and  for  t  in*  y-axis 
motion  f(<f>)  is  plotted  in  Figure  5.3  (f). 

The  Kx  and  Ky  in  the  above  equations  are  constants  given  by  Kx  =  lifioMf/indJ  and 
Ky  =  ^fioM(/4irdy,  where  dx  and  dy  are  the  |)erj)endicular  distance's  of  the'  centers  of  the 
gradiometer  sensor  a.s.s<«rnblies  from  the  x-  and  y-axes.  It  can  lx*  s«'en  that  tlie  six  angular 
variations  in  Figure  5.3  provide  all  the  information  we  need  to  de'scrilx'  the  atigiilar  r<'spon.s«>s 
of  the  configurations  shown  in  Figure  5.4  provided  the  di|H)le  motions  are  restricted  to  the  x- 
and  y-axes. 


5.2.  Rrapoiisc  of  Single-axis  Component  Cradiometers 


5.2.4  Tolal  Field  Gradionieler 

The  response  of  a  total  magnetic  field  gradionieter  under  the  same  two-dimensional  con¬ 
ditions  I  have  been  assuming  can  be  derived  using  the  same  techniques  that  were  used  for 
the  component  gradiometers.  To  illustrate,  1  will  now  derive  the  response  of  a  total  field 
gradiometer  under  the  same  conditions  as  those  a-ssumetl  for  the  simple  paraxial  gradionieter 
(Subsection  5.2.1);  the  configuration  will  be  identical  to  the  one  shown  in  Figure  5.1  (a),  but 
the  two  magnetic  field  sensors  will  now  be  understood  to  Ix'  total  field  magnetometers. 

The  first  step  in  deriving  the  response  is  to  obtain  the  total  field  in  cartesian  coordinates. 
Using  Equations  5.15  for  the  i  and  y  components,  the  total  field  Br  is  obtained  from 


giving 


Bt  — 


K  (iy'*  -I-  5i^t/^  -t- 

(l2 


where  K  =  fioM /Air  as  before.  DilTerentiating,  we  find 

dBr  _  +  5j/^) 

{x^  +  4y'^)^{x^  +  y^)'^ 

Again  using  y  =  x  tan<^,  1  -I-  tan^^  =  cos“^<^,  we  have 

dBr  _  — SK  cos®  (^  ( 1  +  5  tan'^  (]>) 

dx  x-*  (H-4tan2(/.)i 

Thus,  for  dipole  motion  along  the  j/-axis,  the  response  is  given  by 


(.5.  IS) 


(.5.19) 


(5.20) 


(5.21) 


where 


n<i>)  = 


cos®  ( I  -t-  5  tan^  <j>) 
(1-1-4  tan^  (^)i 


(5.22) 


This  last  angular  function  is  the  one  labelled  ‘Total  Field  Gradiometer’  in  Figure  5.2. 

For  dipole  motion  along  the  x-axis,  the  respon.se  of  the  total  field  magnetometer  is  found 
by  substituting  =  cot  d>dy  in  Equation  5.21  and  incor|X)rating  the  resulting  tan^  </>  term 
into  the  angular  function.  We  then  have 


(.5.23) 


so 
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Figure  5.5.  Additional  plot  of  a  response  function  for  a  total  magnetic  field  gradionieter. 


where 


/W  = 


cos*  <f>  tan'*  <^  ( I  +  5  tan^  (t>) 
(1  +  4  tan^^)i 


(5.2i) 


Figure  5.5  shows  the  variation  of  this  final  angular  function  with  angle  <^.  The  data  in  this 
last  figure,  combined  with  the  total  field  data  in  Figure  5.2,  suggest  strongly  that  the  total 
field  gradiometer  is  basically  a  less-directive  instrument  than  a  component  gradiometer.  It  is 
possible,  however,  that  use  could  be  made  of  the  null  in  the  respon.se  (Figure  5.5)  to  obt  ain 
moderate  directionality. 


§5.3  Three-Dimensional  Approach 

Up  to  this  point  I  have  confined  my  derivation  of  gradiometer  responses  to  magnetic  dipole 
fields  to  moderately  simple  dipole/gradiometer  configurations.  This  has  enabled  me  to  obtain 
some  useful  response  data  that,  as  I  will  shortly  di.scuss,  can  be  applied  to  a  number  of  practical 
situations.  Once  the  configurations  are  made  more  general,  the  derivation  of  gradiometer 
respon.se  becomes  difficult  because  of  the  algebraic  complexities  involved  and  the  response 
functions  are  complicated.  I  will  present  the  results  of  a  generalized  two-dimensional  analysis 
in  the  next  chapter,  and  it  will  be  .seen  that  even  in  two  dimensions  the  res{>ons('  functions 
can  become  particularly  complicated.  In  three  dimensions  the  difficulties  involved  in  obtaining 
respon.se  information  are  very  great,  except  for  the  mo.st  simple  configurat  ions  (in  t  heir  simplest 


5.3.  Thm-'Diniensional  Approach 
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Figure  5.6.  Coordinate  sy.stem  for  the  three-dimensional  analysis  of  gradiornet  er  res|)onse. 


forms  these  conhgurations  are  the  two-dimen.sional  ones  we  have  already  analyzed),  and  I  he  ii.s<' 
of  computers  to  obtain  relevant  numerical  data  is  clearly  indicated.  As  a  result,  I  will  confiiu' 
myself  in  this  section  to  a  presentation  of  the  basic  three-dimensional  component  and  spat  ial 
variation  data,  which  is  the  essential  raw  material  for  a  computer  analysis. 


Suppose  our  magnetic  dipole  is  now  located  at  the  center  of  a  three-dimensional  cartesian 
coordinate  system  with  the  dipole  moment  directed  along  the  2-axis,  as  shown  in  Figure  5.6. 
The  three  components  of  the  dipole  magnetic  field  can  be  derived  from  the  two  two-dimensional 
polar  coordinate  components  Br  and  8$  given  in  Subsection  5.2.1  once  account  is  taken  of  the 
fact  that  there  is  no  azimuthal  component  of  magnetic  field  in  Fig»>re  5.6.  The  two  eom{)onen1s 
are  shown,  in  vector  form,  in  the  figure,  and  the  three  corresponding  cartesian  coordinate' 
components  of  the  field  are 

3Kx2 


rS  ’ 


(5.25) 


-  r^) 


where  the  same  K  =  hqM JAir  is  used  as  in  the  previous  section.  Carrying  out  the  appropriate 
differentiation,  the  gradient  matrix  can  be  written 


5S 
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dBr  dBz 

dx  dx  dx 
dBx  dBy  dBz 
dy  dy  dy 

dBr.  dBy  dBz 
-  dz  dz  dz  • 


[yir^  -  5i^)  xir^  -  5y2)  -bxyz 
=  ~  x(r^  -  5y^)  y(3r^  -  hy‘^)  z(r^  -  hy‘^) 
’’  —bxyz  z(r^  —  hy^)  y(r^  -  5z^). 


(5.26) 


The  nine  expressions  in  the  evaluated  gradient  matrix  are  the  raw  material  for  any  three 
dimensional  response  analysis.  It  will  be  noted  that  the  matrix  is  symmetric  and  the  sum  of  the 
three  diagonal  terms  is  zero,  as  is  required  for  the  gradient  of  a  magnetic  field  in  source-free 
region. 

To  illustrate  the  use  of  this  material,  suppose  we  have  a  component  gradiometer  with 
its  center  located  at  the  point  (x,  y,  —\z\),  with  its  two  sensors  both  located  in  the  x,  y 
plane  passing  through  the  center,  and  with  its  sensitivity  axis  parallel  to  the  j/-axis.  This 
configuration  is  deliberately  chosen  for  its  pertinence  to  the  two-dimensional  simple  magnetic 
field  gradiometer  studied  in  Subsection  5.2.1;  if  the  gradiometer  is  visualized  as  being  laid  flat 
on  the  ground,  as  it  was  in  the  original  situation  (Chapter  1),  the  dipole  can  now  be  considered 
to  be  raised  in  the  air  a  distance  \z\,  with  everything  else  unchanged.  Because  of  its  orientation, 
the  gradiometer  measures  dByjdx  which,  according  to  our  evaluated  gradient  matrix,  is  given 
by 


dBy  _  ZKx.  2 

~dl  ~ 


(5.27) 


This  expression  can  also  be  put  in  the  form 


dBy  3Kx  7  .  2 

cos'  —  A  tan^  ( 

dx  cP 


(5.28) 


where  <t)  is  the  angle  shown  in  Figure  5.6  (it  is  the  three-dimensional  equivalent  of  the  (j)  shown  in 
Figure  5.1)  and  d  =  (x^  +  z^)^  is  the  perpendicular  distance  from  the  center  of  the  gradiometer 
to  the  y-ax\s  (it  is  the  three-dimensional  equivalent  of  dx  in  Figure  5.1). 

The  <t>  dependence  in  Equation  5.28  is  the  same  as  that  encountered  in  the  equivalent 
two-dimensional  configuration  (see  Equation  5.5);  there  is  therefore  no  difference  between  the 
three-dimensional  and  two-dimensional  angular  responses.  There  is  a  difference  in  the  distance 
dependence,  however,  and  the  magnitude  of  the  response  is  reduced  in  the  three-dimensional 
case  compared  with  its  two-dimensional  equivalent  because  of  the  presence  of  the  distance 
term  d  (d>x).  This  illustrative  example  shows  that  the  response  curves  given  in  Figure  5. .3 
are  relevant  not  just  to  two-dimensional  dipole/gradiometer  configurations  but  also  to  three- 
dimensional  configurations  as  well. 


To  end  this  section  on  the  three-dimensional  approach,  let  us  consider  what  happens  to 
the  evaluated  gradient  matrix  in  Equation  5.26  when  we  take  z  =  0.  The  gradient  matrix  then 
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applies  to  the  two-dimensional  situations  considered  earlier  in  this  chapter,  and  it  has  the  form 


1C!]2=o  = 


(a;2  +  y^)l  0 


y(y^  -  4x^)  x(x^  -  4y^) 
x(x^  -  4y^)  y(3x^  -  2y^) 


(5.29) 


y(x^  +  j/2)J 


The  last  row  of  this  evaluated  gradient  matrix  gives  us  the  following  equations 


_ 

dz  dz 
dB^  ^  3Ky 

(x2 


(5.30) 


which,  with  appropriate  modifications,  are  the  same  as  those  used  in  Subsection  5.2.3  (General 
ized  Two-Dimensional  Approach;  see  Equations  5.16b). 


§5.4  Applications  of  the  Response  Data 

The  various  two-dimensional  response  curves  and  other  wo-dimensional  response  data 
given  in  this  chapter  may  appear  to  be  too  idealistic  to  be  of  much  use  in  practice.  If  so, 
the  appearance  is  misleading.  In  fact,  the  two-dimensional  data  have  immediate  application 
to  a  number  of  practical  situations  in  both  two  and  three  dimensions,  and  even  when  they 
do  not  have  immediate  application,  as  would  be  the  case  in  a  three-dimensional  situation  with 
arbitrary  orientation  of  the  gradiometer,  they  can  often  be  used  to  give  an  approximate  picture 
of  the  gradiometer  response  pending  a  complete  computation.  One  application  of  the  response 
data  has  already  been  discussed— it  covers  the  response  of  the  simple  gradiometer  described 
in  the  first  chapter.  Another  possible  application  is  suggested  by  the  physical  configuration 
of  the  coaxial  gradiometer:  it  would  be  much  easier  to  tow  a  coaxial  gradiometer  through  the 
sea  than  it  would  be  to  tow  one  of  the  paraxial  varieties.  The  response  functions  derived  for 
the  coaxial  gradiometer  could  then  be  used  to  interpret  the  response  of  the  gradiometer  to 
the  magnetic  field  from  submerged  sources  under  the  following  conditions:  (1)  The  magnetic 
moments  of  the  sources  are  induced  by  the  geomagnetic  field  and  are  therefore  predominantly 
North-South  directed,  and  (2)  the  gradiometer  is  towed  either  in  North-South  or  East-West 
directions. 

The  comparatively  large  minimum  obtained  with  the  coaxial  gradiometer  when  the  dijwle 
motion  takes  place  in  the  z-direction  (Figures  5.1  (b)  and  5.3  (d))  suggests  use  of  the  configu¬ 
ration  for  monitoring.  An  illustrative  use  would  be  the  monitoring  of  vessels  through  a  narrow 
East- West  oriented  strait  such  as  the  Strait  of  Gibraltar  by  means  of  a  North-South  orientxxl 
coaxial  gradiometer  located  on  one  side  of  the  strait.  This  configuration  would  provide  the 
largest  possible  signals.  If  signal  amplitude  was  not  as  important  a  factor  as  directivity,  one 
of  the  configurations  with  a  positive-negative  response  (Figures  5.3  (b)  and  (c))  couK!  be  used. 
Squaring  of  the  output,  or  simply  taking  the  absolute  value,  would  give  a  sharp  minimum  at 
the  time  of  closest  approach. 


(luipU'T  () 

Rotating  Component  Gradiometers 


As  I  pxiinted  out  in  Chapter  4,  the  first  two  magnetic  field  gradiometers  desoril)ed  in 
the  scientific  literature  used  rotating  induction  loops  as  their  basic  magnetic  field  sensors. 
Later  gradiometers  did  not  follow  their  lead,  however,  and  no  gradiometers  with  rotating 
components  appear  to  have  been  constructed  for  many  years.  There  are  compelling  engineering 
reasons  for  avoiding  use  of  a  gradiometer  with  moving  components  when  gradiometers  of 
similar  or  better  sensitivity  are  available  that  have  no  moving  parts,  and  the  early  gradiometers 
had  no  distinctive  advantages  that  would  encourage  their  further  development  once  fluxgate 
gradiometers  became  available.  Nowadays,  with  vastly  more  sensitive  gradiometers  available, 
a  gradiometer  with  moving  components  would  have  to  have  some  particularly  distinctive 
properties  to  be  considered  practicable. 

In  this  chapter  I  discuss  two  magnetic  field  gradiometers  with  rotating  components  t  hat 
might  have  sufficiently  distinctive  properties  for  them  to  be  useful.  In  both  cases  they  differ 
from  the  early  rotating  loop  gradiometers  in  that  their  sensor  assemblies  are  rotated  as  a  whole. 
In  one  of  these  gradiometers  it  is  immaterial  what  kind  of  magnetometers  are  cont  ained  in  t  he 
sensor  assembly,  provided  they  measure  components  of  the  magnetic  field,  are  not  induction 
loop  magnetometers,  and  that  they  are  not  incapacitated  by  motion.  Such  a  gradiometer 
may  have  useful  directional  properties.  The  other  gradiometer  requires  induction  loops  as 
sensors.  Its  response  to  magnetic  dipole  fields  falls  off  either  as  the  inverse  fourth-  or  fifth-power 
of  distance,  depending  on  how  the  loops  are  connected,  but  it  also  depends  on  the  rotation 
frequency.  Thus  its  sensitivity  can  be  changed  simply  by  varying  the  rotation  frequency  and  in 
principle  it  can  be  made  very  large.  Neither  of  these  two  rotating  gradiometers  has  ever  been 
constructed,  to  the  best  of  my  knowlege,  and  they  are  quite  speculative. 

§6.1  Rotating  Component  Gradiometer 

Suppose  we  return  to  the  two-dimensional  .single  axis  paraxial  gradiometer  configuration 
that  we  first  considered  in  Section  5.2.1,  and  we  now  rotate  it  until  its  sensitivity  axis  is  oriented 
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Figure  6.1.  Diagram  showing  the  coordinate  system  and  gradiometer  geometry  that  are  u.sed 
in  the  derivation  of  the  response  of  a  rotating  paraxial  gradiometer. 


at  an  arbitrary  angle  (90P  —  (f)  relative  to  the  axis  of  the  magnetic  dipole,  as  shown  in  Figure 
6.1.  With  the  ultimate  object  of  obtaining  information  about  the  variation  of  the  response  of 
the  gradiometer  as  its  sensor  assembly  is  rotated  about  its  center,  let  us  derive  the  response  of 
the  gradiometer  for  this  arbitrary  orientation. 

Resolving  the  two  vector  components  Br  and  8$  along  the  sensitivity  axis  of  the  grad¬ 
iometer,  we  find  that  the  resolved  component  of  the  magnetic  field  at  the  center  of  the 
gradiometer,  Bg,  is  given  by 

Bg  =  —^{2  cos  0  sin(0  +  V?)  +  sin  0  cos(0  -h  ip)),  (6. 1 ) 

where  K  =  fioM/4n  as  in  the  previous  chapter.  The  response  of  the  gradiometer  is  given  by 


ABg=^-^Ar+^-^A0, 


where,  from  Ek]uation  6.1,  we  have 
dBg  _  ZK 


= - 3“  (2  cos  0  .sin(61  +  <p)  +  sin  0  cos{0  -I-  p)), 

or  r* 


BR  *\K 

{cos0  cos{0  +  p)  —  sin  0  sin(0  -I- 1^)). 

The  incremental  quantities  Ar  and  A0  are  given  by 

Ar  =  t  cos(0  -I-  p), 


A0  =  sin(^  +  p). 
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9  (Degrees) 


Figure  6.2.  Plot  of  the  response  function  /(«^,  45°)  for  the  paraxial  gradiometer  shown  in  Figure 
6.1.  The  dipole  source  is  assumed  to  move  in  the  positive  j/-direction  and,  as  indicated  by 
the  notation  for  the  response  function,  the  sensitivity  axis  of  the  gradiometer  is  taken  to 
be  at  an  angle  ^  =  45°  to  the  i-axis. 


where  the  incremental  change  is  measured  from  sensor  1  to  sensor  2  and  i  is  the  gradiometer 
spacing,  which  is  assumed  small  compared  with  the  distance  r  to  the  dipole. 

Combining  Ek]uations  6.2-6.4  and  simplifying,  we  obtain  the  following  expression  for  the 
response  of  the  gradiometer 

‘iK’f 

(6.5) 

where  the  angular  function  fiO,tp)  has  the  form 

f{0,  ip)  =  —(3  cos  0  sin  2{0  +  ^)  +  2  sin  0  cos  %0  +  p)).  (6.6) 

which  may  also  be  written 

f{0,p)  =  cos  0(7  sin*  0  —  3cos*0)  sin2ip  +  2  sin0(sin*0  —  4cos*0)  cos  2^.  (6.7) 

This  last  expression  for  the  angular  function  is  somewhat  more  convenient  for  computation. 
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These  latter  equations  for  the  response  of  a  paraxial  gradiometer  reduce  to  prt'viously 
derived  equations  when  ip  =  (f,  and  9flP.  They  also  provide  information  about  the  response 
for  any  other  two-dimensional  configuration:  they  represent  another  generalization  of  the 
two-dimensional  approach  I  introduced  in  Section  5.2,  and  they  are  also  relevant  in  a  thrc'e- 
dimensional  dipole/gradiometer  context  because  of  the  applicability  of  the  two-dimensional 
data  to  many  three-dimensional  problems.  To  illustrate  the  use  of  the  new  response  informa¬ 
tion,  I  will  now  derive  the  response  of  the  paraxial  gradiometer  when  p  =  45°  for  dipole  motion 
along  the  y-axis  (Figure  6.1). 

If  V?  =  45°  in  Equation  6.7,  we  have 

/(0,45°)  =  cos(9(7sin2  0-3cos2  0),  (6.8) 

which,  changing  from  0  to  0  to  conform  with  the  work  in  Section  5.1,  can  also  be  written 

/(,^,45°)  =  sin^<^(7tan2(^-3).  (6.9) 

Choosing  motion  of  the  dipole  in  the  positive  y-direction,  without  any  change  in  its  orientation, 
and  taking  r  =  d^cos"*  <j>,  we  obtain 


=  (6.10) 

where  the  response  function  /(<!>,  45°)  is  given  by 

f{4,,  45°)  =  cos'*  0  sin^  (7  tan*  -  3).  (6. 1 1 ) 

A  plot  of  this  response  function  is  shown  in  Figure  6.2. 

The  generalization  I  have  introduced  above  can  be  easily  extended  to  coaxial  gradiometers 
and  thus  the  two-dimensional  response,  or  generalized  two-dimensional  response  (as  discu.ssed 
in  Section  5.2.3),  of  any  arbitrarily  oriented  component  gradiometer  to  a  dipole  magnetic  field 
can  be  calculated.  This  is  the  last  generalization  I  will  make,  since  I  have  now  provided  all 
the  information  necessary  to  compute  the  response  of  gradiometers  to  dipole  fields  in  any  two- 
or  three-dimensional  configuration.  Let  us  now  consider  what  happens  to  the  response  of  the 
paraxial  gradiometer  considered  above  when  it  begins  to  rotate  about  its  center. 

To  investigate  the  rotational  response  of  the  gradiometer,  suppose  we  return  to  Equations 
6.5  and  6.7  and  compute  f{0,p)  for  0  =  (f,  45°,  and  90P,  and  p  =  ujt,  i.e.,  we  allow  the 
gradiometer  sensor  to  rotate  at  angular  frequency  w  (period  T)  for  three  different  values  of 
its  position  angle  0.  If  we  assume  r  is  kept  fixed,  the  physical  situation  now  being  analyzed 
is  equivalent  to  having  a  dipole  source  moving  in  a  circle  around  the  rotating  gradiometer 
sensor  while  maintaining  the  orientation  of  its  moment  in  a  fixed  direction.  This  is  not  a 
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Figure  6.3.  Plots  of  the  function  for  the  rotating  component  gradiometer,  where  <p  is 

allowed  to  vary  with  time  t  according  to  the  relation  ^  =  ut  and  0  is  given  the  three  fixed 
values  5  =  OP  (top),  45°  (middle),  and  9ClP  (bottom).  Also  shown  (dashed  lines)  are  parts 
of  the  variation  of  f{9,  f)  for  0  =  33.2P,  and  63.4°. 


particularly  idealistic  situation  and  the  response  of  the  gradiometer  will  provide  information 
about  its  possible  practical  applications. 

Figure  6.3  shows  the  variation  of  the  response  function  f(0,  ip)  under  the  a&sumed  condi¬ 
tions.  It  can  be  seen  that  there  is  some  variation  in  the  amplitude  of  the  response  functions 
for  the  various  values  of  9,  but  the  change  of  particular  relevance  to  us  is  the  shift  of  the 
maximums  in  the  response  curves  to  .smaller  values  of  t  as  0  increases  from  (f  to  QOP.  For 
example,  the  peak  that  originally  occurred  at  <  =  3T/8  (p  =  135°)  for  9  —  (f  has  moved  to 
t  =  OP  (p  =  OP)  for  9  =  9CP.  Although  there  is  not  a  1:1  relationship  between  the  changes  in 
9  and  the  peak  in  the  response  curve,  the  data  suggest  that  it  should  be  possible  to  track  I  hc 
change  in  bearing  of  a  moving  dipole  source  using  a  rotating  component  gradiometer.  Needless 
to  say,  the  data  from  the  rotating  gradiometer  ought  to  be  proces.sed  by  means  of  a  minicom¬ 
puter  to  take  full  advantage  of  the  theoretical  knowledge  of  the  gradiometer’s  response,  but 
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this  processing  should  not  involve  any  new  techniques  and  need  not  be  discussed  further  here. 
Further  research  on  this  passive  technique  for  tracking  magnetic  dipole  sources  could  lead  to 
new  MAD  instrumentation. 


§6.2  Rotating  Induction  Loop  Gradiometer 

To  introduce  the  concept  of  a  rotating  induction  loop  gradiometer  (RILG),  I  will  once 
again  consider  a  comparatively  simple  example  and  then  generalize  from  it  to  more  involved 
situations.  Let  us  start  therefore  by  considering  the  RILG  configuration  shown  in  Figure  6.4. 
The  RILG  itself  consists  of  two  physically-connected  coplanar  loops  of  N  turns  and  area  A, 
which  are  separated  by  a  distance  £  (measured  between  centers)  in  the  plane  of  the  loops  and 
which  rotate  about  an  axis  passing  through  the  midpoint  of  the  line  joining  their  centers  (the 
baseline).  Different  varieties  of  RILG’s  can  be  constructed  depending  on  the  orientation  of  t  he 
rotation  axis  relative  to  the  baseline  and  the  plane  of  the  loops;  I  will  assume  here  that  the  axis 
is  perpendicular  to  the  plane  of  the  loops.  Thus,  in  the  figure,  the  two  loops  rotate  in  a  circle 
about  the  center  of  the  line  joining  them.  The  magnetic  dipole  source  is  located  a  distance 
r  away  from  the  center  of  the  loop  assembly  and  it  will  be  assumed  to  be  directed  vertically 
upwards,  i.e.,  its  moment  is  directed  upwards  out  of  the  plane  of  the  figure  and  is  parallel  to 
the  axis  of  rotation  of  the  RILG. 

Because  of  the  assumed  geometry,  the  magnetic  field  of  the  dipole  is  everywhere  directed 
perpendicularly  downwards  into  the  loop  plane.  It  has  no  angular  variation  and  its  magnitude 
is  given  by 


which  follows  from  Ek]uations  5.1  when  0  =  The  quantity  K  is  given  hy  K  =  fA-n 
as  in  the  previous  section.  The  loop  spacing  I  is  assumed  to  be  very  small  compared  with  r, 
and  the  loop  areas  A  are  also  taken  to  be  so  small  that  the  magnetic  field  can  be  considered 
to  be  uniform  over  each  area.  With  these  assumptions  the  magnetic  flux  passing  through  the 
coils  can  be  written 

^  =  B.dS  =  -NAB.  (6. 1 .3) 

The  emf’s  induced  in  the  loops  are  given  by  —d^/dt  and  —d^/dt,  which  can  also  be  written 
NAdBifdt  and  NAdBi/dt,  where  the  subscripts  indicate  which  loop  is  involved.  We  now 
need  to  derive  B\  and  B2  in  order  to  compute  the  emf’s. 

If  we  write  ri  and  r2  for  the  distances  from  the  dipole  to  the  centers  of  loop  1  and  loop  2, 
respectively,  we  have,  to  a  first  approximation. 


Figure  6.4.  Diagram  showing  the  geometry  used  in  the  analysis  of  the  response  of  a  rc'tating 
induction  loop  gradiometer  (RILG).  In  the  particular  configuration  shown  the  dipole  source 
is  directed  vertically  upwards  out  of  the  plane  of  the  figure,  which  also  happens  to  be  (1) 
the  plane  containing  the  two  loop  sensors  and  (2)  the  plane  of  rotation  of  the  loop  assembly. 


giving 


Si  = /Cr  +  ^cos(<^- , 
B2  =  Kr~^  ^  • 


(6.15) 


Using  a  binomial  series  expansion  and  omitting  third  and  higher  order  terms,  Flquations  6.15 
can  be  written 


Bi  =  Kr  ^  ^  ^  cos^(<^  -  ip)^, 

I A  3^  ^  3^2  \ 

I  1  +  —  C08(^  ~  4^  vos‘(<i>  -  (p)  1. 


(6.16) 


B2  =  Kr 


Since  the  loops  are  rotating,  we  can  write  ip  =  u/t,  where  w  is  the  angular  frequency  and  t  the 
time,  in  Equations  6.16.  When  differentiated  with  respect  to  time,  and  then  .subst  itu(<Hl  in  (  he 
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expressions  for  the  induced  emf’s,  B\  and  B2  give 

,  mAKCu 
f'mfi  =  — - 


sin(<^  —  uit)  +  ^  sin  2((^  — 

^  iNAKCi^f  .  ,  ,  ^  .  A 

emf2  =  - - I +sin(0  —  a><)  +  —  sin  2((t>  —  u!t)j. 


(6.16) 


If  the  two  loops  are  now  connected  in  series  opposition,  the  output  of  the  rotating  loop 
sensor  assembly  will  be  a  voltage  Vgi  given  by 


,,  3NAKe^  .  ,, 

Vg\  =  - 4 -  Sin(^  -  Ujt), 


(6,17) 


or,  if  the  are  connected  in  series  with  the  emf’s  adding,  the  output  of  the  sensor  assembly  will 
be  a  voltage  Vg‘i  given  by 


3NAKe‘^oj 


sin  2(<p  —  (jjt). 


(6.18) 


The  RILG  can  therefore  respond  to  the  magnetic  field  of  the  dipole  in  one  or  the  other 
of  two  modes.  The  first  mode  can  be  considered  to  be  a  first-order  gradiometer  mode,  since 
the  coils  are  connected  conventionally  in  series  opposition  and  the  output  of  the  RILG  is 
proportional  to  the  inverse  fourth  power  of  the  distance  to  the  source  as  in  other  first-order 
gradiometers.  However,  the  output  of  the  gradiometer  is  far  from  orthodox  in  other  ways,  and 
most  particularly  in  the  way  it  consists  of  a  sinusoidal  variation  with  an  amplitude  that  can  be 
made  larger  by  increasing  the  gradiometer  parameters  N,  A,  and  C,  and  above  all  by  increasing 
w.  Thus  the  sensitivitj  of  a  RILG  can  be  made  large  by  spacing  its  loops  far  apart,  by  giving 
them  large  areas  and  r  lany  turns,  and  by  giving  the  sensor  assembly  a  high  rotation  rate.  It  is 
of  interest  that  the  sensitivity  can  be  varied  continuously  within  wide  limits  by  varying  w:  t  his 
control  could  have  many  practical  advantages.  Note  that  the  dependence  of  the  sensitivity  on 
w  and  (  are  perhaps  more  important  than  the  dependence  on  N  and  A,  because  increasing  the 
latter  quantities  will  lead  to  an  increase  of  thermal  noise  in  the  loops,  whereas  an  increase  in 
the  former  quantities  will  not  lead  to  such  an  increase. 

The  second  mode  can  similarly  be  considered  to  be  a  second-order  gradiometer  mode 
because  of  the  dependence  of  its  output  on  the  inverse  fifth  power  of  distance  to  the  source. 
Its  sensitivity  also  depends  on  u  and  can  therefore  be  varied  comparatively  easily.  It  could  be 
used  where  better  discrimination  against  nearby  sources  of  magnetic  noise  is  required.  Since 
there  .should  be  little  difficulty  switching  the  connections  between  the  two  loop  sensors,  a  RILG 
is  potentially  two  gradiometers  in  one:  a  first-order  version  for  maximum  range  and  a  second- 
order  version  for  better  characterization  of  the  sources  at  closer  distances  in  the  presence  of 
noise. 


Another  dipole/RILG  configuration  that  can  be  analyzed  with  little  difficulty  using  basic 
principles  is  essentially  the  same  as  the  one  we  have  just  studied  but  with  the  important 
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difference  that  the  dipole  moment  lies  in  the  plane  of  rotation  of  the  loops.  For  this  configuration 
the  magnetic  flux  passing  through  the  loops  is  zero,  and  thus  the  respon.se  of  the  RILC  is  also 
zero.  This  may  appear  to  be  a  trivial  example,  but  it  indicates  a  potential  weakness  of  the 
RILC:  it  will  not  always  respond  to  a  dipole  field.  To  ensure  a  response  under  all  circum.stances 
it  may  be  necessary  to  use  two  crossed  RlLGs.  However,  it  should  be  obvious  that  the  exact 
.  conditions  for  zero  response  ''f  the  RILC  are  likely  to  occur  only  rarely  in  practice,  and  if  they 
do  oc'cur  it  is  probable  that  they  will  only  last  for  a  short  time. 

As  a  partial  generalization  of  the  above  results,  I  must  now  admit  to  excessive  simplification 
in  the  use  of  the  distance  approximations  in  Equations  6.14.  More  accurate  expressions  for  t  he 
distances  are 


n  =  r+  -  cos{^  -  tP)  -  3-  sin'^(</>  -  <p), 

z  or 

£  £ 

r2  =  r-  ~  cos(<^  8^  sin^(<^  -  ip), 


(6.19) 


and  when  these  distance  expressions  are  used  in  the  derivation  of  the  RILC  resjmn.se  for  the 
vertically  oriented  dipole,  it  is  found  that 


= - ;:4 - -‘PI  -  ^  2(<ji  -  p)  1, 

NAKe‘^<jj 


^92  =  —5 


sin  2{<j)  —  p). 


(6.20) 


These  more  accurate  expressions  for  the  voltages  produced  by  the  rotating  loop  as.sembly  in  the 
RILC  show  that  the  response  of  the  RILC  connected  in  series  opposition  is  actually  a  mixture 
of  the  first-order  response  derived  previously  with  an  additional  second-order  response.  This 
appears  to  be  a  general  feature  of  the  series-opposition  RILC  response.  However,  it  will  be 
immediately  apparent  that  the  second-order  response  will  be  small  and  probably  negligible 
compared  with  the  first-order  response  under  normal  circumstances  of  use  (i.e.,  r  3>  f)  because 
of  its  additional  i/r  factor.  The  response  of  the  RILC  connected  in  the  series  aiding  mode 
remains  essentially  the  same  second-order  response  as  before,  with  only  a  small  numerical 
change  in  the  magnitude  of  the  voltage  produced  by  the  loop  as.sembly.  This  purely  second- 
order  response  of  the  series-aiding  RILC  also  appears  to  be  a  general  feature. 

My  colleague  D.  M.  Bubenik  has  derived  general  expressions  for  the  first-  and  second-order 
RILC  responses.  The  first-order  term  has  the  form 


CiN X  g).grad  (B.s)), 


where  C\  is  a  constant,  and  u>,  g,  and  s  are  all  unit  vectors  indicating  the  sense  of  rotation, 
the  direction  of  the  line  joining  the  two  sensors,  and  the  direction  of  the  .sensitivity  axes  of  (he 
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loops.  The  second-order  term  has  the  form 

C2NAijj£^  g .  grad  ((w  X  g) .  grad  (B.s)), 

where  C2  is  another  constant  and  the  same  unit  vectors  are  used  as  described  above.  The.s(' 
general  terms  can  be  used  to  investigate  the  response  of  the  RILG  under  different  conditions 
from  those  I  have  discussed  above. 

The  response  of  the  RILG  differs  intrinsically  from  the  responses  of  the  other  gradiometers 
considered  in  this  report.  Unlike  the  other  responses,  it  consists  of  a  periodic  variation  with 
time  (usually,  but  not  always,  a  purely  sinusoidal  variation)  and  it  is  the  amplitude  of  this 
variation  that  relates  to  the  properties  of  the  dipole.  From  the  operational  point  of  view  this 
could  have  some  advantage:  when  no  dipole  sources  are  in  range  the  output  of  the  RILG  should 
be  flat  (assuming  that  the  largely  known  gradient  of  the  geomagnetic  field  has  been  removed  by 
real  time  data  processing-  essentially  a  trivial  procedure  with  modern  microcomputers),  and 
it  is  only  when  a  source  comes  in  range  that  any  periodic  variation  will  be  seen.  The  distance 
from  the  source  to  the  RILG  and  the  dipole  moment  are  both  combined  in  the  amplitude  of  the 
variation  and  they  probably  cannot  be  separated  without  additional  information.  Directional 
information  is  available,  even  though  this  may  not  be  immediately  obvious  from  the  expressions 
I  have  given  above.  The  directionality  comes  from  the  angle  (^>  in  Equations  6.17  and  6.18,  for 
example,  which  implies  that  the  direction  to  the  source  can  be  obtained  from  the  times  of  the 
zeros  or  of  the  maximums/minimums  in  the  periodic  variation.  There  may  be  some  ambiguity 
in  the  direction  deduced  from  this  information,  but  this  problem,  like  the  others  I  have  pointed 
out,  can  probably  be  eliminated  by  further  study  and  modifications  to  the  measurement  system. 


§6.3  Conclusion 

As  I  noted  at  the  beginning  of  this  chapter,  the  Rotating  Component  Gradiometer  and 
the  RILG  are  more  speculative  gradiometers  than  the  others  1  have  considered  in  this  report. 
However,  I  should  also  note  that  they  are  speculative  at  this  stage  entirely  becau.se  they  are 
new  and  have  not  been  studied  to  the  same  extent  as  other  gradiometers.  They  are  presented 
here  partly  vO  show  that  gradiometers  departing  from  the  conventional  form  are  possible,  and 
partly  because  they  contrast  with  the  conventional  gradiometers,  thus  providing  a  fresh  view 
of  the  properties  of  the  conventional  instruments.  Above  all,  the  rotating  gradiometers  are 
described  to  show  that  new  ideas  and  concepts  are  possible  in  the  field  of  gradiometer  studies 
and  to  encourage  further  research  on  gradiometers. 


(’hiiptcT  7 

Discussion 


This  report  cites  many  applications  of  magnetic  field  gradiometers  in  areas  of  particular 
importance  to  the  United  States:  energy  resources,  defense,  and  health  care.  Despite  lhes<' 
applical  ions,  t  here  is  a  notable  lack  of  articles  in  the  recent  scientific  literature  describing  basic 
restvirch  relevant  to  magnetic  field  gradiometers.  There  are  a  few  articles  describing  t  he  n'sidts 
of  basic  rest'arch  pertaitiing  to  superconducting  gradiometers,  but  even  in  this  case,  where  t  he 
opportunities  are  so  great,  the  scientific  lit.erature  suggests  that  basic  research  is  not  well  sup¬ 
ported.  It  is  particularly  disturbing  to  find  so  few  mejisurements  of  the  spat  ial  uniformity  of 
natural  magnetic  noise  with  magnetometer  arrays  and  gradiometers,  since*  such  measurements 
are  es.sent  iai  to  the  informed  application  of  magnetic  gradiometer  and  diffen'iil  iai  magnetomet  er 
techniques.  It  is  surprising  to  find  so  few  attempts  to  measure  natural  magnetic  graelients,  par- 
ti<*ularly  gradients  associated  with  ionospheric  currents,  either  with  high-sen.sit,ivity  supe'rcon- 
ducting  gradiometers  or  with  long-baseline  total  n<*ld  gradiometers.  It  was  pointed  out  many 
years  ago  by  Chapman  and  Nelson  [1957)  that,  measurements  of  these  natural  gradients  could 
provide  imjxjrtant  information  about  the  overhead  electric  current  systems;  tin*  measurernc'tils 
would  also  hel[)  characterize  a  source  of  gradient  noise.  Despite  the  revolution  currently  tak¬ 
ing  place  in  the  use  of  small  computers  for  real-time  computation,  t  here  is  litt  le  evidence*  of 
this  revolution  having  a  significant  impact  on  present  gradiometer  techniques.  In  the  fut  ure 
I  would  expect  to  see  small  computers  incorporated  into  the  racks  of  electronic  e(|uipm(*nl 
at  each  gradiometer  installation,  or  even  included  in  the  gradiometer  sen.sor  as.sembli(*s  them- 
.s<*lves,  with  the  gradiometers  being  given  greatly  enhanced  capability  as  a  r(*sult .  In  short . 
there  are  many  opportunities  for  innovative  basic  and  applied  research  either  specifically  oji 
gradiometers  or  in  relat  ion  to  them.  and.  in.steaci  of  mer(*ly  noting  that  fact  and  leaving  it  to 
the  reader  to  work  out  what  the  opportunities  really  are.  I  will  devote  the  remainder  of  this 
chapter  to  a  listing  of  the  ar<*jis  where  I  perceive  a  n<*ed  for  further  res(*arch.  or  wh(*re  I  l>eli(*v(* 
r(*search  could  provide  new  results  aiul  capabilities. 
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7.  I)is<'ussi()n 


§7.1  Suggested  Research  Areas 


7.1.1  Total  rich!  MagiK'lomrtcr  Improvement 

The  first  art'a  wliore  research  appears  to  f)e  needed  is  a  broad  one  involving  all  asp('(  ts  of 
total  field  magnet omt'ters:  their  basic  physics,  iinprovcanents  to  their  performance,  and  their 
use  for  nu'iisurernents  in  the  natural  environment.  I  have  deliberately  given  first  listing  to 
this  area  of  needed  research  because  of  its  particular  relevance  to  energy  resource  development 
and  defense,  as  I  will  now  describe  in  greater  detail.  Alt  hough  the  total  field  inst  ruments  do 
not  have  the  high  nominal  sen.sitivity  of  superconducting  miignetometers  and  gradiometers  at 
rest,  they  an*  still  the  most  sensitive  instruments  of  their  kind  to  have  been  used  routinc'ly 
on  aircraft,  and  they  are  likely  to  remain  the  most  sensitive  for  sonn*  time  to  com('.  since  tlu' 
pract  icability  of  t  he  airborne  operation  of  superconducting  magnetometers  and  gradiometers 
has  yet  to  be  den\onstrated.  I  single  out  this  airborne  capability  because  of  its  pertinence  to 
geophysical  pros{)ecting  and  to  submarine  localization.  When  considered  wit  hin  t  heir  wider 
contexts  of  energy  resource  development  and,  in  particular,  of  antisubmariiu'  warfare  [Wit, 
1981],  these  two  activities  have  great  strategic  significance  at  the  present  time.  Hecause  of  t  his 
significance,  it  would  be  expected  that  there  would  be  a  continued  high  level  of  re-search  on 
total  field  magnetometers  and  gradiometers.  Unfortunately,  the  level  of  research  in  t  he-  U.  S. 
not  only  fails  to  meet  this  expectation,  but  it  is  at  so  low  a  level  that  a  major  failure-  e)f  the- 
researeh  management,  process  in  the  U.S.  is  indicated.  Consider  the  following  elata: 

The  journal  Physics  Abstracts  lists  the  abstracts  of  pract  ically  all  se  ientific  papers  publisheel 
in  the  worlel  on  what  is  very  broadly  classified  as  physics  topics.  In  view  of  the  pre-ssure  on 
scientists  throughout  the  world  to  publish,  the  number  of  article's  (and  thus  the  number  of 
abstracts)  on  a  partie-ular  physics  topic  is  likely  to  be  a  good  measure-  of  the  level  of  re-seare-h 
activity  involving  the  topic,  bet  us  now  look  at  the  numbe-r  e)f  abstrae-ts  of  iirticle-s  publishe-el 
during  the  feair-ye-ar  interval  1977  1980  that  describe  re-search  e-ither  em  or  involving  total 
fie-lel  magnetometers.  My  .source  for  this  information  is  the  i.ssue  of  Physics  Abstracts  that 
specifically  covers  the  indicate'd  four-year  interval,  and  I  find  that  of  the  232  magne-tome-ter 
references  listed  there  are  28,  or  12%  of  the  whole,  that  involve  work  e-it,her  on  or  wit  h  t  otal  field 
instruments.  Of  these-  28.  2r)  were  references  to  articles  by  Soviet  scientists  (89% of  the  total 
fie-lel  magne-tometer  articles),  2  by  .Japanese  .scienti.sts.  and  1  by  .scientists  in  the  U.  S.  (1%  of  the- 
article-s  on  total  field  magnetometers;  0.4%  of, the  total  number  of  article-s  on  magnet e)meters) 
The-  message  of  the.se  elata  is  unambiguous  and  it  clearly  supports  my  argument  that  more- 
research  is  ne-eded  on  total  field  magnetometers  in  the  U.  S. 

It  is  intere-sting  to  list  some  of  the  titles  of  the  Soviet  article-s:  (1)  “Knhance-d  aee  iiraey 
of  the  helium  Mz  magnetometer”  [Aleksandrov  and  Yakobson.  1979],  (2)  “Alkali-he-liiim  mag- 
ne-toneete-r”  [Blinov  et  ai,  1979],  (3)  “Sele-ction  of  passbanel  of  ji  magne-tome-te-r  whe-n  tTie-asiiring 
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tho  iiuluclioii  of  a  moving  magnetic  body”  \Chechurina  et  ai,  1976],  ('I)  "High-accuracy  T-MP 
fu'ld  proton  magnetometer”  [Maksimovskikh  and  Shapiro,  1976],  (5)  “Evaluation  of  opt  imum 
conditions  for  towing  a  (proton)  magnetometer”  [Shreyder,  1977],  (6)  “Helium  magnetometer 
with  pulsed  optical  pumping”  [Yakohson  and  Aleksandrov,  1978].  These  exampl(?s  illustrate  the 
breadth  of  the  Soviet  research  effort;  the  remaining  19  artich's  mostly  describe  either  detailed 
studies  of  the  biusic  processes  involved  in  the  op<'ration  of  the  total  field  magnetometers  or 
specific  improvements  t.hat  have  been  achievcHl  in  their  performance. 

The  Soviet  work  specifically  on  magnetic  field  gradiometers  is  much  less  comprehensive 
than  the  r(>search  on  total  field  magnetometers  that  I  have  cited  above  (the  few  available  Soviet 
articles  on  gradiometers  mostly  contain  rather  academic  theoretical  studi(!s  [e.g.,  Reznik,  1967: 
Roze,  1973,  1978]),  but  sensor  improvement  is  probably  the  foremost  requirement  for  improved 
gradiometers  and  the  Soviet  program  of  U)tal  field  magnetometer  improvement  could  as  well 
s('rv('  as  a  model  for  improved  total  field  gradiometers. 

7.1,2  Innovative  Total  Field  Gradiometer  C'onfigiiration.s 

In  addition  to  further  research  on  the  basic  physics  of  total  field  magnetometers,  as 
recommended  above,  it  is  also  desirable  that  a  fresh  look  be  given  to  the  br't.ter  utiliitatinn 
of  total  field  magnetometers  in  gradiometer  configurations.  It  is  well  recognizt'd  t  hat  total 
field  magnetomeU'rs  are  not  as  .sensitive  as  the  most  recent  supercondtjciing  models,  and  as 
a  result  it  is  commonly  assumed  that  total  field  gradiometers  must  be  much  h'ss  sensitive 
than  superconducting  gradiometers.  This  assumption  is  not  nece.s.sarily  correct,  however,  since 
gradiometer  sensitivity  does  not  depend  solely  on  the  sensitivity  of  the  individual  sensors:  the 
spacing  betwecti  the  sensors  is  also  an  important  factor.  Ek'cause  of  their  need  to  b<'  kept  at 
li(|uid  helium  temperatures,  .superconducting  gradiometer  sensors  are  necessarily  limit  ed  t  o  sizes 
that  can  be  iaserted  in  a  dewar  this  compactness  is  generally  cited  as  an  advantage,  as  indeed 
it  is,  but  it  limits  the  sensitivity  that  can  be  achieved.  Total  field  gradiometers  are  not  similarly 
limited  and  any  baseline  can  be  used  that  is  consistent  with  gradiometer  operation.  When  t  he 
combination  of  magnetometer  sensitivity  and  baseline  is  considered,  total  field  gradiome«ers  are 
not  inferior  to  t  he  superconducting  varieties  in  gradiometer  sensitivity.  To  illu.strate,  suppose 
we  have  a  sujierconducting  magnetometer  that  is  1000  times  more  sensitive  than  the  l>est 
available  total  field  magnetometer,  and  that  two  of  these  superconducting  magnetometers  are 
combined  to  produce  a  gradiometer  with  a  baseline  of  30  cm,  which  is  representative  of  the 
Ionge*st  baselines  that  have  been  used  up  to  the  present  time.  To  match  th<'  sensitivity  of  t  his 
superconducting  gradiometer,  a  gradiometer  constructed  from  the  total  field  magnetometers 
would  need  to  have  a  baseline  of  1000  X  30  cm,  or  300  m.  Such  baselines  are  quite  feasible,  and 
have  already  be(’n  used.  In  fact,  there  appears  to  be  no  reason  why  much  longer  baselines  cannot 
be  utilized,  since  natural  geomagnetic  field  fiuctuations  are  likely  to  be  relatively  uniform  over 
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distances  mucti  greater  than  300  m  (but  measurements  of  this  uniformity  are  reciuired.  as 
detailed  below). 

For  practical  reasons,  the  total  field  gradiometers  that  have  b('en  developed  for  airboriK' 
use  have  typically  had  a  vertical  configuration,  with  the  two  magiu'tometers  being  towed  om* 
above  the  other  on  a  single  cable.  Another  possible  configuration  that  could  be  develop<'d  is  a 
synthesized  horizontal  gradiometer,  which  wouhl  nMiuire  two  separate  total  field  magnetometers 
to  be  towed  at  approximately  the  same  altitude  by  two  separate  aircraft  .  Large  spacings  would 
be  possible  with  this  configuration,  giving  it  a  high  gradient  sensitivity  and  making  it  mor(' 
effective  for  surveying  large  areas  than  the  vertical  gradiometer.  The  important  [)art  of  this 
suggestion  is  the  requirement  t  hat  the  spacing,  altitude,  and  responses  of  t  he  magnetonu'ters  be 
monitored  in  real  time  and  the  gradiometer  res|K)ns<'  of  the  combination  be  synthesized,  using 
a  small  computer,  either  on  one  of  the  aircraft  or  on  the  ground.  The  spacing  and  altitude  can 
be  obtained  by  incorporating  transceivers  and  altimeters  in  the  magnetonu'ter  housings  and 
th('  data  fed.  iti  real  t  ime,  to  the  computer.  Since  one  of  the  important  results  to  be  obt aim'd 
during  a  surve'y  is  the  ident  ification  and  localization  of  sources  of  magnetic  field,  tfu're  is  likely 
to  be  a  transition  from  operation  in  the  gradiometer  mmle  (distance  to  the  source'  large  re'lative' 
to  the  magnetometer  spacing)  to  the  differential  magnetome'ter  mode'  (distance'  to  the  source' 
comparable  to  magnetometer  spacing)  during  many  surve'ys.  The  computer  proce'ssing  e)f  the' 
magne'tometer  elata  e-an  take  this  transition  into  account,  but  it  is  a  further  topie-  that  ne'e'ds 
to  be  stuelie'el. 

7.1.3  (’onipute'r  I'lilization 

By  using  the  t  echniejues  1  have  described  in  this  report,  it  is  pejssible  to  compeit  e  t  he  response 
of  any  cemipeme'tit  or  te)tal  fielel  gradiometer  to  magnetic  dipole  se)ure'es.  The'  e-ornputat  ie)ns 
are  ne)t  exte'tisive  and  there  is  no  rea.son  why  a  small  computer  could  ne)t  be  pre)gramme'el 
to  e-arry  out  the  e-ejinputations  during  a  gradiome'ter  survey  and  thus  enable  the  e)pe'rate)r  to 
compare  t  he  ae-t  ual  measurements  with  theoretical  predictions  in  real  time.  This  is  just  one 
of  many  possible  computer  applications  in  gradiometer  measurements.  In  the  preceding  sub¬ 
section  I  introduced  the  idea  of  using  a  computer  to  synthesize'  a  gradiometer  response'  fre)m 
the'  me'asure'rnents  luaele'  by  twe)  physie-ally-separate'el  total  field  magnetometers.  ('e)mi)uters 
e-oulel  alse)  be  iise'el  te)  remewe  a  re'sponse  to  the  ge^omagnetic  field  graelient  or.  more  ge'nerally. 
to  e'lirninate'  se)ure'e's  e)f  gradient  noise  by  ada|)tive  teH-hniejnes.  In  the  use  of  multi-axis  sujM'r- 
e'e)nelucting  graelie)niete'rs,  computers  are  essential  fejr  the  proe'cssing  of  the  elifferent  streams  of 
elata  so  that  the  grad iome'ters  can  perfe)rm  in  the  “tracker-classifier”  me)de  (Sectiem  l.r)).  which 
may  well  be  erne  e)f  the'ir  most  important  applicatiems.  Clearly,  the'  siicee'ssful  applicatie)n  e)f 
inagne'tic  lie'lel  graelieune'ters  is  likely  to  hinge  on  the  e'llicient  utilization  of  e'omi)ute'rs.  Tenlay. 
many  small  anel  e'asily  pe)rtable  e-ornputers  are  available  that  are  ieleal  for  graeliomete'r  use' 
and  it  is  me)st  eh'sirable  that  tlu'V  be  ine-ejrfmrate'el  inte)  graelie)me'te'rs.  Meale'rn  ceanputational 
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techniques,  particularly  adaptive  techniques,  also  need  to  l)e  studied  with  a  view  to  improving 
gradiometer  performance.  Since  there  is  little  evidence  in  the  recent  scientific  literature  for 
the  use  of  computers  to  analyze  gradiometer  data  in  real  time,  or  to  carry  out  the  higher  level 
procedures  I  have  just  described,  I  recommend  computer  utilization,  and  the  accom|)anying 
use  of  modern  computational  techniques  (including  adaptive  techniques),  as  another  area  for 
further  gradiometer  research. 


7.1.4  Measurements  with  .Superconducting  Ciradiomelers 

Despite  their  great  promise  for  magnetic  field  gradient  studies,  few  superconducting  grad- 
iometers  have  actually  been  used  as  scientific  instruments  for  the  measurement  of  natoral 
gradients  either  gradients  of  the  geomagnetic  field  itself  or  the  gradients  produced  in  the 
geomagnetic  field  by  other  natural  sources.  One  exception,  which  also  indicates  the  kinil  of 
rt'search  that  I  would  recommend,  are  the  measurements  made  by  Gillespie  et  al.  [1977]  and 
Podney  and  Sager  [1979],  first  on  the  noise  in  a  superconducting  gradiometer  at  a  remote 
location  and  then  on  the  magnetic  gradients  produced  by  internal  waves  in  t  he  sea.  Apart 
from  the  valuable  contributions  tht^se  measurements  could  make  in  a  variety  of  geophysic'al 
disciplines,  the  measurements  would  help  establish  a  foundation  of  expertise  in  t  he  scient  ific 
community  that,  in  the  long  run,  is  likely  to  ix'  the  most  important,  factor  in  the  successful 
application  of  superconducting  gradiometers. 


7,1.5  Superconducting  (iradiometer  Development 

As  f  have  mentioned  many  tinais  in  this  refwrt,  one  of  the  key  reejuirements  in  a  com¬ 
ponent  gradiometer  is  accurate  alignment,  of  t,he  component  sensors.  No  matter  what  kind  of 
component  gradiometer  is  involved,  there  always  appears  to  be  some  limit  to  the  degree  of 
alignment  that  can  be  achieved  and  continual  research,  and  new  ideas,  are  needed  for  furt  her 
improvement.  Superconducting  gradiomet<‘rs  are  most  severely  affected  by  this  need  for  align¬ 
ment,  or  balancing,  because  of  the  great  sensitivity  of  the  individual  sensors;  it  is  primarily  t  he 
degree  of  imbalance  that  limits  the  gradiometer  sensitivity.  There  is  already  a  substantial  level 
of  research  on  the  balancing  of  superconducting  gradiometers  [e.^.,  Jaworski  and  Crum,  1980], 
which  1  will  as.sume  will  continue  in  the  future.  However,  there  is  an  important  qiu'st  ion  con¬ 
cerning  balancing  that  arises  when  the  possible  applications  of  superconducting  gradiometers 
are  considered:  how  well  can  such  a  gradiometer  maintain  its  balance  when  airborne,  or  when 
used  on  a  .series  of  flights?  There  are  other  questions  that,  arise  when  the  airl)orne  use  of 
superconducting  gradiometers  is  contemplated,  all  involving  the  motion  of  components  of  t  he 
gradiometer  {e.g.  the  liquid  helium  in  the  dewar)  or  of  the  aircraft,  and  the  way  thes('  mo¬ 
tions  may  degrade  t  he  gradiometer  measurements.  Research  on  the  effects  of  mot  ion  on  a 
su|M'rcon(luct ing  gradiometer  is  therefore  desirable,  particularly  tlu*  effects  of  aircraft  motion. 
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7.  Discussion 


7.1.6  Uniformity  of  Geomagnetic  Noise 

In  several  places  in  this  report  I  have  suggested  that  natural  geomagnetic  field  fluctuations 
are  likely  to  be  uniform  over  distances  ranging  from  hundreds  of  meters  to  hundreds  of 
kilometers.  This  suggestion  is  in  agreement  with  the  few  measurements  that  have  been  made 
\e.g.,  Zelwer  and  Morrison,  1972]  and  it  is  one  of  the  basic  reasons  for  the  usefulness  of  magnetic 
gradiometers  and  differential  magnetometers.  However,  there  is  an  urgent  need  for  better  data 
on  the  spatial  uniformity  of  the  fluctuations.  In  particular,  the  variations  of  the  uniformity 
with  the  frequency  of  the  fluctuations  and  with  location  on  the  earth’s  surface  are  not  well 
established.  The  mathematical  quantities  that  provide  a  measure  of  the  uniformity  arc  the 
coherence  and  phase  angle,  where  the  coherence  measures  the  similarity  of  the  waveforms 
recorded  by  two  sensors  and  the  phase  angle  is  a  measure  of  the  difference  between  the  times 
of  arrival  of  the  waveforms  at  the  sensors,  and  measurements  of  the  coherence  and  phase 
angle  need  to  be  made  at  representative  locations  over  the  earth’s  surface  for  a  variety  of 
sensor  spacings  in  the  North-South  and  East- West  directions  and  for  a  variety  of  frequencies. 
With  some  hesitation  I  will  add  that  these  measurements  also  need  to  be  repeated  at  different 
times  during  the  day  (to  detect  diurnal  variations),  at  different  times  during  the  year  (to 
detect  seasonal  or  annual  variations),  and  at  different  times  during  a  solar  cycle  (for  solar  cycle 
changes).  I  recommend  the  study  by  Zelwer  and  Morrison  [1972]  as  a  model  for  the  more 
extended  gradiometer-oriented  study,  and  the  report  by  Banerjee  et  al.  [1982]  also  can  be 
recommended  for  its  particularly  relevant  discussion  of  the  coherence  of  geomagnetic  variations. 

7.1.7  Rotating  Gradiometers 

As  we  saw  in  the  preceding  chapter,  the  rotating  induction  loop  gradiometer,  or  RILG, 
has  a  response  to  magnetic  dipole  fields  that  varies  either  with  the  inverse  fourth  power  of 
distance  from  the  dipole  or  as  the  inverse  fifth  power  of  distance,  depending  on  the  way  the 
two  induction  loops  are  connected.  This  variability  of  response  is  of  some  interest  by  itself; 
when  it  is  also  recognized  that  the  response  varies  linearly  with  the  frequency  of  rotation,  the 
RILG  begins  to  look  as  if  it  might  have  some  potential  for  gradient  measurements.  It  is  unlikely 
that  its  sensitivity  will  approach  that  of  superconducting  gradiometers,  but  it  is  the  RILG’s 
other  properties  that  are  distinctive  and  which  need  to  be  studied  further  to  see  if  it  has  any 
unique  application. 

7.1.8  Reexamination  of  the  Source-Free  Assumption 

I  have  discussed  the  source-free  assumption  at  some  length  in  Appendix  A.  The  problems 
that  arise  when  the  assumption  is  made  for  gradient  measurements  in  sea  water  or  other 
conducting  media  on  the  earth’s  surface  need  to  be  studied  further  in  their  appropriate  context. 
For  example,  how  do  the  problems  affect  the  “tracking  and  classification”  function  of  a  five-axis 
superconducting  gradiometer  when  it  is  immersed  in  sea  water? 


7.2.  ('onriiision 


7/ 


§7.2  Conclusion 

To  conclude  this  report,  I  would  like  to  return  briefly  to  my  comments  in  the  introduct  ion 
about  the  great  advances  that  have  been  made  in  magnetic  field  gradiometer  technology  over  t  he 
last  two  decades.  I  hope  I  have  clearly  explained  in  the  body  of  the  report  how  these  advances 
came  about  and  how  I  believe  the  momentum  can  l>e  maintained.  Only  a  small  number  of 
scientists  and  engineers  were  involved  in  the  advances  and  even  today  there  is  every  reason  to 
expc<‘t  that  the  next  improvements  in  gradiometer  |)erformance  will  come  about  as  the  result 
of  the  work  of  a  single  gifted  student  or  professional  researcher.  Unfortunately,  gradiometer 
research  is  difficult  to  conduct  without  gradiometers  and  the  most  sensitive  gradiometers  are 
costly  to  acquire  (particularly  multi-axis  superconducting  gradiometers).  However,  1  have 
r(*peatedly  emphasized  how  gradiometers  and  gradiometer  research  impact  upon  three  areas  of 
crucial  importance  to  the  U.  S.  and  Western  nations  generally  (energy  resource  develo|)ment . 
(h'fense.  and  health  care)  it  is  to  be  ho()ed  that  this  fact  will  become  recognized  at  what  I 
have  termed  the  “research  management”  level  in  the  U.  S.  so  that  funds  will  be  made  available 
to  our  contract  monitors  and  grant  administrators  to  assist  in  gradiometer  purchases.  We  can 
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Appendix  A. 

Validity  of  the  Source-Free  Assumption 


The  Maxwell  equation  linking  the  magnetic  field  at  a  point  with  the  conduction  currcnl 
density  J  and  displacement  current  density  dD/dt  at  that  same  point  is 


curlB  =  //o(J  +  ^)- 
ot 


(Al) 


( 


where  it  is  assumed  that  the  medium  is  non-magnetic  (/x  =  fio).  As  discussed  in  Section  2.3. 
in  a  region  where  J  =  dD/dt  ~  0  (described  as  ‘source  free’),  this  equation  can  be  writ  ten 


curl  B  =  0, 


(A.2) 


and  three  terms  in  the  gradient  matrix  become  redundant.  Under  such  circumstances  only  five 
measurements  of  spatial  variations  of  magnetic  field  components  need  be  made  to  determine  the 
gradient,  of  the  magnetic  field  completely.  This  apparently  purely  academic  result  is  import  ant 
in  magnetic  field  gradiometer  studies  because  its  assumption  forms  part  of  the  theoretical  basis 
for  a  method  of  measuring  the  position  and  moment  vectors  of  a  magnetic  dipole  source  [  e.g., 
Wynn  et  ai,  1975]. 

Unfortunately,  the  source-free  assumption  (i.e.,  the  assumption  that  J  =  dH/dt  =  0), 
involves  an  idealization  that  may  fail  in  practice.  Two  specific  conditions  for  this  failure  to 
occur  are  (I)  use  of  a  highly  sensitive  gradiometer,  which  has  the  capability  of  measuring 
very  small  spatial  rates  of  change  of  the  magnetic  field  components,  and  (2)  operation  of  (  his 
gradiometer  in  a  conducting  medium  (such  as  the  sea)  at,  large  distances  from  the  source  of 
magnetic  field,  i.e.,  at  locations  where  the  gradients  of  the  magnetic  field  are  small.  These 
conditions  are  most  likely  to  apply  when  a  sensitive  gradiometer  is  being  used  in  the  sea 
to  measure  the  nearly  undetectable  gradient  produced  by  a  far  distant  source.  Under  less 
demanding  conditions,  the  source-free  a-ssumption  may  not  necessarily  fail  but  it  will  be  a  poor 
assumptioii  and  measurements  based  on  it  will  be  inaccurate.  The  purpose  of  this  appendix  is 
to  show  how  the  problem  with  the  source-free  a.ssumption  arises  and  to  indicate  some  general 
theoretical  criteria  for  establishing  its  validity. 
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Before  we  examine  the  conditions  under  which  a  region  can  be  considered  to  l)e  source 
free,  there  is  an  important  clarification  that  needs  to  be  made  concerning  terminology.  As  1 
have  stressed  in  the  first  two  lines  of  this  appendix,  Equation  A.l  links  the  magnetic  field  at  a 
point  with  quantities  at  that  point  representing  sources  of  magnetic  fields.  However,  the  .source 
quantities  should  not  be  viewed  as  being  the  source  of  the  magnetic  field  at  the  j)oint,  for 
otherwise  we  would  undoubtedly  have  B  =  0,  and  not  curlB  =  0,  whenever  we  specify  J  = 
dD/dt  =  0.  Thus  we  can  have  a  region,  described  as  source  free,  in  which  J  =  dD/dt  =  0  at 
every  point,  but  where  electric  and  magnetic  fields  exist  that  are  produced  by  sources  outside 
the  region.  Note  that  no  statement  is  made  about  the  size  of  the  region;  it  can  be  of  any 
■size.  It  is  only  in  such  a  region  that  we  can  write  curlB  =  0  and,  as  we  have  seen,  it  is  only 
for  magnetic  fields  in  such  a  region  that  the  gradient  matrix  contains  only  five  independent 
terms.  The  difficulty  that  arises  with  these  source-fre<>  regions  is  a  practical  one.  It  is  simply 
not  feasible,  except  under  the  most  exceptional  circumstances,  to  monitor  J  and  dDjdt  with 
great  accuracy  throughout  a  nominally  source-free  region  where  gradient  rneasurenu'tits  arc- 
being  made  and  to  determine  that  the  two  source  quantities  are  zero  at  all  times.  Whether  it  is 
because  of  this  difficulty,  or  because  of  the  perception  that  J  and  dD/dt  are  negligible  under  all 
conditions  at  low  frequencies,  it  is  not  common  practice  for  the  source  terms  to  be  measured, 
and  they  are  usually  assumed  to  be  zero.  As  I  will  now  show,  this  .source-free  assumption  may 
not  be  valid  when  small  gradients  are  to  be  measured  in  such  a  region. 

At  low  frequencies  in  a  conducting  medium  such  as  sea  water  it  is  common  practice  to 
neglect  the  displacement  current  term  in  Equation  A.l  in  comparison  with  the  conduction 
curti..i,  term.  This  gives  us  a  convenient  point  at  which  to  start  our  review  of  the  .source-free- 
assumption.  If  the  time  variations  of  the  fields  and  source  terms  are  all  assumed  to  be  harmonic 
with  angular  frequency  uj,  Equation  A.l  can  be  written 

curl  B  =  /Lio(J  -I-  z'wD) 

=  //(((ctE  iujeE),  (A..‘i) 

where  c  is  the  permittivity  of  the  medium  and  cr  its  conductivity.  It  follows  that  the  displace¬ 
ment  current  term  in  Equation  A. 3  can  be  neglected  if  a/ux  ^  1.  For  sea  water,  which  has 
a  conductivity  close  to  4  S/m  and  a  permittivity  of  SOcq,  this  condition  implies  that  the  di.s- 
placement  current  can  be  neglected  in  comparison  with  the  conduction  current  for  frequencies 
very  much  less  than  890  MHz.  For  other  common  conducting  media  on  the  earth's  surface,  the 
conductivity  is  generally  substantially  smaller  than  the  conduct  ivity  for  sea  water,  but  typically 
no  less  than  10“*S/m,  and  the  permittivities  of  these  media  will  probably  be  1  2  orders  of 
magnitude  smaller  than  the  exceptionally  high  p<-rmittivity  for  .sea  water.  The  comparison 
frequency  for  these  media  is  therefore  likely  to  Im-  1  2  orders  of  magnitude  smaller  than  890 
MHz.  It  is  clear,  however,  that  if  the  frequencies  of  interest  are  restricted  to  the  ultra-low  and 
extremely-low  frequency  bands  (i.e.,  to  frequencies  less  than  about  3  KHz)  the  displacement 
current  term  will  be  negligible  in  compari.son  with  the  conduction  current  term. 
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The  problem  in  obtaining  a  source-free  region  for  the  gradiometer  measurements  can  now 
be  discussed.  Suppose  we  take  a  relevant  situation  and  assume  that  magnetic  field  gradient 
measurements  are  to  be  made  in  the  sea  at  frequencies  less  that  3  kHz.  There  is  no  doubt 
that  the  displacement  current  term  can  be  neglected  for  this  situation,  but  the  conduction 
current  term  remains  and  it  may  be  substantial,  depending  on  the  frequency.  The  currents 
included  in  J  may  be  induced  by  acoustic  or  infrasonic  waves  propagating  through  the  sea,  by 
surface  and  internal  waves  or  other  sea  motions,  or  by  magnetic  field  fluctuations  originating 
above  the  sea  surface  (e.g.,  by  geomagnetic  pulsations),  to  name  just  three  general  clas.ses  of 
possible  sources.  Recognizing  the  problem  obtaining  J  =  0  in  the  sea,  or  other  conducting 
medium  on  the  earth’s  surface,  suppose  we  move  our  gradiometer  up  above  the  sea  and  proceed 
to  make  gradient  measurements  in  the  air,  which  is  largely  non-conducting.  The  conduction 
current  term  can  now  be  neglected,  but  the  displacement  current  term,  although  small,  cannot 
be  discarded  automatically,  since  it  is  the  only  .source  term  remaining.  For  this  situation. 
Equation  A.3  has  to  be  reexamined  to  see  when  the  displacement  current  can  be  neglected. 

Expanding  Equation  A.3,  with  =  0,  we  have 


rdB,  _ 

V  dy  dz 


\  dBA./dBy  dB^\  . 

j  ■  + (.37  ■  ^  j-' +  laT  ■  a;- J 


Thus,  for  the  displacement  current  to  be  completely  negligible,  the  following  three  condit  ions 
should  all  apply: 

dBy  _ 

»  (jJHQiEx, 


dBz 

dx 


»  IjJHOtEy, 


(A.4) 


dBx 

dy 


Note  that  simultaneous  measurements  of  the  electric  field  components  and  of  three  spatial  rates 
of  change  of  the  magnetic  field  components  are  required  to  determine  the  applicability  of  these 
conditions. 


How  important  are  the  above  conditions?  Few  electric  field  data  are  available,  but  the 
limited  measurements  that  have  been  made  indicate  that  the  displacement  current  conditions 
are  easily  satisfied,  whereas  the  equivalent  conduction  current  conditions  are  not  as  easily 
satisfied.  For  example,  suppose  the  magnetic  field  gradient  is  to  be  measured  in  a  ‘source-free 
region’  at  a  frequency  in  the  Pc  1  pulsation  frequency  range  (0.2  5  Hz).  For  convenience,  I  will 
assume  that  the  frequency  is  1  Hz.  Troitskaya,  [1961]  quotes  amplitudes  of  0.01  0.1  mV/km 
for  the  horizontal  electric  fields  induced  in  the  earth  by  these  pulsations  during  quiet,  days  at 
middle  latitudes,  and  amplitudes  in  the  tens  of  mV/km  on  disturbed  days.  These  electric  fields 
are  tangential  to  the  earth/air  boundary  and  they  should  therefore  apply  in  the  air  just,  above 
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the  earth’s  surface.  If  we  suppose  it  is  Ex  that  is  being  measured  and  that  there  is  a  mild 
Pc  1  disturbance  in  progress,  the  maximum  value  of  Ex  can  be  taken  to  be  lOmV/krn  or  10 
/tV/m,  giving  dByjdz  =  ujfioeEx  =  ijJ^iQtoEx  —  7.0  X  10“*°  pT/m.  This  spatial  variation  of 
magnetic  field  is  much  lower  than  the  noise  level  of  0.1  0.3  pT/m  (1  Hz  bandwidth)  that  ha.s 
been  measured  at  1  Hz  with  a  sensitive  single-axis  superconducting  gradiometer  [Gillespie  et 
al,  1977].  This  calculation  suggests  that  the  displacement  current  source  term  can  certainly 
be  ignored  in  the  air.  Since  the  conduction  current  term  is  also  very  small,  we  can  conclude 
that  the  source-free  approximation  is  valid  in  the  air.  The  situation  is  quite  different  in  the  sea 
or  in  the  earth. 

For  conducting  media,  the  source-free  criteria  equivalent  to  those  above  for  a  non-conduct¬ 
ing  medium  are 

>  PooEx, 

BTi 

>  HoaEy,  (AJi) 


where  the  displacement  current  density  has  been  neglected.  Once  again  using  Troitskaya’s 
electric  field  amplitude  of  10  /xV/m  for  the  Pc  1  pulsations  occurring  during  a  day  of  mild 
activity,  an  amplitude  that  is  likely  to  be  equally  applicable  in  the  sea  or  in  the  earth, 
we  obtain  values  for  pocrEx  of  50  pT/m  for  sea  water  (cr  =  4.0  S/m)  and  0.13  pT/m  for 
moderately  dry  earth  materials  (<7  =  10“^  S/m).  The  former  of  these  spatial  variations  would 
be  easily  measurable  with  a  superconducting  gradiometer,  and  the  latter  probably  also  could 
be  measured  at  1  Hz  [Gillespie  et  al.,  1977].  The  sea  water  spatial  variation  is  quite  large, 
and  if  the  computation  were  to  include  the  electric  fields  generated  by  motions  of  the  water 
in  the  presence  of  the  earth’s  steady  magnetic  field  it  could  well  be  considerably  larger.  These 
computations  aiC  order-of-magnitude  only,  and  they  only  apply  over  a  limited  frequency  range, 
but  similar  results  would  probably  be  obtained  for  other  frequencies  in  t  he  ultra-low  range  (f 
<  5  Hz). 


Up  until  this  point  we  have  been  considering  low-frequency  magnetic  (and  electric)  fields 
and  sources.  Let  us  now  briefly  consider  the  situation  at  dc.  The  displacement  current  term 
dD/dt  is  necessarily  zero  for  /  =  0  and  it  is  probable  that  the  conduction  current  term  J 
will  be  close  to  zero,  since  there  are  no  obvious  mechanisms  for  maintaining  steady  currents  in 
the  conducting  media  encountered  on  the  earth’s  surface.  Under  these  circumstances,  it  might 
well  be  thought  advantageous  to  use  dc  measurements  of  the  various  terms  occurring  in  the 
magnetic  field  gradient;  it  certainly  appears  that  these  measurements  could  be  made  in  regions 
satisfying  the  source-free  assumption  and  that  all  five  independent  terms  in  the  gradient  matrix 
could  be  measured  satisfactorily.  There  are  .some  difficulties,  however.  The  designation  dc  is 
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to  some  extent  relative  in  practice,  and  there  exist  a  variety  of  slowly  varying  electromagnetic 
phenomena  in  the  geophysical  environment  that  would  appear  to  be  dc  during  the  courst'  of 
a  ‘dc’  measurement.  In  particular,  there  are  well-known  slowly  varying  ionospheric  current 
systems  that  produce  electromagnetic  fields  on  the  earth’s  surface  and  below  (for  example,  the 
electric  current  system  producing  the  quiet-day  solar  daily  variation  [Chapman  and  Bartels, 
1940]),  and  these  electromagnetic  fields  will  appear  to  be  dc  over  intervals  of  several  minutes 
up  to  several  tens  of  minutes.  For  measurements  in  the  air,  the  fields  can  be  ignored,  but  in 
the  sea  or  other  conducting  media  on  the  earth’s  surface  there  is  likely  to  be  a  significant  J 
and  the  source-free  assumption  will  not  be  applicable. 

Two  new  circumstances  also  arise  at  dc.  First,  some  sources  of  interest  are  variable  they 
are  either  rarely  at  rest  or  they  have  varying  magnetic  dipole  moments.  Unless  these  variations 
take  place  very  slowly,  a  dc  approach  is  not  applicable.  The  second  circumstance  is  likely  to  be 
more  important:  at  dc  the  earth’s  own  magnetic  field  becomes  a  significant  factor.  As  shown  in 
Appendix  B,  the  earth’s  steady  magnetic  field  has  a  gradient  and  the  spatial  variations  included 
in  the  gradient  are  substantial.  Whether  these  circumstances  are  relevant  in  practice  or  not 
will  depend  on  the  particular  situation. 

The  limits  on  the  applicability  of  the  source-free  assumption  are  now  reasonably  clear. 
There  certainly  seems  little  doubt  that  the  assumption  is  valid  in  the  air,  or  other  non¬ 
conducting  medium.  However,  in  the  sea,  and  other  conducting  media  typically  met  on  the 
earth’s  surface,  it  appears  that  the  source-free  assumption  will  not  generally  be  valid,  although 
there  may  be  certain  frequency  ranges  (yet  to  be  determined)  where  the  assumption  is  not 
a  bad  approximation.  There  is  another  problem  with  conducting  media  that  should  also 
be  mentioned.  It  is  obviously  not  good  experimental  practice  to  immerse  a  magnetometer 
or  gradiometer  directly  in  the  sea,  for  example.  Thus  some  form  of  container  is  required. 
Superconducting  instruments  already  have  their  sensing  assembly  enclosed  in  a  dewar,  but 
the  dewars  are  not  usually  designed  for  immersion  in  the  sea  and  an  additional  container  is 
necessary.  These  containers  change  the  distribution  of  current  density  around  the  .sensors 
and  reduce  it  to  zero  in  their  immediate  vicinity,  thus  distorting  the  measurements  a  typical 
example  of  the  measuring  instrument  altering  the  quantity  it  is  supposed  to  measure  during 
the  measurement  process.  This  problem  pertains  to  any  gradiometer  irumersed  in  the  sea 
in  a  container;  an  additional  problem  arises  if  the  container  is  in  motion,  for  example,  if  it 
is  being  towed  behind  a  vessel.  In  this  latter  case  the  turbulent  flow  of  the  sea  around  the 
container  in  the  presence  of  the  earth’s  magnetic  field  will  generate  electric  currents  that  will 
produce  spurious  signals  in  the  gradiometer.  Joseph  et  al.  [1981]  and  Joseph  and  Thomas  [1981] 
have  studied  these  container  problems  and  proposed  corrections  for  motion  of  the  container 
through  the  sea.  However,  the  distortions  produced  simply  by  the  presence  of  the  container, 
without  motion,  cannot  be  avoided  and  they  must  be  taken  into  account  when  interpreting  the 
gradiometer  measurements. 

In  summary,  it  appears  that  the  source-free  assumption  will  be  valid  for  gradiometers 
operating  in  the  air  at  frequencies  less  than  3  kHz.  However,  for  gradiometers  operating  in 
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such  conducting  media  as  the  sea  or  the  earth  the  source-free  assumption  is  unlikely  ever  to  be 
entirely  valid  in  practice  and  under  some  circumstances  it  could  be  substantially  incorrect.  The 
containers  necessary  for  operation  in  tha  conducting  media  also  introduce  perturbations  into 
the  measurements.  The  implications  of  these  problems  for  gradient  measurements  in  the  sea 
or  earth  are  unfavorable,  and  it  appears  unlikely,  in  particular,  that  the  technique  proposed  by 
Wynn  et  al.  [1975]  for  deriving  the  location  and  strength  of  a  magnetic  dipole  source  will  be 
effective  in  the  sea  except  at  very  short  ranges.  I  might  conclude  by  observing  that  experimental 
data  are  badly  lacking  in  this  area,  and  the  right  experimental  results  could  alter  the  negat  ive 
picture  I  have  drawn  for  measurements  in  conducting  media. 


Appendix  B. 
Gradient  of  the  Geomagnetic  Field 


The  earth’s  steady  magnetic  field  approximates  to  a  dipole  field,  with  the  dipole  located 
at  the  earth’s  center  and  tilted  at  an  angle  of  about  11.4°  to  the  rotation  axis.  The  moment  of 
this  dipole  is  about  0.815  X  10^^  A  m^,  and  the  field  calculated  from  it  is  commonly  referred 
to  as  the  geomagnetic  field.  We  can  use  this  dipole  field  model  to  estimate  the  gradient  of  the 
earth’s  magnetic  field,  provided  it  is  understood  that  is  is  an  imperfect  model  on  the  earth’s 
surface  and  that  there  will  be  regional  and  local  anomalies  where  the  gradient  is  likely  to  Ix' 
larger  than  our  estimates. 

Using  a  spherical  polar  coordinate  system  based  on  the  geomagnetic  dipole,  as  shown  in 
P’igure  B.l,  the  geomagnetic  field  can  be  written 

1)^ 

- ■ 

47rr'’ 


iBE)0 


fioMsinO 

4wr^ 


where  M  =  0.815  X  10^^  A  m^  is  the  dipole  moment,  (BB)r  is  the  radial  component  of  the  field, 
and  {Be)0  is  the  angular  component  (there  is  no  azimuthal  component).  Since  MKS/SI  units 
are  used  in  this  report,  the  units  of  the  magnetic  field  quantities  in  Equation  B.l  are  Tesla. 
The  orientation  of  the  components  relative  to  the  dipole  axis  is  shown  in  Figure  B.l;  the  figure 
also  shows  how  the  geomagnetic  dipole  points  toward  the  south. 

As  a  check  of  the  above  expressions  for  the  geomagnetic  field  and  of  the  geomagnetic  dipole 
moment,  consider  the  magnitude  of  the  geomagnetic  field  Bg  at  the  poles  (0  =  OP,  180P). 
Denoting  the  eai  th’s  radius  by  /?£,  we  have 
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which  gives  \^E\poie  —  6.3  X  10“®  T  (i.e.,  0.63  Gauss)  if  Re  is  assumed  to  have  a  value  of 
6370  km.  This  value  for  the  polar  geomagnetic  field  is  consistent  with  the  total  magnetic  field 
amplitudes  usually  shown  for  the  polar  regions  [e.g.,  Dessler,  1965]. 

Because  the  geomagnetic  field  is  a  vector  field  it  would  normally  be  expected  to  have  t  hree 
independent  parts  to  its  gradient.  However,  because  the  field  is  symmetric  about  the  dipole  axis, 
it  only  has  two  parts,  grad  {BE)r  and  grad  {BE)e,  in  our  adopted  spherical  polar  coordinate 
system.  These  gradients  can  be  written 


grad  {BE)r  = 


d{BE)r  .  1  d{BE)r 

dr  ^  r  dO 


9, 


grad  {Be)b=^ 


d(BE)e  .  1  d(BE)6  x 

dr  ^  r  86 


(B.li) 


where  r  and  6  are  unit  vectors  in  the  direction  of  increasing  r  and  6.  We  therefore  have  four 
gradient  terms  to  consider. 


Using  equations  B.  1,  the  four  terms  can  be  expressed  in  the  following  functional  forms: 

d{BE)r  —dfJnM cosO 

dr  iirr'^  ’ 
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d{BE)$ 
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Table  B.l  gives  representative  values  of  these  spatial  variation  terms  on  the  earth’s  surface 
{r  =  Re  =  6370  km)  for  0  =  OP,  180P  (polar),  45°  and  135°  (midlatitudes),  and  90P  (equatorial). 
These  values  show  that  the  steady  magnetic  field  has  substantial  gradients  even  in  the  absence 
of  the  regional  and  local  anomalies  that  are  known  to  exist. 


First  consider  the  polar  fields  at  the  earth’s  surface.  Combining  the  previous  equ.ations 
and  the  data  in  Table  B.l,  these  fields  have  the  following  properties: 


{BE)r  = 


2pioM 

AttRe^ 


=  6.3  X  10“®  T, 


{Be)«  =  0, 


grad  {BE)r 
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(B.5) 


grad  {Be)0  =  g  =  1^5.0  »  pT/m, 
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Figure  B.l.  Geometry  of  the  geomagnetic  field.  Note  the  orientation  of  the  geomagnetic 
dipole;  the  geomagnetic  field  is  directed  out  of  the  earth  at  the  south  geomagnetic  pole 
(r  ==  Re,  0  =  OP)  and  into  the  earth  at  the  north  geomagnetic  pole  (r  =  Re,  0  =  180P). 


where  the  upper  and  lower  signs  in  the  gradient  terms  apply  to  the  north  and  south  polar  fields 
respectively.  Because  {BE)e  —  0.  the  polar  magnetic  fields  are  entirely  radial,  or,  equivalently, 
vertical  to  the  earth’s  surface.  The  gradient  of  the  vertical  component  of  the  geomagnetic  field 
is  also  entirely  radial  at  the  poles,  but  its  direction  is  antiparallel  to  the  direction  of  the  field. 
Thus,  at  the  north  geomagnetic  pole  the  field  is  directed  into  the  earth’s  surface,  whereas  the 
gradient  of  the  vertical  component  (the  only  component  of  the  field  at  the  pole)  is  directed 
outwards.  Note  that  the  gradient  of  the  horizontal  component  of  the  field  is  not  zero  at  (he 
poles,  even  though  the  component  itself  is  zero.  Howt  ver,  the  magnitude  of  the  gradient  of  t  he 
horizontal  component  (5  pT/m)  is  substantially  smal'er  than  the  magnitude  of  the  gradient  of 
the  vertical  component  (29.7  pT/m). 
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Appendix  B.  Gradient  of  the  (Jeoinagnetir  Field 


TABLE  B.l.  Some  representative  numerical  values  of  the  four  terms 
comprising  the  gradient  of  the  geomagnetic  field  in  the  geomagnetic- 
dipole-based  spherical  polar  coordinate  system  (Figure  B.l).  The  units 
are  pT/m,  and  the  values  apply  on  the  earth’s  surface  (r  =  6370  km)  at 
the  geomagnetic  south  pole  {Q  —  OP),  southern  midlatitudes  (0  =  45°), 
the  geomagnetic  equator  (0  =  90P),  northern  midlatitudes  (0  =  135°), 
and  the  north  geomagnetic  pole  (0  =  180P). 


0 

d{BE)r 

dr 

1  a(BE)r 
r  dO 

d{BE)e 

dr 

1  d{BE)e 

r  dO 

OP 

-29.7 

0 

0 

+5.0 

45° 

-21.1 

-7.0 

-10.5 

+3.5 

90P 

0 

-9.9 

-14.8 

0 

135° 

+21.1 

-t-7.0 

+  10.5 

-3.5 

ISflP 

-t-29.7 

0 

0 

-5.0 

The  equatorial  field  at  the  earth’s  surface  can  also  be  obtained  from  the  preceding  equations 
and  the  data  in  Table  B.l;  it  has  the  following  properties: 


(BE)r  =  0, 


iBE)e  =  =  3.4  X  10-'"  T, 

AttRe^ 

grad  (Be )r  =  6  =  -9.9  0  pT/m, 

AttRe 

grad  (Be)$  =  f  =  —14.8  r  pT/m. 

AitRe 


(B.6) 


Because  {BE)r  =  0,  the  equatorial  geomagnetic  field  is  entirely  horizontal.  The  gradient  of 
the  radial  component  of  the  field  is  also  horizontal,  even  though  the  radial  component  it.s<4f 
is  zero.  The  gradient  of  the  angular  component  (Bg)®  is  radial  (i.e.,  vertical  to  the  surface) 
and  directed  into  the  earth;  its  magnitude  of  14.8  pT/m  is  larger  than  the  magnitude  of  9.9 
pT/n(  for  the  gradient  of  the  radial  component,  but  the  difference  between  the  two  magnitudes 
is  much  less  than  the  corresponding  difference  for  the  polar  gradients. 
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